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A GENERAL SOLUTION FOR THE RECTANGULAR AIRFOIL 
IN SUPERSONIC FLOW* 


BY 
JOHN W. MILES** 
University of California at Los Angeles 


Summary. The potential on a rectangular airfoil due to an arbitrarily prescribed 
motion at its surface is obtained by an operational solution of the linearized equations 
and subsequent comparison with the known solution in steady flow. It is shown that 
the result can be extended to more general planforms with the aid of the Lorentz trans- 
formation. Other methods of solution are noted. 

1. Introduction. The problem of unsteady motion of a rectangular airfoil in super- 
sonic flow has been treated in closed form by Goodman [1, 2], Miles [3, 4], Rott [5], 
Stewartson [6], and Stewart and Li [7, 8]. The results of refs. 1-6 are in mutual agree- 
ment, but those of Stewart and Li are believed to be incorrect due to their (in our 
opinion?) erroneous conclusion that Evvard’s “equivalent area’’ concept is applicable 
to non-stationary flow, the derivation [9] of which has been criticized elsewhere [10]. 

The reason for the addition of the present paper to this already voluminous literature 
is to present a solution that is valid for an arbitrarily prescribed motion of the airfoil. 
In principle, the solution of ref. 3 is sufficiently general in virtue of Fourier’s theorem 
(with reference to the time dependence) and the possibility of expanding velocity dis- 
tributions of practical interest in powers of y", but the following solution has the ad- 
vantage of exhibiting the dependence of the potential on the prescribed velocity in a 
more explicit form. Moreover, the end result affords an immediate extension to plan- 
forms with oblique edges. Finally, the result permits a precise extension of the equivalent 
area concept to non-stationary motion of a rectangular wing, albeit this interpretation 
does not permit further generalization to other planforms. 

2. Formulation of the problem. We consider a quar.er infinite airfoil whose pro- 
jection on the plane z = 0 is bounded by X = —Ut/l and y > o in the fized, right 
handed, dimensionless (with / as the characteristic length), Cartesian coordinates 
(X, y, 2), U being the flight velocity and ¢ the time. On the basis of the usual assumptions, 
the linearized equation for the velocity potential is the wave equation 


xx + Puy + Pee = Orr (2.1) 
T =at/l (2.2) 


*Received December 26, 1951. 
**This work was carried out while at Auckland University College, New Zealand under the Ful- 
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where a is the sonic velocity and 7 a dimensionless time. In the end results we shall 
exhibit ¢ as a function of the moving coordinates (x, y, z), related to the fixed coordinates 
(X, y, z) by the Gallilean transformation 


x 1 M X 
- (2.3) 
t 0 1 IT’ /a 
but in carrying out the analysis it is expedient to introduce the modified Lorentz trans- 
formation 
g ] M\/X 
=p" , b=(M-1)" (2.4) 
T M ] T 
WE, ¥, 2, 7) = $2, ¥, 2, DY (2.5) 


Transforming the original wave equation (1), we obtain the new wave equation (with 
the positions of the axial and time variables reversed) 


Vere + Vos a Wes sed Vee (2.6) 
which is conveniently regarded as the canonical form for non-stationary flow at super- 
sonic speeds. [If in (4) we take 8 = (1 — M’)'”, M < 1, we have the classical Lorentz 


transformation, and (1) remains invariant with respect to the position of the space 
and time coordinates. The application to subsonic airfoil theory has been discussed by 


Kussner [11].] 
The boundary condition to be invoked at the airfoil is 


d, = —w(z, y, b); z>0,y>0,z2z=0 (2.7) 
where w is the prescribed velocity, positive down. In the transformed coordinates, we 
write 

y, = —v(é, y, 7); é>0,y>0,z=0 (2.8) 

vf, Y, 1) = w(z, y¥,) = w([Bé, Y; (1/Ba)(Mé aay 7)] (2.9) 


In consequence of the symmetry of this boundary condition with respect to z = 0, the 
potential itself is anti-symmetric thereto, and we may restrict our consideration to the 
half space z > 0, provided that we pose the additional boundary condition 


y=0; §>0,y<0,z2=0 (2.10) 


In addition to the local boundary conditions (8) and (10), it is necessary to impose 
the requirement that y behave asymptotically as a disturbance originating at the wing 
and vanishing identically forward of the zone of action of the wing, toward which end 


it suffices to write 
yv=0; &<0 (2.11) 


together with the Sommerfeld finiteness and radiation conditions in the original co- 
ordinates (X, y, z, JT), although the latter may not be necessary in all cases, 
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To complete the formulation, we note that the perturbation pressure at the upper 
surface of the wing is given by 


Pp = —pol(d.)x = —pas'(My+¥,); §€>0,y>0,2 = 0+ (2.12) 


We now pose the following problem: given w, find ¢ on z = 0+ or, equivalently, 
find a solution to (6) satisfying the boundary conditions (8), (10) and (11), together 
with the appropriate conditions at infinity. 

3. Reduction to steady flow problem. Let W be the transform obtained by posing the 
time dependence exp (ixr) and taking a Laplace transform with respect to &, viz. 


exp (ix7)W(s, y, 2, x) = &fy} = [ e WE, y, 2, 7) dé (3.1) 


with a similar representative, V, for v. Transforming the boundary value problem of 
the preceding section, we have 


Vv, +, —’v =0, V=s+ (3.2) 
v¥,=-V; y>0,z=0 (3.3) 
Vv =0; y¥<0,2=0 (3.4) 


We now consider the “quasi-steady” problem obtained by setting A = s (x = 0) 
in (2) [but « # 0 in (3)], which is tantamount to neglecting y,, in (2.6), thereby reducing 
the boundary value problem to one in steady flow. The solution to this reduced problem, 
which is obtained most directly by Evvard’s method [12], is designated by ¥ and given 
at the wing (all subsequent potentials also are specified at z = 0+-) by 


E-u) + 


a ew 2 2 -1W2 
vou 0t,)=2f def (e-w-Y- MTG, dn 5) 


1(&—m) 
the integration being carried out over the trapezoid bounded by the Mach lines passing 
through (£, y), the reflection of one of these [7 = y — (€ — u)] in the side edge, and the 
leading edge. 

Taking the Laplace transform of ¥ , we have, by a slight extension of the Faltung 
heorem, 


Ty tg ” 
Vo = Le{Yo} = eed! [ : [? — (y — 0)*)'v(u, 2, 7) an} (3.6) 
“li—s9 
Carrying out the indicated operation on v and substituting \ for s in the operation 
with respect to &, we have 


V(s, y, O+, x) 


© 1 E+y (3.7) 
= [ exp [—(s? + amet [. | [ — yy — »)*J'? VG, 2, an} dt 


where the phase of (s? + «’)'” is chosen to ensure a positive real part. 
4. Solution for harmonic time dependence. To effect the inverse Laplace transform 


of (3.7) we introduce the theorem [cf. p. 123 of ref. 13, a source hereafter denoted by MO] 


*/2 
SFE +e =fO+ [feos F Jueesin de (4.1) 
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the application of which yields 


4 ne oe ‘ ; 
y=Yt+ - | dé I 06 Jolx(E — w) sin 6] dy 


. [(é aap u)” cos" i= (y ~~ n) |] “U(u, n, 7) dy 


os8—y! 
u)cos 7, 


where the second term appears as a correction on the ‘“‘quasi-steady”’ solution. 
An alternative expression for y that separates out the “strip theory” contribution 


y, = | Jolx(E — uw) ]v(u, y, 7) du (4.3) 
can be obtained from (2) after integrating by parts and noting that 
lim | ((é — un) cos’ 6 — (y — n)*J] n(n, 0, 7) dn = mE, y, 7) (4.4) 
6-4/2 « t—n)cosd—1 
whence 
, OF ec a . 
Y,=yY¥,-—- | ag | Jo[x(—é — pw) sin 6] du 
(4.5) 
o—# J ry 
0 > ? 2 2 
| [((é — uw)” cos 6 — (y — n)] “v(m, 2, 7) dn 


o6. 


u)cost—y 


In the special case where v is assumed to be independent of y the » integration may 
be carried out explicitly, the result being 
ei i Qo. 
y=y,-- | do | Jo[x(é — uw) sin O]v(y, 7) 5p Sin” [y/(E — pw) cos 6]'”? du (4.6) 
where the arc sine is to be replaced by 7/2 when its argument exceeds unity. Changing 
the variable of integration yields 


y=%- : } du Jo[x(u* — y° sec* yg)" w(E — wu, 7) de (4.7) 
The last result is essentially in the form given by Stewartson [6]. [We remark that it 
was the form of Stewartson’s result that suggested the present approach, although the 
use of (8) in connection with supersonic airfoil theory is due to Magnaradze [14, 15] 
and has been applied previously to the rectangular wing by Galin [14], but without 
much progress, since he found it necessary to introduce Fourier series in & (cf. ref. 15).] 
5. Arbitrary time dependence. (Generalizing the result (4.2) we have, after Fourier 
transformation (MO119) and convolution, 


© af a(E—p)cosd+y 


dé | du 


/ 1 (§—p) cosO—y 


; _ . = 
Y= vt | [( — u)* cos’ 6 — (y — »)*}'” dy 


(5.1) 
a at+(E&—p) sind 
= | [((€ — uw)’ sin’ 6 — (r — 3°", 0, 9) a 
00 / r—(E—p) sind 
In connection with the ¢ limits of integration, we remark that the maximum and mini- 
mum values of ¢ permitted in v are ~ and r — M(é — y), corresponding to — ~ and 
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present time after transformation back to the physical variables. Both of these limits 
being outside of the domain in which the inverse transform of J» [x( — yw) sin 4], is 
non-vanishing, we choose the limits as shown, although v may vanish over some part 
of the latter region, as in transient problems. 

A more convenient expression for y, obtained by introducing a trigonometric variable 
in place of ¢, is given by 


ax/ a (E—-p) cosO+y 


| J a 3 9 9 27y- 
y=y+ = | dé | dx [ du | [((é — u)’ cos’ 6 — (y — )*)'” 
J /0 /0 “| (&-p)cos0—y! 


(5.2) 
ra] } 
° 39 Ul n, T os (E a ht) sin 6 cos x] dn 


Numerous, alternative forms may be obtained by additional changes of variable and 
by integration by parts. 

Finally, upon substituting the physical variables zx, t, @ and-w from (2.2, 2.3, 2.4, 
2.5, 2.9) we obtain 


az aB-*(r—p) +4 


l ‘ / 2 2 21-1/2 
bo =F J de | l(c — uw) — By — VT” 
r Jo J |B-*(2—-u) -y 
(5.3) 
- wy, n, t — (M1/B*a)(x — w)] dy 
1 r/2 a? ar 8-*(z—-p)cos6+y - . : iain 
dé=a+3 | dé | dx | du [ [(w — »)’ cos’ 6 — By — »)*)'” 
7 ) “0 “0 « |B" (2—p) cosd—y!| 


(5.4) 


< wiu, 7, t — (1/B%a)(M + sin 6 cos x)(x — p)] dn 


6. Oscillating wing. For the important special case of a wing undergoing the oscil- 


lating motion prescribed by 

w(x, y, t) = Riiw(2, ye’*'} (6.1) 
we have only to choose 

x = kM /8, k = wl/U (6.2) 


in (4.3) and (4.5) and return to the original variables, whence 


( 


?,. = B RIN | exp [iwt — i(kM?/s’)(x — w)]Jo[(kKM/B’)(x — p)]w(u, y) du} (6.3) 


7 


do | | exp [iwt — i(kM?/8’)(x — w)]Jo[(kM/8")(a — yw) sin 6] du 


o=>-odo, — 1 pd | 
(6.4) 


a 


- »B-*(z—p) cosO+y 
} 9° 9 9 9 — 9 
= | [((x — u)* cos’ 6 — B(y — n)*)*w(u, 0) dn 
c J 1g-3 


z—y) cos6—y| 


7. Extension of Evvard’s “equivalent area” concept. It is of some interest to note 
that in (4.2) et. seg. the domain of the (u, 7) integration is bounded by the pseudo Mach 
lines n = y + (E — pn) cos 6, together with the reflection of » = y — (€ — u) cos 8, 
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namely 7 = ( — u) cos @ — y, in the side edge. This interpretation furnishes an ex- 
tension of Evvard’s ‘equivalent area” concept to unsteady flow, but the extension is 
not valid for oblique or curved side edges, since, in general, the equation for the reflected 
Mach line would be of the form 7 = f(t, & — uw) — y, rather than simply f(é — ») — y, 
and the Faltung theorem [ef. (3.6)] would not be applicable. 

8. Oblique edges. The extension of the foregoing results to wings having arbitrarily 
prescribed supersonic leading edges and streamwise side edges is trivial, since it is 
necessary only to circumscribe a fictitious rectangular wing and set w = O over those 
portions that are not included in the original planform. [However, it should be remarked 
that the Mach lines from the opposite corners may intersect on the wing only if it is 
rectangular (cf. ref. 3).] To extend the results to a straight, subsonic leading edge that 
is adjacent to a supersonic leading edge and does not interfere with any other subsonic 
edge, we may apply the Lorentz transformation 


, 


é’ ] m 
= (1 — m’)'” , m>l1 (8.1) 
y’ m ] y 


w¥r 


under which (2.6) remains invariant. Unfortunately, it is then no longer possible to 
obtain results for the spanwise integrals of the potential that are comparable in sim- 
plicity to those for the rectangular wing (cf. refs. 3, 5, 6.). 

The’ use of (1) in connection with supersonic wings in steady flow is well known 
(cf. refs. 16-18), but the application to unsteady flow problems seems not to have been 
noticed previously, perhaps due to the fact that the differential equation usually is 
written for ¢(z, y, z, t), rather than y(&, y, z, 7). 

9. Other methods. It is rather curious that of all the methods that have been applied 
to the problem of diffraction by a half plane, the classical counterpart of the supersonic 
rectangular wing, none [e.g., Poincaré’s original treatment, Sommerfeld’s celebrated 
application of multivalued integrals, Lamb’s use of parabolic coordinates, Magnus’ 
solution of the integral equation, the Wiener-Hopf technique applied by Copson and 
Schwinger, etc.; cf. ref. 19] has proved as effective as the methods derived especially 
for the wing problem, notably Busemann’s use of homogeneous solutions (‘“‘conical 
flows’’) and Evvard’s method. (We remark that each of these methods had antecedents 
in the work of Bateman [20].) This is at least partially due to the different foci of interest 
in the two situations, viz. the potential on z = 0 in the wing problem and the solution 
at a distance (particularly near the boundary of the geometric shadow) in the diffraction 
problem, but it is of interest to note that the method of conical flows furnishes an elegant 
approach to the problem of pulse diffraction [21, 22, 23], while the application of Evvard’s 
method to case of an arbitrary incident wave has received attention in a recent paper 
by Friedlander [24]. 

The foregoing remarks notwithstanding, Lamb’s method has proved rather at- 
tractive for the special case of no spanwise variation of w, an application discussed 
recently by Rott [5]. This method is, in fact, applicable to more general spanwise dis- 
tributions [25], but, in the form presented by Rott, it is much less direct than the present 
method, from which it differs fundamentally in prescribing the boundary data on the 
wing tip Mach cone, rather than the wing proper. However, if the potential is desired 
only at the wing the boundary data may be prescribed there. Thus, a solution to (3.2) 
in the polar coordinates (p, ¢) that vanishes on z = 0, y < 0 (g = =), gives a null value 
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of the normal derivative on the wing (g = 0) and vanishes properly at infinity is given by 
Vv, = F(s, x)p'e~™ cos (v/2). (9.1) 


This solution evidently is appropriate to the case where v is independent of y. More- 
over, at large distances from the edge ¥ must reduce to the two dimensional solution, 
viz. [ef. ref. 6 or (3.2) and (3.3) with ¥,, = 0] 


lim v = ¥, = A7'v(s, ye". (9.2) 


yo 


Integrating (1) subject to the conditions that Y must vanish at y = 0 and satisfy (2) 
asymptotically, we have on g = 0 


Vv = d7~'v(s, x) erf [(Ay)'”?] (9.3) 


the inversion of which yields the result (4.7). 
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REFLECTION OF WAVES FROM VARYING MEDIA* 


BY 
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Abstract. Formulae are found for the coefficient of reflection from varying media of 
a type encountered in physics. These are applied approximately for some general classes 
of media, and exactly for some specific cases. Many media which would normally be 
expected to be highly reflecting are shown to be completely transparent to certain 
waves at least and, in some cases, to a whole spectrum of waves. The results are con- 
sidered both for electromagnetic (or other classical) waves and for mass waves. 

I. Introduction and Summary. Many problems of wave propagation through varying 
media introduce the equation 

2 

5+ oa)f = 0 (1) 
for solution. This equation was studied some years ago by Lord Rayleigh [1] and others 
in connection with classical wave problems. The advent of quantum mechanics has led 
to a more thorough investigation of the solutions of (1), centered around the method 
developed by Brillouin [2], Wentzel [3] and Kramers [4]—the B.W.K. method. More 
recent work along these lines has been carried out by Kemble [5, 6], Langer [7, 8, 9], 
and Furry [10]. Exact solutions, in terms of hypergeometric functions, can be found if 
g(x) is of the proper form; this method of attack was developed by Eckart [11] and 
Epstein [12]. In all these treatments, great attention has been paid to solutions obtained 
when g(x) is positive for some ranges of x and negative for others. Such solutions are 
of particular interest in quantum mechanics, in some oblique incidence problems, and 
in problems of electromagnetic wave propagation through a conducting medium such 
as the ionosphere. 

In the present paper we will deal with functions g(x) which are everywhere real, 
finite, and continuous, have a constant positive value for large values of | x |, and lesser 
values in a region of variation about z = 0. 

Solutions are obtained in terms of a new variable r = r(x) which is itself a solution 
of a Rieccati equation.t Exact formulae, also in terms of r, are developed for the re- 
flection coefficient of a medium. The curve of r(x) is traced by inspection of the de- 
termining differential equation. A surprisingly large amount of information about r can 
be found in this way, but the procedure is only qualitative. 

By this method, then, we obtain an approximate evaluation of the exact solution, 
rather than an exact evaluation of an approximate solution as is found in the B.W.K. 
method. This difference of approach leads to a difference in applicability. The B.W.K. 
method is of prime importance in deriving numerical results in a particular given 
problem, although it has been used to infer some general results. The present method 


*Received Sept. 12, 1951. Received in revised form March 14, 1952. 

**Present address: Fitzwilliam House, Cambridge, England. 

+See the paper by S. A. Schelkunoff [16] for comments on the occurrence of the Riccati equation in 
wave transmission problems and, indeed, for a number of interesting and pertinent remarks on the 
whole problem of propagation under varying conditions. 








C. 0. HINES [Vol. XI, No. 1 


10 


leads to the construction of problems to which it can give an exact solution, predicts some 
general results, and indicates some surprising (?) results which can be obtained if the 
medium is of the proper form. In particular, it is shown that many media which might 
be expected, by inference, to be highly reflecting can be completely transparent to one 
frequency (or energy) at least. This is particularly true of media in which g(z) is negative 
in two separated regions; in this case a whole spectrum may be transmitted. 

The formulae are developed in the symbolism commonly used in electromagnetic 
wave propagation, but the results are also considered in terms of mass waves. 

II Electromagnetic Waves. In the electromagnetic wave problem, an equation of 
the form of (1) comes from eliminating time, t, from 


oy n? ay 
—— oo a eG (2) 
Ox € ot 
by setting 
Y = fixe’. (3) 
One thus obtains 
df a? oe 
— + kn? f = 0. (4) 
dx : 
Here n? = n*(zx) is the effective permittivity, c is the vacuum speed of light, and k is 
the (real positive) circular wave number in vacuum. 
We shall here be concerned only with media for which n’ = n*(z) is real, finite, 


continuous, and such that 


R= I for iFi>ra>d, 
(5) 


for ties 3 


lA 


2 
n 


These conditions are imposed primarily for convenience of discussion, and most may be 
removed if proper modifications in the development are made. 
Introduce r = r(x) by defining 


f = exp E [ r(s) as| (6) 


Fas 
where 2, is an arbitrary real constant. Substituting this into (4) we get 


2_ 2  tdr ~ 
atedied k dx 7) 


as the equation determining r(x). We shall allow only continuous solutions of (7) so 
that f and df/dzx will be continuous. Three solutions of (7) will be of particular interest 
to us; they will be denoted r*, r-, and 7’, and are characterised by the following conditions: 


* = ] for 2S 2; (8) 
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r=-l for 2S —h, (9) 
r=1 for «r>%. (10) 

Writing 
ri =r, + ir (11) 


where 7, and 7, are real, equation (7) yields 





. 1 dr, 
V=n—mnt+ : >, (12) 
ld 
0= 27, —- k 7, (13) 


since we are taking n” to be real. Equation (13) gives immediately 





1 d 
To = og, oe |" | (14) 
by which we obtain 
[rds = fr) ds + & tog | HO | (15) 
J xe ; a "7 2k * | r,(2o) | 
For convenience, take 7 = —z2, so that, by (8) and (11) 


1, (Xo) = r,(—2,) = aS (16) 


and hence we can write 


oe are oy 
f= Fale? E r _n@ as| (17) 
as a solution of (4). It will be shown in Part IV that 7,(z) is always positive, hence the 
absolute value bars may be dropped. 
The expression 


Ft frente! = WAG) exp E [. r(s) ds — iet (18) 


with F* an arbitrary complex constant, will then be a solution of (2). In the region to 
the left of the variation (ie. for z < —z, , by convention) this represents a wave of 
amplitude F* and speed c moving to the right. In the region to the right of the variation 
(cx > 2,), it must represent, in general, two waves, one moving to the right and the 
other to the left, each with speed c. 

It can be seen by inspecting (12) and (14) that if r, + ir, is a solution of (7), so too 
is —r, + ir, ; this is, in fact, simply r~. Hence we have a second solution of (2) namely 


ee rn oe a ae 
F-f-e = ar exp E [ r,(s) ds ict | (19) 


| —z6 


Remarks similar to those above apply to this solution also. 
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The two solutions obtained are independent, the Wronskian of the system f*, f~ 
having the constant value —2zk, hence the general solution of (2) is 


V=(F fi +Ffje* (20) 


with F* and F’ arbitrary complex constants. 


Another solution of (2) is 


‘Sho ee E | ri(s) ds — tet |. (21) 


/ 
V 7,(v) “ze 


Here we have taken x) = x, and?” = r{ + ir}, with r{ and r} real. This solution repre- 
sents, in the region to the right of the variation, a single wave moving to the right. 
III. Formulae for Reflection Coefficients. In obtaining reflection coefficients, we 
consider solutions of (2) which give a single (transmitted) wave leaving the region on 
one side or the other. If this occurs on the left we must have F* = 0, and the complete 
solution is then given by (19). If the single wave moves off on the right, then the ratio 
F~/F* must be such that the left-moving components on the right cancel out. In this 
case, however, we can equally well give the complete solution by (21). For equivalence, 


Pept + Fp = F's’ (22) 


and, differentiating this with respect to x and dividing out 7k, 


From these we obtain the complex amplitude reflection coefficient 


i r —?7f 
Rs— = -—— 
f 7 out aiee 
(24) 
(r, — rt) + ar. — 73) ; F 
=— 1 aT? =~ exp | 27k | r,(s) ds 
(—r, — 73) + ir. — 12) Faia 
which has modulus 
(r, — ri) + u(r. — 75) ror 
R|= e - (25) 


(r, + 73) + air. — 72) 


Since F is a constant, the evaluation may be made at any z. 
If the variation is symmetrical about x = O—i.e.,if n?(x) = n°(—zx)—then it may 
easily be shown that 


ri(z) = r,(—2): r(x) = —7r(—2). (26) 
In particular, at zc = O we get rj = r, and r} = —r, . Substituting this in (25), we get 
R,| = [11 +7nrOmO)]'? (27) 


for a symmetrical medium. 
It may be noted that formulae (22) to (25) would be valid even if n’? were complex 
in the region |z| < 2, 
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IV. Approximate Curves for r* (x). No general solution of the Riccati equation is 
known, so we are apparently no further ahead. However, the sort of solution to be 
expected can be estimated, to a certain extent, by inspection of the equations (12) to 
(14), and this task will be undertaken in the present part. The development is necessarily 
lengthy, and has only been outlined here. However qualitative it may be, it does lead 
to some general results and does provide a guide for further extension, as illustrated in 
succeeding parts. Moreover, since the required curves are shown to be smooth, at least 
over the regions of interest, it appears that electronic computers might be used to ad- 
vantage in this sort of approach to the problem. (Contrast the actual wave solution 
which is, in the main, oscillatory.) The operator D will be used for d/dx throughout 
the discussion. 

We shall be concerned here primarily with fairly slowly varying media, having 
simple curves of n*(x), in which n* is negative over some range, say % <2 <2.A 
representation of such a medium is shown in Fig. 1. Media in which n’ is slowly varying 
but never becomes negative will be classified as Type I, and will receive only scant 


treatment, near the end. 








- x 
—Xa Oi Xa 
Fig. | TYPICAL MEDIUM TO BE CONSIDERED 
We start from the region x < —x, , where 7; = 1 and r, = 0, and proceed to the 


right. Since r* is continuous, no change can occur in either r, or r. until one occurs in 
Dr, . Since n’ starts to decrease, we see from (12) that Dr, must become negative. Then 
D*r, and Dr, become negative, and r? decreases from unity. So long as 3(Dr,)* > 
2r,D’r, , however, r; must stay greater than n’, as may be seen on eliminating r, from 
(12) and (14). This gives, for the start of the rj and r, | r. | curves the form shown in 
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Fig. 2. (We plot r, | 7. | so that the magnitude of r; and the signs of r. and Dr, may be 
shown.) 
Before proceeding farther, we might first note two general properties of the r; curves, 
based on (12) and (14): 
i) r, will never become zero (if 2, is finite) for this would require r, — — ©, and 
hence n® — — ©, a case which has been excluded from the present paper. 
| 
| 





} or dor 
| Xo fe lt, | 


Fig. 2 STARTING CURVES OF f° AND girl 


ii) if ever Dr? becomes zero after having been negative, so too will Dr, , hence at 
that point r, = 0 and Dr, > 0, and so n® > r; . Hence, so long as r; remains greater 
than n’, it will not have a zero derivative. 

For slowly varying media we may expect slowly varying r curves, in which case 
we will have 7; = n’ until very close to x, , where n* vanishes. Such a result is, in fact, 











Fig. 3 POSSIBLE LATER CURVES 1, AND «Ir! 
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employed in the B.W.K. method where, in effect, r, = n (for x < 2) is derived from 
the first term of a series solution of (12) and (14), and then employed in (18). Except 
in the limit of vanishing variation per wave length (k —«) we cannot know that this 
result must be obtained, but the reverse certainly is true: if r, and Dr, are very slowly 
varying, then n’ is too, and rj ~ n”. Later curves of rj and r, |r. | may then be like 
those in Fig. 3 or, more likely if the medium is slowly varying, like those in Fig. 4 or 5. 














Fig. 5 POSSIBLE LATER CURVES OF 1f* AND In 


At any rate, when we come to x, we have n® = 0, r; > 0. Media having Dr, > 0 
at 2, will be classified as Type II while those having Dr, < 0 there will be classified as 
Type III. There does not appear to be any general way of continuing the curves past 
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x, for media of Type II, which we shall therefore drop from our discussion for the 
present. Type III can, however, be treated generally, as follows: 

As before, there will be no point in the region xz, < x < x, where Dr; = 0 = 7, 
since the property (ii) above applies to Type III media there. This means, then, that 
Dr, < 0 and r, < 0 throughout this region. Moreover, if the medium is sufficiently 
slowly varying, 7; will not differ from n* by very much in the region to the left of 2, , 
and hence will be very small at, and to the right of, x, . In fact, it can be shown that 


|r, dx < 4/2k (28) 
Sin 
for any medium of Type III. (This is easily deduced from the fact that the graph of 
the real part of f must be curved away from the z-axis in regions of negative n°.) From 
this we see that r,; must become quite small between x, and x, if k(x, — 2,) is reasonably 
large, no matter how large r, is at 2, . 

Proceeding to the right, we see that if —r; is less than n° by more than the (small) 
amount r; , then Dr, must be positive and —r> must be increasing. So long as n° is not 
increasing, then, —r; cannot be much less than n’ for long. Once —r; becomes greater 
than n° — r; it must remain there, at least until this latter starts to increase (near 2,, , 
the point of minimum n’). All that time Dr, will be negative, and so —r} will be de- 
creasing. If the gap between —r; and n” — r; becomes large, so too will —Dr, , thus 
forcing the gap to diminish (assuming fairly slow variation of n’). Figs. 6 and 7 show 


the type of curves to be expected. 











Fig. 6 POSSIBLE CURVES OF rf? AND f,Ir,1, SHOWING 
-r? BECOMING GREATER THAN n° 


All these considerations hold especially well when the variation is sufficiently slow. 
In the limiting case k © we will have r, — n, r, ~ 0 for z < x, andr, — 0,7 — 
—|n!| fora, <x < x,. As we relax this condition of slow variation, the considerations 
above are modified. Inspection shows that this leads to curves such as those sketched 
in Figs. 8, 9, 10, and 11. They are, of course, only rough drawings meant to indicate 
the sort of curves to be expected; they are not mathematically accurate, and “slow 


variation” is not defined. 
In most cases of slow variation, the curves can be extended as far to the right as 
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Fig. 7 POSSIBLE CURVES OF f° AND alr, 
SHOWING -f BECOMING GREATER THAN n° 


xq Without too much error. After this point, however, it is impossible to continue by 


inspection, for 7, and its derivatives become large. However, we do know that 


exp it | r’(s) as| = Ac’ + Be*™ (29) 


for z > x, . Here A and B are complex constants and, in our case, n = 1. From this 
can be derived 
oo io a | Oe ea 
rige=n— Hoga ge et (30) 
Ae" + Be 


In general (i.e. B # 0) this gives oscillatory curves for r, and r, . 


=== x 











Fig. 8 f° AND pI, CURVES FOR A VERY SLOWLY VARYING MEDIUM 
OF TYPE IL 
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Fig. 9 1? AND fglr,1 CURVES FOR A SLOWLY VARYING MEDIUM OF 
TYPE I 














Fig. 10 7 AND «,l,1 CURVES FOR A SLOWLY VARYING MEDIUM 
Or TPE if 


Actually, there is no need to describe the curves past z,, , since the formulae (25) 
and (27) can be used at z,, itself. We can trace the curves for r{’ and rj | ri |, working 
to the left of x, in a similar way, to the region of the minimum of n’. This gives a com- 
plete set of curves for the medium, such as that shown for a typical case in Fig. 12. 
Other media will have curves differing in detail from these, but we may take it as a 
fairly general result that media having a single fairly wide region of negative n’ will 
have 7, and r{ quite small, and —r, and r; approximately equal to | n |, at z,, . 
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——— 


Fig. 11 6° AND «1,1 CURVES FOR A FAIRLY RAPIDLY 
VARYING MEDIUM OF TYPE IL 


This last statement will be true even for media of Type II, since it can be shown 
for them that ; 


k r, dx < m/k. (31) 


“2b 


If, then, n® stays negative for long enough, r,; must become very small.“Along with 











Fig. 12 CURVES OF rf girl, ri AND igi FOR A FAIRLY SLOWLY VARYING 
MEDIUM OF TYPE I ON BOTH SIDES 
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this, Dr; must become negative, and the results derived above for Type III media will 
become applicable. 

In Type I media having n* slowly varying, we expect the approximation 7, = n, 
r, = 0, to be fairly good throughout the entire range, at least if the minimum n° is not 
too close to zero. 

V. Approximate Reflection Coefficients. i) Type I media, minimum n*® > 0: So long 
as the approximations r; = n = r{, To = 0 = rj are good, (25) gives | R| = O. If these 
approximations break down, as will happen in rapidly varying media, the resulting 
reflection coefficient may become appreciable. 

ii) Type I media, minimum n* = 0: It can be seen that, under certain circumstances 
(eog.7, = —,7 = 7s, Dr. = 0 = Dri at z,,) we will obtain| R| = 1 V2. This is a 
result frequently quoted for this case (e.g. Rydbeck [13]), but as can be seen, it is by 


no means general. An example will be given (Part IX, Fig. 20) in which it completely 


breaks down. 


iii) Media having a broad region of negative n°: We will have, at x, , —T2 > 7 > 0, 
r, > rl > 0, and hence |R| = 1. We thus get the expected result, that such media 


are almost completely reflecting. 
iv) Symmetrical media: On integrating (14) and introducing the result into (27) we 


obtain for symmetrical media 


— a + ge 
hi = [ — r2(0) exp q tk | ni Ts dx | ‘ f< 


This gives a form which may be compared with the common (approximate) formula 


R = [ + exp ( — 2k | | n ax} 33) 


which becomes, for symmetrical media 


Re =|) +exp{—4k | |n az} | P 


This was originally developed by Gamow [14] by considering discontinuous media, and 
more recently by the “good path” method [6] [13]. If n° is sufficiently negative for a 
sufficient range of x, both methods give | R,| ~ 1, as in (iii). However, for rapidly 


varying media which are not too thick, (32) [or in general, (25)] and (33) may give 


w 
bo 


quite different results. 

v) Thin media of Type II: The reflection coefficients for such media may differ 
considerably from the “good path” result, and so somewhat surprising values may occur. 
An example is given in Part IX, Fig. 21. 

VI. Accurate Reflection Coefficients. Regardless of the accuracy or inaccuracy of 
the discussion in Part IV when applied to general cases, the formulae (25) and (27) 
can be used for specially constructed media as follows: 

Use any function p,(x) which equals 1 for « < —z, , and remains positive and has 


a continuous first derivative everywhere. Define 


ld : 
p2(x) = Ie ym log Pi (34) 
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and construct a medium for which 
2 2 2 ld 25 
n = pi — pot k dy ®? (35) 


using any specified / 
For this medium, at the frequency kce/27, we will have 


ry * Des To = Pe (36) 
exactly. To satisfy the conditions (5) on n, p, would have to have a certain form for 
x > x, . This difficulty can be avoided, however, by using another function pj which 


equals 1 for « > x, , and remains positive and has a continuous first derivative every- 
where. The medium to be used will have n* given by (35) for x < x, and by a similar 
relation in p{ , ps for x > 2, , x, being a point where the two relations give the same 
value. (The existence of at least one such point is a further condition, though we may 
now remove the restrictions on p, to the right, and those on p{ to the left, of x, .) 

We should then be able to derive the exact reflection coefficient for the medium so 
constructed, at the frequency kc/2x. Using the discussion of Part IV as a guide, one 
could obtain a specified form of medium. For media with large regions of negative n’, 
it is easier to choose p, and p{ over part of the range only, and choose p, and p} over 
the remainder (x, < x < 2,). Examples are given in Part IX. 

It should be pointed out that this technique is accurate in all cases. The variations 
may be as rapid as we please, and n° need not be restricted by the conditions imposed 

this paper (see, for example, Fig. 26, Part IX, where n ¥ 1 on the right). Even com- 
plex n-’s may be obtained, if proper modifications are made. There are certain dis- 
advantages to the method, however: 

i) The reflection coefficient is only known for the frequency used in the construction. 
For nearby frequencies we may expect little change, but we have no way of estimating 
the amount. In some cases (see Part VIII) the frequency is very critical in determining 
reflection. However, if we choose functions to give curves as in Fig. 12, we can be fairly 
sure that for higher frequencies the curves approach the limiting values mentioned 
before (aside from dispersion effects). 

ii) The exact form of the medium, i.e. 


n* = n*(x), is not known until after the con- 


struction is made. However, with proper judgment, we could construct almost any 
medium desired, at least approximately. 
iii) In general, dn” /dzx will be discontinuous at x = 2, . This is not a serious fault, 
however, and need not occur if care is taken in selecting the p functions. 
VII. Non-Reflecting Media. One of the primary values of the method just developed 
is that it shows the existence of non-reflecting media (for any given frequency) even 
is considerably negative over a considerable distance. To see this, we note 


when n™ 
that if p, is taken so it — 1 for x > x, , then we may take p! = p, , hencer, = r{, 
1 r, , and hence | R| = 0. Using the great freedom of Hides of p, allowed, we can 


construct a great variety of rn for which the reflection coefficient is zero (for the 
frequency used in construction). This applies to rapidly varying media as well as slowly 
varying ones. The possibility of non-reflecting media was noted by Rayleigh [1], based 
on similar reasoning, but not developed by him nor applied for cases of negative n’. 

If we choose p,; so that p. is appreciable when p, = 0 we can get a medium with n® 
appreciably negative (and for as great a distance as we wish) which has the property 
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of non-reflection (for the frequency used). Several non-reflecting media are sketched in 
Figs. 13, 14, 15, and 16, the last two being of this type. 
All non-reflecting media of the types considered in this paper must have 

nto 

| r,(2) dx = 0 (37) 
in order that 7,(z.) = 7,(—2,). Since regions of negative n* usually arise due to | rz | 
becoming large, we can expect that a non-reflecting medium having one region of 
negative n’ (due, say, to — 72 becoming large) will have another one (due to +7, becoming 


——__ 





| 
} 
| 
| 








Fig. 13 A NON-REFLECTING MEDIUM (FOR A CERTAIN FREQUENCY). 7° CURVE 
PRESCRIBED 











14 A NON-REFLECTING MEDIUM (FOR A CERTAIN FREQUENCY). ° CURVE 
PRESCRIBED 


Fig. 
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Fig. 5 A NON-REFLECTING MEDIUM (FOR A CERTAIN FREQUENCY) WITH TWO 
REGIONS OF NEGATIVE n*. rf, CURVE PRESCRIBED 


f 











Fig. 16 A NON-REFLECTING MEDIUM (FOR A CERTAIN FREQUENCY). WITH A 
REGION OF NEGATIVE n°. 15° CURVE PRESCRIBED 


large) or else one in which n’? =~ 0 for a long distance (due to r, being positive but small 
for a long distance, , growing only slowly). The two cases are illustrated by Figs. 15 
and 16. In the case of very thin regions of varying n’, in which 7, and r, may vary 
rapidly, these results need not follow (as, for example, Fig. 21, Part IX). 

For media which do not satisfy (5), equation (37) and the results derived from it 
need not hold in eases of non-reflection. This is illustrated in Part [X, Fig. 26. 

VIII. Reflection from Media with Two Regions of Negative n°. General media 
having n’ negative in two separated regions cannot be treated by formulae (25) and 
(27) directly, because of the difficulty of extending 7, , 72 , 7; , and r; curves to give 
values at any one point. Useable formulae for this type of media can be obtained by 
the following extension of the previous methods: 
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Introduce two more solutions of (7) 
p (x) = pi + tpn ; p (x) = —p, + tpz (38) 
where p, and p, are real, such that 
DP =|ni, Po = dp./dx = 0 at r=2 (39) 
where 2, is the point at which n” is a maximum (>0) between the two negative regions. 
Curves of p, and p, may be drawn by the same reasoning as used for r, and 7, , at least 
if | n(x.) | is not too small. This gives curves such as those shown in Figs. 17 and 18. 
By introducing corresponding solutions of (1), and proceeding as before, we get the 


complex reflection coefficient 


7 _ ‘ _ Oy}. [22 il _ > 
R = exp | 2 | r, i | A exp | = J; p, dx) a, (40) 
/-z¢ C exp [—27k 2? p, dx] — D 











Fig. 7 CURVES FOR A MEDIUM WITH TWO REGIONS OF NEGATIVE n° 


where 


B= (rt — p )soe,(P  — 1) a2, 


C= (Pr — Pp )sos.(P — 1)ses. 


/ 


Da (tf — p )eos,.(P — eas, 


and 2, and z, are any two values of x. For purposes of evaluating (40) we would take 
Z, near z,, , the point of minimum n’ on the left, and z, near z, , the point of minimum 
n® on the right. 
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aa 











Fig. 18 CURVES FOR A SYMMETRICAL MEDIUM WITH TWO REGIONS OF NEGATIVE n° 


n slowly varying media, with sufficiently deep minima of n?, we will have —r. >", , 
I ; ; } , 


Po > p, atz,,andr; >r{, —p. > p, atx, . Taking 27, = z,, , 72 = 2, , as We Shall from 

now on, we find A = B = C = D and hence, in most cases, | R | = 1, the expected result. 
However, if two conditions are fulfilled, namely 

A|=/|B (41) 

ga — 2k | Pp, dx = gg — 2qr (42) 

then we will have, from (40), | R | = 0. Here g, = phase angle of A, gs = phase angle 

of B, and g = 0, +1, +2,---. 

If the medium is symmetrical, | A | = | B| for all frequencies, as can be shown by 

choosing x; = —2, and comparing the quantities involved. The second condition, (42), 

will probably be fulfilled for a whole succession of frequencies. If the medium has curves 


like those of Fig. 18, where —r.(z,,) >> 7(2,,) ete., then ¢, = gs = nearly independently 
of frequency. Thus | R | will vanish every time 2k [7* p, dx increases by 27 as k varies. 

It can further be shown that, for such a medium, the band width Ak of passed waves 
(those having | R| < 1/V 2) is given by 


Ak | p, dx = tm, | : (42) 


(p2 — 12)” 


ri 


On the assumptions made in developing this formula, Ak will be a very small quantity. 
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These media act as very narrow band pass filters, then. For such a medium, a graph 
of | R | against 2k {2? p, dx would appear much like that sketched in Fig. 19. 

For media having the minimum of n° not very deep, the maxima of | R | would not 
be so close to 1, nor would the dips in the curve be as narrow. In a plot of | & | against 
k, some distortion of the graph will occur. This comes from the change of p,(x) with k, 
the changes in yg, and gy, with k, and a possible dispersion n° = n*(k). 

Although these results have been derived for symmetrical media, there is no reason 
to believe that non-symmetrical media need behave much differently. Possibly | R | 
would not become zero in most cases, even for selected frequencies, but it would probably 
become small for a sequence of frequencies, with intervening values near 1. Certainly 
we can construct non-symmetrical media for which there is no reflection, and in such 
cases we might expect a marked decrease in | R | for a whole sequence of frequencies. 
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Fig. 19 REFLECTION COEFFICIENT FOR A SYMMETRICAL MEDIUM 


The condition (42), and the very nature of media with two regions of negative n’, 
remind one of the theory of the Fabry and Perot interferometer. This suggests as the 
explanation of perfect transmissions through these media a process of multiple reflection 
with subsequent interference. It further suggests that the phenomenon will occur to 
some degree even if the media are somewhat absorbing. Comparison should also be 
made with the problem of periodically varying media, as treated in Mathieu’s equation, 
where a sequence of bands may be passed. (See Brillouin [15].) 

IX. Accurate Reflection Coefficients in Particular Cases. Fig. 20a shows the macro- 
scopic appearance of a symmetrical medium having minimum n’ = 0. This was drawn 
on the basis of an assigned 7; curve for | z| > 10, using k = 100 (i.e. Aya. = -063 units). 
This figure actually shows n” for three media which differ almost inappreciably from 
one another in the region | z| < 3 and are identical outside this region. For |x| < 10, 
n’ is derived by specifying r, (which is continuous, at z = 10, with the r, curve outside). 
The three media, with the corresponding r, curves, are shown over the range —6 < 
x < 0, on a greatly enlarged scale, in Figs. 20b, c, d. They have values of | R | equal 
to 1/2, 0, and .92 respectively. It is evident from this that very slight changes in 
the form of n’ can change the reflection coefficient remarkably, and that 1/2 is by 
no means a general result for these media. 
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Fig. 21 shows a thin, fairly rapidly varying medium, having n’ negative in one 
region only, and having zero reflection for the frequency used in construction. The 
curve of 7, (not shown) was specified, and k taken as 1 (A... = 6.3 units). In coming 
from the left, it can be seen that this medium is of Type II and that Dr, remains positive 
1.0 
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throughout the region of negative n’. In this case, [ii r, dx = 0.34 which is less than 
w, in agreement with Eq. (31). The reflection coefficient modulus, based on the good 
path method, is 0.79, which is in obvious disagreement with the correct value, ‘zero. 

Fig. 22 shows, on the right, one half of a symmetrical medium, and on the left, the 
curves of 7; and 7, | r. | from which it was constructed. Here k = 10, \,.. = 0.63 units. 
The reflection coefficient for this wavelength has modulus which differs inappreciably 
from 1 when calculated by either (32) or (33). The transmission coefficient has modulus 
approximately equal to «°°; the good path method gives « ”’. 





Fig. 23, shows similar curves for another medium. Here, k = 1 and ,,. = 6.3 units 
for the construction. For these waves, | R | = 0.954, while the good path method gives 
|R | = 0.998. The medium is, of course, rapidly varying and “narrow’’ relative to the 
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wavelength used. Figs. 21 and 22 may be compared to see how closely the r curves 
“stick to’’ the n curves in the two cases, “fairly slowly” and “fairly rapidly’’ varying. 
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Figs 24a and b show curves for a medium such as occurs in ionospheric work. The 
value of k used was 1000, corresponding to dy. = 0.0063 units. For an isonosphere 
120 km. thick, this corresponds to d,,. = 63 m, or a frequency of 4.8 mc/s. Reflection 
is essentially complete here. With such large values of k (for media of this thickness) 


there is practically no deviation of the r; and r, | 72 | curves from the limiting forms 
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mentioned previously, and almost identical curves would be obtained with any k of 
this or greater order of magnitude. In the case of the ionosphere, of course, the medium 
is dispersive and there would be a change of the n’ curve as k is changed. Inter-modal 
coupling must also be included before a complete solution is obtained in this problem. 


FI.O 




















4644 *&4/4 5 6 
Fig. 240 
~~ £ +1.0 - 
a --- ff p* ond r! Ps 
s 
\. — tle! p,lp,land gini-os Ff 
\ j 
\ r0.6 / 
\ a / 
\ iy 
+0.4 / 
* \ / 
\ \ / 
\ 0.2 / 
\ ~ \ / 
\ + / 
\ a. ak 
T T . 2 7. . t T Tt u ' 7 ’ 
6 -5 -4 -3-2 -I \ I 2 3 4 5 6 
~ poz 
\ 
sae 
Koa. 
Fia. 24b 
Fig. 25 shows, on the right, half of a symmetrical medium having two regions of 
negative n: on the left are curves used in deriving this medium. A value of k = 10 
-—. ‘ F1.0 
~ um , 
\ s== fond pF fos 
\ —-—,lr,Jand p,!p,! 
\ +06 
\ 
\ Lo.4 
‘ 
\ ~ -0.2 
\ 7 
\ Yh ' 
6 7 4 BL 4 i 2 4 . 
“~~ 
O02 











30 C. O. HINES [Vol. XI, No. 1 


was used (A... = 0.63 units), but similar 7; and r, | r. | curves would be expected for 
any k near this. Here [7* p, dx = 0.36, and at 2,7, = € °°", rg = 0.4, p, = €°™, 
and p, = 0.4. These approximations will be fairly good for all k near 10. Then g, = 
¢e = 7 as stated before. There will then be no reflection of waves having k = 8.7, 
17.4, --- (= 2.8 gm with g = 1, 2, ---). The band widths passed are given by Ak = 
4¢~*-**/(0.8)7(0.36). For k = 8.7 this is Ak = 1.1 X 10°° and for k = 17.4 it is Ak = 
7 X 10°"*. These are obviously very narrow bands. 

Fig. 26 shows a medium which has n* negative only once, and which has n° = 1 
on the left and n’ = 0.003 on the right. There is no reflection for the waves used in 


construction, which had k = 1 and X,,, = 6.3 units. 
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X. Mass Waves. In quantum theory, an equation of the form of (1) comes from 


eliminating ¢ from 
. OW h ow , - 
th — = —-—— 3 + Vi(a)¥ (43) 
at 2m Ox 
where h is Planck’s constant divided by 27, m is the mass of a particle, and V(z) is a 
potential energy field which varies only with the x coordinate. We set 
iEt/k (44) 


WwW = f(xe 


where £ is to be the (constant) total energy of the particle. This gives 





os = [E — V(a)]f = 0. (45) 

The forms of (45) and (4) are identical; the solutions will be identical if we set 
k’n? = [E — V(x)]2m/h’. (46) 
The factorization of the right hand side into k* and n* parts may be done arbitrarily, 
but the &? part must be independent of x if our previous equations are to hold. If V = 0 


for |z| > zx, , then the most convenient association is 

k? = 2mE/h’; n = 1— V(x)/E (47) 
since this gives n? = 1 for |x| > x, as before. The other restrictions previously placed 
on n’ correspond to V being a real continuous function of x which is non-negative be- 
tween —<2, and z, . This is then, the potential barrier of quantum mechanics. The cases 
of maximum V being less than, equal to, or greater than, FZ, correspond to minimum 
n* being positive, zero, or negative, respectively. 
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The results of all the previous parts of the paper may now be applied to the present 
problem. One modification only is important; we are usually interested in the behaviour 
of the solutions as Z is varied, keeping m and V fixed. This corresponds not only to a 
variation in k but also to a variation in n’—i.e., the media considered are dispersive. 

Physically, the results of this paper when applied to mass waves indicate that streams 
of particles of the proper energy can get through some potential barriers without any 
reflection. In the case of double barriers, a whole spectrum of such beams could pass. 
It may be possible to apply this result to obtain sharp ranges of velocities for such 
experiments as Thomson’s e/m determination. 

Since the radial wave equation can be put in the same form as (45), the methods of 
the present paper can be applied in radial problems. For example, exact eigenstate 
conditions can be established, in terms of r(x), for diatomic molecules. The potentials 
occurring in such problems do not lead to the type of n’ curves considered in this paper, 
however, so the matter will not be pursued here. 

XI. Remarks. The results of this paper are in general agreement with those of other 
methods. However, by showing that unexpected reflection coefficients may occur, they 
contain a warning against using the formulae of the other methods in cases where they 
have not been proven applicable, such as media with two regions of negative n’, and 
media having roots of n° in the complex plane near those on the real z-axis (see [15]). 

This investigation was carried out at the Radio Physics Laboratory of the Defence 
Research Board, Ottawa. The problem was proposed by Mr. J. C. W. Scott, whose 
suggestions and encouragement are much appreciated. 
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CONNECTION FORMULAS BETWEEN THE SOLUTIONS OF 
MATHIEU’S EQUATION* 


BY 
GREGORY H. WANNIER 
Bell Telephone Laboratories, Murray Hill, New Jersey 


Abstract. The problem of connecting the various types of solutions of Mathieu’s 
equation is solved by the introduction of a new parameter © which is a function of the 
the two equation parameters a and gq. This quantity ® is introduced and enclosed be- 
tween two very close analytic limits in section 2. In sections 3, 4, 5 precise definitions 
are given and information is collected for the three main types of functions which are 
to be connected. Section 6 contains the connection formulas. Section 7 reviews the 
status of knowledge achieved. Section 8 is an appendix on integral equations which 
are more general than those developed earlier in the text, but which appear to be of 
no use for the main purpose of this paper. 

1. Introduction. A variety of different types of solutions have been written down for 
the Mathieu differential equation 


d’f 
7 u + (a — 2q¢ cos 2x) f(x) = 0 (1) 
aa 


a general solution of which we shall call me x. The special solutions proposed stress 
different qualitative features of this general solution and suggest themselves in different 
types of applications. By general principles, there must exist a connection formula 
between any three such solutions. It is the purpose of this paper to write down these 
connection formulas for some of the more important solutions of equation (1). 

In carrying out this program we shall assume q real, and once it is assumed real it 
may be assumed positive because the transformation 


will reverse the sign of g. We shall express this sometime by writing 
q = kK’. (2) 


The parameter a will also be assumed real. The treatment is particularly designed for 
positive a; this becomes important in the discussion of the next section. 

2. On the function ke y. Equation (1) contains three distinct real equations, one of 
which is (1) itself. The second is obtained by setting 


which yields 


d’f 
dz 


— (a — 2q cosh 2z)f = 0 (3) 


*Received March 25, 1952. 








34 GREGORY H. WANNIER [Vol. XI, No. 1 


and the third by setting 
e= a + 1y 


which yields 
d°f 


35 — (a + 2¢q cosh 2y) f = 0. (4) 
dy si 


The study of equation (4) for f as a function of the real variable y will occupy the rest 
of this section. 
Whenever | y | is large equation (4) can be approximated closely by 


df sivi « 0 a 
> — = (¢ 
dy? (a + qe dj )) 


which admits as solutions the modified Bessel Functions 
fly) = K,::(ke”’), fly) = Larvlke”’). (6) 


The solutions (6) give the asymptotic character of the solutions of (4); they show us 
that there must be one solution vanishing as K,:/.(k e”) for positive y. We shall call it 


ke y; we shall normalize it by the prescription that 


exp (—k¢ 
(k e’)'”” 


for y> 0. (7) 


ke y~ 


We now assume that in equation (4), we have the inequality 
a+ 2q> 0. (3) 


The solutions of (4) are then always curved away from the axis. This means for ke y 
that it stays positive throughout, and that it cannot vanish exponentially for large 


negative y. Hence we must be able to write 


> [}: , + D/ ~ ae 5 
key ~ : —_ a : — K)] ior y «Kx 0 (9) 
(ke) 


with & (a, k) real. 

We shall now show that under the restriction (8), (a, k) is a smoothly varying 
function of a and k, that is, it does not partake in the oscillatory character of the param- 
eter 8 introduced by Floquet’s theorem. We shall establish this fact in the remainder 
of this section by enclosing ® between narrow limits, neither one of which shows any 
oscillation. These limits are reproduced in Figs. 1 and 2; they furnish incidentally a 
close numerical approximation to when desired, The structural discussion of Mathieu’s 
equation will then be resumed in section 3. 

We start out with the approximate value of @ which is obtained by application of 
the Jeffreys’ or WKB method’ to equation (4). We set 


key = A(y)e"*“ (10) 





1H. Jeffreys. Proc. Lond. Math. Soc. 23, 428 (1924). 
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LOWER BOUND RESULTS FROM THE JEFFREYS APPROXIMATION. 
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FIG.2 UPPER BOUND ¢ OF $(a,Q). THE BOUND IS NOT PROVED 
TO THE LEFT OF THE WAVY LINE, BUT PROBABLY HOLDS THERE 
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and set 


(48) = | 99 eosh 2 (11) 
dy = &@-+ “£q COSN zy. \ ) 
This yields for A the equation 

PA A dS *S 

ff gS .. 4S nw 0. (12) 

dy dy dy dy 


The standard solution is obtained by neglecting the first term in (12). We find in this case 


exp [4@, — [3 (a + 2k? cosh 2n)'”” dy] 
ke > ae. em — ag! Meer key mes 9, ‘ ( 13) 
(a + 2k cosh 2y) 


The term 44, enters into the exponent in order to satisfy the asymptotic requirement 
7); it obeys the relation 


5 Po = lim (a + 2k? cosh 2n)'’* dyn — 2k sinh y |. (14) 


Equation (14) then represents the Jeffreys approximation to ® as defined in (9): 
$= g. (15) 


®, is easily evaluated in terms of complete elliptic integrals; we find 
a — 2k’ fa — 2k°\'” =e 
6, = 2——=353 D (S— 3 (16a) 
(a + 2k) a+ 2k°/ 
fora < 2k 


? —aJ,,| (2k? — a)” (2k? — a)'” or 
= Maafy{@e=a) fara" ay 


The two expressions (16) do not indicate a break at a = 2k’, for they are analytic con- 
tinuations of each other. 

The approximation (16) to ® is shown in Fig. 1; we shall now prove for it that the 
approximate identity (15) is an inequality, that is, that 


fora = 2k 


> f. (17) 
To prove this let us write equation (12) in the form 


GS. o,94dS_ FA _ og 


dy” ~~ dy dy . dy’ 7 


Subtracting (dA/dy)’* on either side we get 


d » aS aA) (d4)" 

f(g 4) = -|— ). 

dy ( , dy , dy dy <6 
Integrating this inequality between —o and +, and observing that A(dA/dy) 
vanishes at either end, we get 


dS) _ (4248) 
(4 is) 4 dy! _. iadiane 
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If the Jeffreys’ approximation were correct this difference would be zero; actually, the 
inteiaihan 7), (9), (10), (11) and (14) yield 


' is exp i] [2 Po} 

A(+o0)~ (dS/dy)” (18a) 
¥ exp {—3@ + 9} ; 

At ) (dS dy)’ ba ‘ (18b) 


These expressions reduce the inequality to 


which is equivalent to (17). 
In order to gain an upper limit for &, we transform equation (12) by the substitution 





; ~ tis (19) 
which yields the equation 
8 5 anit = 
and hence the inequality 
dB dS a8 


- — 2B — — =; < 0. 
dy ° dy dy « 


Multiplying with e ~” and integrating we get from this 


~28 wee 
¢ B+ [- e* au 


“9 


Returning to A by (19) we get 


l dA Ss da -2 as 
oS 4 s7° dy > 0. 
iw? / ca? 


Integrating the second term by parts we get this in the form 


d dS\'? l os [f -2s d (ES/av") 
ay E (28) | T 9° | © dy \ dS/dy sl 


v 


With the help of the equations (18) this is finally transformed into 


ach [Tate ay [7 ern & (Slee) 
b— d < af 4 dy dS/dn dn. (20) 


The inequality (20) is the desired upper limit and could be evaluated by numerical 
methods. However, we shall proceed instead to majorize the double integral by explicit 


analytic expressions. 
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In discussing (20) we observe first that the expression in the integrand 


d (fS/as" _ 4q(a cosh 2n + 29) (21) 


dn \ dS/dn/ ~~ (a + 2q cosh 27)” 
is positive everywhere and vanishes strongly at either infinity; the inner integral is 


therefore majorized by replacing e~**‘” by its largest value e-**’. This procedure is 
of no use at this stage as it gives a divergent result. However, if we integrate the outer 


28 


integral by parts according to the scheme 


2S(y) 2S(y 
e dy ~ | d(e = 
| d ; ) 2(dS/dy) 
then we verify by this method that the integrated out part vanishes. The expression 


(20) then takes the form 


“as fi Boal (£8 ay’ 
® ay q i } | ds dy dy dS, dy dy 


ee" Pee ase YY sn oe (es ie) 
——C (ae + reel 1. > 
T 4 | . (ds dy)? ” | dn \ dS/dn ” 


The first integral is positive, as can be seen by another integration by parts: 


a+ as dy’ 


we: 
(dS ay) Y 


b- 4 <i] 
(22) 


bot. eee enc. esi @ (d’°S/dn’ 
+ - | - fl "ea | e 25 — [= : a dn: 
(dS/dy) Ie dn \ dS/dn 


The sign of the second integral is not immediately obvious, because the integrand 
changes its sign with d’S/dy’. An easy way to obtain an upper limit of @ is to replace 
the integrand by 0 in the range in which it is negative; we find thus 


1 6*” (dS/dy’)” Lt" €al sec me (é8 a 
— ¢ ” einen ee 2 9 v) 2 -_ pad f 
® Po < 4 | » (dS/dy)” dy + 4 | dS/dy) Pr dy / dn \ dS/dn dn 


Having now the integrand positive throughout, we can majorize it by the trick dis- 
cussed above of replacing the exponent S(n) by S(y). Observing that 
d° S/dn- 


lim “F§/dy = 1 


we get 


1 6*° d&’S/dy’ 1 f** d&S/dy | d’s a | 
— G, - — + - ae = ren LE 
oS | I. (dS/dy) dy 4 I, (dS/dy)” dS/dy | “Y 


l 


= [ ° (a? S/dy’)’ — 
(dS/dy),<0 


] 
~ 4 Jo (dS dy)* dy + 4 
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By the observation that 


|@S| _ | ds 
| dy? | ~ | dy| 


the first term is seen to be smaller than the second. We get thus with (11) the inequality 
& < &, (23) 
where 
&, = & + = (24) 
2(a + 29) 


The upper limit (24) is shown in Fig. 3. This figure makes the limit appear rather close, 

















FIG.3 UPPER BOUND $5 OF $(a,qQ). THIS BOUND IS 
EASILY PROVED FOR ALL VALUES OF @ AND q. 


but a detailed study for small a and q shows deviations. We shall now establish a tighter 
limit by showing that the second integral in (22) is always negative and can be dis- 
carded, provided we introduce the restriction 


a < 5¢. (25) 


Introduce the abbreviation 


ny _— ,2S8(y) r 72S (0) dad (8/ds") 
E(y) =e / é dn \ dS/dn dyn. 
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Because of (21), this function is positive everywhere and vanishes at +o. It is seen 
to obey the differential equation 
dE 


9 dS | 
dy ‘: dy ss oe 


where 
Bue ™ ayes S (3 ww) 
‘lim = O78 dy) dy \ dS/dy /° 


In an E,y-plane as shown in Fig. 4, the differential equation above defines a slope at 
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FIG.4 DIAGRAM FOR THE STUDY OF THE AUXILIARY FUNCTION E(y). 
THE SHORT LINES INDICATE SLOPES, THE SOLID CURVE CONNECTS 
POINTS OF ZERO SLOPE. THE CURVE E(y) |S SHOWN DASHED. BY 
CONSTRUCTION, ITS VALUE FOR A POSITIVE Y IS SMALLER THAN 
IT iS FOR -y- 


at every point. The curve E = E£,,;,, divides this plane in two parts. Above this curve, 
the slopes are positive, below negative. The curve E = £,;,, is always crossed with 
zero slope. The inequality (25) is needed at this point because if it is satisfied then we 
find from (21) that £,;,, has a slope whose sign is opposite to that of y. When (25) does 
not hold then E develops two humps as shown in Fig. 5. In the case of Fig. 4, E starts 
out by being 0 for + ~, rises for positive decreasing values of y, but must stay below 
the line E = E,;,, because of the slope requirement; thus £ still rises as y becomes 
negative. At some negative y, E reaches its maximum as it crosses the line FE = Lin ; 
then this same slope requirement forces E to stay above the latter curve while going 


to zero. The result of this behavior is that we have for all y 
0< E(-|y!) < E(+/ y)). 
The second integral in (22) now takes the form 


ana << ae 


or, for symmetry reasons 


l eh as dy” sia 
‘ 7320 ‘ > a —_— 1(— ) t ; 
te) | _. (dS dy)” (B Y) Ki Y) 5 dy 
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This last expression is negative because the curly bracket and d’*S/dy* have opposite 
sign, and thus the second integral in (22) is proved to be negative. 

The proof given breaks down if the inequality (25) is reversed. The reason for this 
is shown on Fig. 5; as soon as £;,, develops a minimum at y = 0 instead of a maximum 
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FIG.5 DIAGRAM ANALOGOUS TO FIG.4 FOR THE CASE WHEN 
NEQUALITY (24) IS REVERSED. THE FUNCTION E(y) 1S SHOWN 
IN DASHED OUTLINE. ITS MODE OF CONSTRUCTION NO LONGER 


IMPLIES THAT E(iyl) 2 E(-Iy)). 


the method of constructing E outlined above may give E a positive slope for small | y |. 
Thus the needed inequality does not hold everywhere. Nevertheless it still does hold 
in ‘the greater part of the interval. Numerical checks actually indicate that the second 
integral in (22) is always negative. 
We have thus found, at least under the restriction (25), and probably everywhere, 

that 

& < 9%, (26) 
where 

1 f** (d’S/dy’)’ 
ss lecttcataall [. (aS/ay* 








or explicitly 


fora => 2k’ 

(2 _a— 2 a ) o| (2 _ a") 

“ P 2)1 Drew 2/ .2\3/2 ye - 
(a + 2k 3(a + 2k”) a + 2h (27a) 


1 1 Ila — 2 "| 
+ 3@ 4 2k)” K| (2 Tor) | 





G, 


II 





fora < 2k’ 


(27b) 





2k? — a a oe — ad 
+ ( t rar + a/?L 2% 
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This upper limit is shown in Fig. 2, together with the line (25) above which the theorem 
(26) has not been proved. These established limits determine fairly closely the behavior 
of (a, k°). 

It is possible without difficulty to come closer to the numerical value of ® than Figs. 
1 and 2 do. For this purpose, equation (12), which is linear in A, may be solved in suc- 


cessive approximations beginning with 


dA,dS ads 
2 es 7~=0 
dy dy dy 
which yields (13); successive additive corrections A, , A, , A; , -+- are then obtained 
from the recursion system 
: dA, dS ads "A: ae 
dy dy '~ "dy dy 


If this calculation is carried out up to A, , and if we take the corrected A in the form 
A, exp A; /A, rather than in the more obvious form A, + A, , then a corrected © results 


which is exactly the arithmetic mean of the two limits ®, and ®, 
© = 1(@, + 4,). (28) 


Thus a good approximation to the value of © is obtained by taking the arithmetic mean 
of the readings on Fig. 1 and Fig. 2. 

3. The Lindemann-Stieltjes Functions. The following three sections define and dis- 
cuss the solutions of (1) between which connection formulas are to be established. These 
discussions contain a good deal of information which is already available but which has 
to be combined with the new material to yield the desired results. 

The Lindemann-Stieltjes functions’ are based on the symmetry of the equation (1) 
about the points x = 0 and x = 2/2 which are regular points of the equation. There 
must thus exist an even and an odd power series solution about either one of these 
points. We introduce the following definitions 


I. ce (x; a, q) shall be the even function about x = 0; its value at z = 0 shall be 1. 
II. se x shall be the odd function about x = 0; its derivative at x = 0 shall be 1. 
III. dex shall be the even function about z = 2/2; its value there shall be 1. 

IV. te x shall be the odd function about x = 7/2; its derivative there shall be —1. 


ce and se always form a linearly independent pair, as do de and te. Whenever the periodic 
Mathieu functions ce,x or se,x exist they are respectively identical with the generalized 
functions ce x and se x defined here. Similarly, cé2, or sé2,,; , When existing, are identical 
with de x, Cé2,.; OF 82, With te x. The simplest realization of these functions is by the 
power series method: 


ce x has a series in even powers of sin z. 


= oo = in odd 
ma ~ even ie ™ 66:2, 
loge % & “6&6 Odd sc ‘4 
2E, T. Whittaker and G. N. Watson. A course of modern analysis. Cambridge University Press, 


Fourth Edition, Section 19.5. 
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The coefficients are obtained directly from (1) through term by term recursion with 
the first coefficient equal 1. However, the radius of convergence of these series is only 1, 
and this excludes the possibility of getting connection formulas by this procedure. 

As these functions are constructed relatively easily from others not having their 
symmetry we will not investigate their structure further, but get it indirectly from 
the study of the other types. 

4. The Mathieu analogues to the Hankel functions. We will now use the function 
ke y of section 2 as an auxiliary function to define solutions of (1) having given asymptotic 

a, hex, by the 


’ . : . . 1 2) 
character. We define four such functions, namely he'’a, he~’x, he 


supplementary prescriptions 


(7 ' 
he (z + iy) = ke y, (29) 


he” (~z + iv) = key, (30) 
) = key, (31) 


3 T . 
he' = — 


he“ (: — iv) 


y being taken as real initially. We now drop this restriction and continue the functions 
outside their defining lines. As equation (1) is free of singularities for finite z, equations 
(29)-(32) hold then everywhere and the four functions he“'’z are related to each other 


~~ 
~ 


Y, (32) 


by symmetry operations of equation (1). 

The argument used to establish (6) can be repeated along any line parallel to the 
imaginary axis. As soon as we are sufficiently far away from the real axis, the general 
solution me x will behave as 
exp [+2zk cos 2] (33) 

[2k cosax]'”* * a 





mer~ 


lowever, a particular asymptotie form 


l , , 
= 172 (A exp [27k cos x] + B exp [—2zk cos 2}) 
2k COS & 
can define a given function only within a strip of limited width parallel to the imaginary 
axis. One way to see this is by observing that the above expression is formally periodic 
in 2 with period 47, while by Floquet’s theorem equation (1) has generally no such 


solution. 
In order to obtain this range we start out by proving the integral relation 


| exp [—2k sinh y cos u]me udu = A ke y (34) 


i is an arbitrary solution of (1), ke y the special solution of (4) defined by 


" here irl 
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(7), and A is a number whose value depends on the choice of me u. As integrals of the 
type (34) have been discussed in the literature* we need not dwell upon the formal 
steps necessary to prove (34); we have only to find out for what range of y this particular 
combination of limits and functions is chosen correctly. From the asymptotic formula 
(33), it follows that the integral exists as long as y is real and not negative. However, 
in the formal steps necessary to prove that the integral obeys equation (4), factors such 
as cosh y-cos wu appear which demand that the exponential produce convergence; hence 
the integral (34) defines some solution of (4) only as long as y > 0. Our next observation 
is that the solution is always the same regardless of the choice of me u. This is so because, 


for symmetry reasons, we have that 
| exp [—2k sinh y cos u]lse u du = 0. 


Thus, only one linearly independent solution of (1) is left in (34), producing always the 
same solution of (4) on the right hand side. That this function is just ke y is seen by 
evaluating (34) for large y by the saddle point method. The saddle point is at the origin, 
which permits us to write 


exp [—2k sinh y cos 7v] ~ exp [—2k sinh y(1 + 32°)]. 


Hence 


exp [—2k sinh y cos u]me u du 
v—i@o 


a+@ 


~ i me(0) exp [—2k sinh y] | exp [—k uv" sinh y] dv 


— exp [—2ks sinh 2 
(k sinh y) 


This asymptotic behavior is the same as the one of the definition (7) and hence equation 
(34) is proved for positive y. 

We now use equation (34) to continue the function ke y analytically outside this 
original range of definition. We start out by permitting values of y slightly off the positive 
real axis. As long as this deviation x is less than 7/2 this goes without difficulty; for the 
convergence producing factor in (34) is changed from its previous value to 


exp [—2k sinh y cos x cos u] 
which will produce convergence as long as cos x > O and sinh y > 0. Using the definition 


(29) we arrive thus at the equation 


atria 


| exp [—2ik cos x cos uJmeu du = A he x (35a) 


which holds as long as 
g(x) > 0, (35b) 
0 < R(x) < m. (35¢e) 


tN .W. McLachlan. Theory and Application of Mathieu Functions. Clarendon Press, 1947, Chapter X. 
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This range of definition can be extended further if we allow the path of integration 
to be deformed. The imaginary axis along which the above integral is running is one 
of a series of valleys along which it could be taken to produce a convergent result. The 
other valleys are separated from this one by distances 27, 42, etc., as shown in Fig. 6. 
These valleys have a width x, as shown by (35c); they are separated by ridges of the 
same width for which the integral diverges. Now as ®(x) changes the valleys of exp 
[—2ik cos x cos u] = K(2, u) shift. Set z + 2 + iy, u— u + w then we get for the real 
part of the exponent in K(z, u) 


exp [—2k (cos x cosh y sin wu sinh v + sin z sinh y cos u cosh v)] 





FIG. 6 PERSPECTIVE DRAWING IN THE U-PLANE SHOWING THE 
ABSOLUTE VALUE OF THE INTEGRAND AND THE PATH OF INTEGRATION 
FOR THE INTEGRAL (35). THE FACTOR ME u IS NOT INCLUDED AND 
x MUST LIE ON THE WAVY LINE NOT FAR FROM THE REAL AXiS 


which, for v positive and large, becomes approximately 
exp [—ke" (cos x cosh y sin u + sin x sinh y cos u)). 


Clearly, if y is also large, the exponent contains sin (x + u) which is kept at its maximum 
value by setting 


r+u=-. (36a) 
2 
If y is not large then the bottom of the valley is given by 
tanh y tan x tanu = 1. (36b) 


The movement of x as a function of u does not differ essentially between (36a) and 
(36b). As @(x) decreases from 2/2, R(u) increases by an amount which is essentially 
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The movement becomes gradually more jerky as y becomes 


ceases to function for y = 0. Similarly for v negative and large we find 


equal to this decrease. 
nd 


smaller and 


for y large 


7 oF 
(3/8) 


and for general y 
tanh y tan xtanu = —1. (37b) 


The movement of the valleys is thus in the opposite directions on the two sides of the 


real axis. This means that for decreasing ®(2) the path of the integral (35) deforms 


as shown in Fig. 7. 


























FIG.7 DEFORMATION OF THE PATH IN THE U-PLANE 
FOR THE INTEGRAL (35). SUCCESSIVE NUMBERS ARE FOR 


R(x) DECREASING, STARTING FROM Z. 


We have thus found, for the entire half plane for which g(x) > 0, the solution of 
as, along the line parallel to the imaginary axis and passing through 7/2, 
- the solution is given to us in the form of an 


(1) which h 
the behavior prescribed by (29) and (7 
integral 


— Dik: COS X COS U | me uU du 3a) 
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vhere 
u, obeys (37), (38b) 
Us» obeys (36), (38e) 
d(x) > 0. (38d) 


Fig. 8 illustrates the path for the special case R(x) = 0. 

Formula (38) will now be used to find the range to which the asymptotic formula (7) 
is applicable. If we decrease ®(x) from 2/2 to some smaller value the two ends of the 
path displace as shown in Fig. 7; this displacement leaves undisturbed the location u = 0 
of the saddle point while at the same rotating its orientation clockwise; this rotation 





FIG.8 PATH FOLLOWED IN THE U-PLANE BY THE’ INTEGRAL (38); 
THE CASE SHOWN IS WHEN x LIES ON THE POSITIVE IMAGINARY 
AXIS: THE SADDLE POINT OF THE INTEGRAL LIES AT U=0. 


shows up in the passage from Fig. 6 over Fig. 8 to Fig. 9. For large 9(x), the angle a 
of the saddle with the z axis is found to be 


a=-~+=R(2) (39) 


The integral (38) yields then 


exp [—2ik cos x — in/4] 
* yi72 ~ 


(40a) 


i¢ 


(2k cos x 
ith the restriction 


I(x) > 0 (40b) 
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ray 
@ 


However, when &(x) drops below —7/2 the situation is altered. The two valleys are 
now so far removed that the path has to proceed over three saddles, lying at —7, 0, x. 
This situation is illustrated in Fig. 9. Of the three contributions, the one at the origin 
retains its analytical form (40a) but the two others will add to it and thus invalidate 
it. A similar modification must be applied when ®(x) increases from +72/2 beyond 
+ 3/2x. The asymptotic expansion (40a) is therefore valid within the range 


3 
R(x) S57. (40c) 


ES) 
IIA 





THE U-PLANE BY THE_INTEGRAL (38) 
ER R(x) HAS DROPPED SOMEWHAT BELOW -Z. THE ORIGINAL 
POINT REMAINS AT THE ORIGIN, ROTATING CLOCKWISE; IN 


FIG.9 PATH FOLLOWED IN 


rs 
ON, TWO NEW SADOLE POINTS, AT -77 ANO +77 MAKE 
A 


PPEARANCE. 


In the same way we obtain from equation (30) 





. exp [+27k cos x + in/4] 
hens (2k cos x)'”” (41a) 
if 
3. < az < +! (411 
97 HDS 9™ )) 
I(x) > 0, (41c) 
and from (31) 
hee ~ SP [= 2th cos x — tx/4) (42a) 
(2k cos x) 
if 
3, <= R(x) Ss +5 TT, (42b) 


d(x) K 0, (42c) 
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and finally from (32) 


(*) _exp [+27k cos x + tx/4] 
. (2k cos x)'”” (48a) 








he 


r 5 G(2) < +35, (43b) 


I(x) K 0. (43 ¢) 


This defining range gives to (40) and (41) a common domain of existence and also 
to (42) and (43). In order to give a common range to three functions we use equation 
(9). We get from it and (29) that, for positive large y 


a(t _ .\ _ exp [2k cosh y + #] 
” ( iv (2k cosh y)'” 





while from (42) and (43) 


a(n. exp [2k cosh y — ix/2] 
he ( iv) (2k cosh y)\? ’ 





(4) T =e exp [—2k cosh y] 
” (; iv) (2k cosh y)'” 





We know in addition that a universally valid linear relation must exist between these 
three functions. In the range considered, the first two are asymptotically large, while 
the third is asymptotically small. Hence, the relation must read 


hea = ie* hex + (unknown factor) - he z. (44) 


Equation (44) is half a connection formula; we shall now see that the other half can be 
picked up from consideration of Floquet’s theorem. 

5. The Floquet Function. According to Floquet’s theorem, there exists at least one 
solution of (1) which is multiplied with a constant factor when we apply the transla- 


tional symmetry operation of equation (1) 
Z—->Z+ fF. 


This constant is usually written in the form e***. It is known that e‘** oscillates back 
and forth from the unit circle to the real axis, having alternately positive and negative 
sign in the latter case. Whenever e**® is real we shall use the supplementary definition 


ev = +e’? (45) 


the sign being chosen so as to make the real quantity e"’ positive. The lines along 
which e‘*® = +1 have been studied intensively; a reproduction of the published results 
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is shown in Fig. 10.°. Figure 11 shows a reproduction of a map of 6 published 
by McLachlan. This type of information is still rather incomplete at this time. 
We now define as je x or more specifically je* 


obeys 
e (x = "je" x. (46) 


x the solution of equation (1) which 


A second, generally independent, solution je~ x is then obtained through 
(47) 


which vields the identity 
js 18) 


will be left for later disposal. A very important property 


The normalization constant 
je also equals ce,x, and similarly for se,x 


ae equals CE, Lt, 


of je x is that wheneve! 
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w 
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" FOR THE MATHIEU FUNCTIONS 
IN SHOW THE COMBINATION OF 
VALUES YIELDING PERIODIC SOLUTIONS. INFORMATION CONCERN- 


THE TWO PHASE ANGLES £ AND 7 'IS ADDED. 


FIG.10 “BUTTERFLY DIAGRAM 
@-q PLANE. THE LINES 


THE 


ING 


ir instance 8. Goldstein, Trans. Camb. Phil. Soc. 23, 303 (1927); E. L. Ince. Proc. Roy. Sox 
Also the textbook of McLachlan. loc. cit. 


‘Soo f 


Edinburgh 46, 20 (1925); 46, 316 (1926); 47, 294 (1927). 
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FIG. 11 DIAGRAM SHOWING CURVES OF CONSTANT # IN THE 
@-Q PLANE (McLACHLAN). INFORMATION OF THIS TYPE IS 
MUCH MORE RESTRICTED THAN THE TYPE SHOWN IN FIG.10 


The asymptotic behavior of the Floquet function was first obtained by Dougall.* 
We shall derive it more quickly by our integral equation method developed in the last 
section. The pair of equations to be established is 


B’ je" z= | exp [—2zk cos x cos ulje u du, (49a) 
B je x= | exp [— 27k cos x cos u]je” u du. (49b) 
6 


“> is defined by (38c). The resultant path is shown in Fig. 12 for the case R(x) = 0. 
5J. Dougall, Proc. Edinburgh Math. Soc., 34, 4 (1916); 41, 26 (1923); 44, 57 (1926). 
°The path is identical with the one for the generalized Bessel integral for Bessel functions. The 
integral permits thus derivation of the Bessel expansion by substituting the Fourier expansion under 


the integral sign. These relationships suggested to the author the symbol je for that type of function. 
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It is obvious from the discussion preceding equation (38) that as long as g(x) > 0 
the two integrals (49) exist and define two solutions of Mathieu’s equation. To complete 
the proof we have to show only that they also obey Floquet’s theorem. This is seen 
as follows. Both terminals of the path lie on the positive imaginary side of the real axis; 
therefore if ®(x) increases from an initial value, 0 say, the valleys in which the path 
terminates move in the sense contrary to 2, in accordance with equation (36). When x 
has been increased by 27 the path has been shifted without distortion by an amount 
. This shift leaves the kernel K(x, u) of the integral equation invariant, but multi- 
plies the factor je” wu in the integrand of (48a) with e**'’; hence the integral has been 
multiplied with this same factor and thus obeys Floquet’s theorem with the factor 
*‘® The same procedure establishes (49b). B* and B™ are constants which will be 


ou Dap 
milk 


? 
discussed below. 

The asymptotic expansion of je” x and je 2 is obtained from (48) by applying 
the saddle point method discussed in the previous section. The two saddle points lie 
at —7z and 0, as shown on Fig. 12. Each saddle point furnishes one of the exponentials 





FIG. 12 PATH OF INTEGRATION FOR THE INTEGRALS (49) 
WHEN x LIES ON THE POSITIVE IMAGINARY AXIS 


(33); they can be combined into a single term because we obtain the relative magnitude 
of the terms from (46) or (48). The resultant expressions are thus found to be 


1/2 l " 
B* je” x ~ 2 je 0 exp (4778) ( - cos (24 cos & — : + 1 13), (50a) 


k cos x 


1/2 | e 
B je” x ~ 2 je 0 exp (— 4x8) ( * ) cos (2 cos x — ; - 1 78), (50b) 


k cos x 


with the restriction 
I(x) > 0. (50c) 
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In the real direction the formulas (50) are also limited. The domain of validity is ob- 
tained by the methods used to establish (40c). We find 

51S R(x) S +7. (50d) 
We have, however, in this case the exceptionally favorable situation that the formulas 
(46) and (48) just supplement (50d) so as to furnish the asymptotic expansion of je x 
for all values of ®(x). One simple way to express this is by saying that under the con- 
dition (5Uc) alone je* x is asymptotically equal to J_s(ke"‘*), and je” 2 equals J..(ke™'*). 

The equations (50) would yield the connection formulas between the functions je 
and he were it not for the two undetermined constants B* and B~. This indeterminacy 
can be partially removed by inspection, as follows, 

(a) As the B’s in (50) are factors in the asymptotic expansion of the same function 
in two different regions of the complex plane, their ratio can never be zero or infinite. 

(b) When je x « se,x both B’s are zero; this follows from symmetry considerations 
on (49). The formulas (50) do not lose their meaning however, because je 0 also vanishes. 
In these equations symmetry demands that the ratio B°/B™ approach —1 as je ap- 
proaches se,, . 

(c) When je x « ce,x, the B’s in equation (50) cannot vanish, because ce 0 does 
not vanish. Symmetry or the equations (49) demands that the B’s be equal. This means 
B’/B” = +1. 

(d) Suppose now we are in a region where e‘*” is real. We then enter into equation 
(1) with the substitution suggested by (45) 

mex  e p(x). 
The resultant equation in p(x) is real; its periodic solution must be real because other- 


wise there would be two linearly independent periodic solutions. Hence je z is a real 
function of x, and for real je 0, the two expressions (50) must be conjugate complex. 


This gives 


BR 


a eo (51a) 


where 6 is given by (45) and y is some real number which is integer at the two limiting 
lines of the region and increases (or decreases) by 1 as we proceed from the line ce, z 


to 8é,, Z. 
(e) Now let e'”® be on the unit circle. We then enter into (3) with the substitution 


mex = e**p(z). 


The resultant equation in z is real. If there is to be only one periodic solution p(z) it 
must be real along the z direction, and hence the asymptotic expansions (50) must be 


real. It follows that 
BY = pag? tore (51b) 


where c is some real number which vanishes for integer 6. The undetermined sign is 
fixed in each of the separate regions of real 8 being + between ce,, and séz,,, and — 


between Céo,4; and sé2, . 





f) W S Dp S ormation by writing 
B 
2) 
J 
where 6 . ¢ ) I 1 I h are hxed p to 1 even integ i the 
sig Che ( he boul \ es ol ig. 10 is shown in Ta [. These 
lines divid he 
| f ) l l! Ss kind 
ven 
id 
ld 
a-q-plaie into ‘wi ‘ l “oaps”’. Betwer these two there Is a reciprocal ehavio1 
of 8B and y. In the 1 s, 8 is real and changes by l; y does the same thing in the gaps. 
Inversely, the real part of 8 is fixed in the gaps, and it has an imaginary part vhich 


varies; this behavior is d iplic: ted by y in the wings, where it has a variable 1 winary 
part zc introduce 51b). This variation is exhibited in Fig. 10. 

6. The Connection Formulas. We start out by writing a formal connection mula 
with the help of the parameters 8 and y of the last section. We dispose of the normaliza- 


tion factor by setting in accordance with (52) 
> 1 1 ] 
B~ exp [—3ir(8 + y)] B™ exp [4ir(8 + y)] — 
e 0 = >) 
9/9. - 9 ar 


This reduces (50) te 
1 


exp (— 477} ( -, © 3 . 
E : ma cos | 2k cosx — - xB), ota) 
2k Cos 2 } 


1 
ener eran . 1 : 
r~ a = “ cos (ox cos xz — =e x8). 54b) 


2/ COS I } a 


The domain of validity is given by (50): it is contained in the larger domains (40) and 
(41). We may therefore write down the connection formulas 


et Li exp (—4imry) [exp (—4irB)he' a2 — exp (+3i7B)he~' x], (55) 


je x = 4 exp (+3iry) [exp (+4i28)he'x — exp (—4i1B)he "gi. (56) 


Fu-her, reversing the sign of x with (29), (30), (31), (32) and (47) 


‘4 1, exp (S017) exp simB he’”’s — exp (— SimB)he . rl. (57) 


. . ° 4 -O 
ex = 31 exp (—3iry) [exp (—3irB)he”’'x — exp (ZimB)he”'z}. (58) 


‘ : » . 1 
Solving (55) and (56) for he’’’x we get 


hex = cosec rBlexp {4ix(y — B)}jex — exp {—}inly — B)}je-z]. (59) 
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Eliminating je* 2 and je x from (57), (58) and (59) we get a relation between he” x, 
he x and he'® x which must be identical with the previously derived incomplete equa- 


tion (44). Comparing coefficients, we find thus 


.sin ry : 


ee =e, (60) 
sim 7p 
and the mpleted connection formula (44) 
i 1 7 > 
he''x =tehe” x — (te cos 7B + COs Try)he oe. (O1) 
With equation (60), the formal connection formulas (55)-(59) become actual ones, 
th coefficients expressible in terms of 3 and ®, which, in turn, are known functions 
of a and y. By the use of the subsidiary definitions (45) and (51) we can show up (60) 
as an equation between real quantities. Inside the wings of Fig. 10 we get 
ae = 
sinh rw =e |sin x6 |, (62a) 
sinh rb = e€ | sinzy |. (62b) 
Actual 0) contains a little more than (62); in the wings for instance, it tells us for 
each one of the two 6’s on which side of the real axis the corresponding y is to be found; 
similar sign is determined in the gaps. 

Equa 61) effectively terminates the search for continuation formulas because 
tells us how to continue the simple exponential asymptotic behavior (40) on the other 
le of the real axis. A large number of other formulas are derivable from the ones 

; ed. auch 
— : : 1) , . { (es 
exp | Lime — )|exp (Stry)he Z + Cxp \ — dimy)he zi, (63) 
, , «6 - : eS 
] x= —tehe?« + (ie cos 78 — cos ry)he”' 2, (64) 
ad so torth 


We now come to the connection formulas for the Lindemann-Stieltjes functions. This 
sk has two stages of difficulty. By symmetry alone we can write down relations such as 


cexr« je x aa je Be 


cex x he Ye + he x, 


dex « exp (—}ixB)je*x + exp (+4178) je” 
dex « he’ a+ he x, 


(4 
tex « hex — hex. 


In order to make connection formulas out of these proportions it is necessary to have 
quantitative information for the points z = 0 and x = 7/2. This paper contains no 
such information for the point x = 0. For the point x = 2/2, on the other hand, the 
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necessary results were obtained incidentally in section 2. We will pursue this only to 
the zero stage of the Jeffreys approximation in which we get from (13) 





1h 
ke 0 = (a a » 3 avi (65) 
ke'0 = —exp (44,)(a + 2k’)'. (66) 

Hence we may write more precisely 

] (4) 1) a 
dex = > kk G (he x + hex), (67) 
t a) 4) 2O 
lex = 2 | ke’0 | (he ’x — hex). (68) 


The connection formulas which exhibit the asymptotic properties of the Lindemann- 
Stieltjes functions follow from the ones above by application of (61). We find 


cex « (1+ ie*)he'?x — (ie® cos 78 + cos ry)he z, (69) 
sex « (1 — ie*)he x + (cos ry + ie” cos rB)he” zx. (70) 
From (40) and (41) it is evident that these two functions always are of the form 


i ay on ae 71s 
a a 5 COS (24 cos zr — | + 1) (71a) 


with 
oxy _ second coefficient (71b) 


—F = — OE 


first coefficient 


Perusal of (60) shows that the y so defined is always real (not alternating as in (54)). 
Simple phase shifts of 0 or 7/2 result from (71) when 


cos 78 = +cos ry = +1 


in the combination circumstances warranted by Table I; the shifts are then identical 
with the ones in (54). The remaining relation of the pair (69) and (70) gives us then 
the asymptotic expansion of the Mathieu function of the second kind whose phase 


shift comes out to be given by 
tan y = (—)7"eO""®, 
The same formula (61) yields for (67) and (68) 


] 


2 ke 0 [—ie*he' x + (1 — cos ry + ie” cos rA)he zx], (73) 


dex = 


ez= 9 "es Ti) 7 [i e*he 2 + (1 + cos ry — ie® cos rf)hez]. (74) 


These equations yield a real phase shift only in the special circumstances warranted 
by Table I and equation (54). 
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The special information (65) and (66) regarding the point x = 2/2 also produces 
new information about the Floquet function at this point. Using (63), we get 


T Bs 1 = 
i: —~ exp (3 imp — w) he 0- cos 5) TY; (75a) 
cot rc 4 <a = 
je’ 5 = —exp |, inB — )-| ke’0 ‘sin 5 TY. (75b) 


From (75), the connection formula between je, de and te is readily derived. 

7. Concluding Remarks. This study is based on the notion that the Floquet param- 
eter 6 is a known function of a and q. That this is partly a convenient fiction is seen from 
Fig. 11. It is a surprise that the results of this paper do furnish some new information 
concerning 8. When we traverse the gap between wings from one bounding curve to 
the other, the exponential damping constant b is related to y by (62b). On such a path 
y changes from one integer to the next and | sin y | passes therefore through its maximum 
once. We get therefore the relation for such a path 


Max [e® sinh xb] = 1. (76) 
This relation was checked from the Figs. 1, 2 and 11 for a stretch where both quantities 


are known. The result is Fig. 13 which confirms the prediction (76). In general, e® varies 
sufficiently slowly so that (76) can be used to determine a rough upper limit for b. 
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FIG.13 PLOT OF e® SINH 77b AGAINST & FOR Q=0.65. 
ACCORDING TO EQUATION (76) THIS PRODUCT 1S BOUNDED 
BY THE VALUE 1 WHICH IS ACTUALLY ATTAINED AT SOME 
POINT. THE CRUDENESS OF THE GRAPH REFLECTS THE 
LIMITED INFORMATION OF FIG.11. 


Notwithstanding this small bit of supplementary information concerning 8, there 
remains the fact that the two natural parameters for the connection formulas are 6 
and y, and that their dependence on the equation parameters a and gq is erratic and 
not expressible in closed form. It is not likely at this stage that an analytic relation will 
ever be found connecting 8 and y to a and q. For this reason, this paper proceeds instead 
to find slowly varying and easily determined functions of a and qg, from which 8 and y 
can be determined by analytic means. One such parameter is , whose determination 
was carried out in section 2, and for which the analytic relation to 8 and y is equation 
(60). It is obvious that the accomplishment of the program calls for another such 
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qr 
i? 2) 


parameter, to replace 6. This parameter has not yet been found. It is interesting to 
note as a possibility that 8 + y and 6 — y are simpler in their behavior than either 
one of them alone. In the meantime, the connection formulas of this paper must be 
used in conjunction with whatever published information is available concerning 8. 

8. Appendix on Integral Equations. The following integral equations are new, to 


my knowledge, but proved to be of no use in deriving connection formulas. They may, 


however, be useful in the hands of others 

Let z, ¢, w be complex numbers whose real parts are x, &, uw and whose imaginary 
parts are y, 7, v. Then the formulas are 
B jez je& = je 0 | exp [—27k cos z cos ¢ cosw — 2ksinzsin ¢sin wljewdw, (77a) 
B*je*z je § = je 0 | exp [—27k cos z cos § cos w — 2ksinzsin {sin w]jewdw. (77b) 
The formulas are generalizations of (49) to which they reduce for the case ¢ = 0. The 


path and the range are best discussed in two stages. If 
= (), (77¢c) 
then we need 


y—n> 0, (77¢’) 


and the path is exactly the one shown in Fig. 12. If 


x — é arbitrary, (77d) 
then we need 
2 4 rT, 
y—n> QO, cosh (y — n) > 3? (77d’) 
and the abscissa u of the terminal valleys is given by the generalization of (36b 
tan (x — &) tanh (y — n) tanu = 1. iza’*) 


We set out immediately to prove d and will get c as a special case. We follow earlier 
proofs quite closely. The formal equivalence of the integral to a product of solutions 
of the Mathieu equation is found in the literature.* Although the integral exists when 
the exponents cancel, we need a negative real part in the exponent in order to implement 
the formal steps. This exponent reads 


exp 2ik cos z cos ¢ cos w — 2k sinzsin ¢ sin w + 2k cos w. 


The last term arises from the contribution (33) of the Floquet function; the difficult 


{ 


sign is the positive one. By an obvious transformation this becomes 


exp ike cos (2 — ¢) — tke COB (2 “Tr () T+ 2k cos w. 
The second term is of no importance because v is to be positive and large. Introducing 


+ 


real and imaginary parts this becomes 


exp = tke’ e eos (x — £) cosh (y — n) — isin (a — &) sinh (y — n)] + the 


Ss 
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and the exponent’s real part is 


R(ex] —ke" {cos (x — &) cosh (y — 7) sin u 
(77d’’) 
+ sin (2 — £) sinh (y — 7) cosu — sin u]. 
We now introduce the choice of u indicated by (77d”’). This means 
, cos (x — &) cosh (y — 7) 
sinu = ———— - ot 
(cosh” (y — n) + cos (x — §&) — teed 
sin (2 — &) sinh (y — 7n) 
COS = eo a ee ae _ — ac <pamias 756 
(cosh” (y — n) + cos’ (x — &) — 1)'” 
We thus get 
R(exp) = —ke'| (cosh” (y — n) + cos (rx — &) — i 
cosh (y — ») cos(x — §) 
(cosh” (y—t+ cos’ (x — &) — | a 
—ke’ . 
= _ - —;———. ; [eosh* (y — 7) 
(cosh” (y — n) + cos” (a — &) — yi [ (y * 
+ cos’ (xr — &) — 1 — cosh (y — n) cos (x — &)]. 
As stated in (77c), the curly bracket is positive when cos (x — £) = 1 (and also for 


; in the general case we form the perfect square 


cos (x — &) = 0); 
[eos (x — &) — 4 cosh (y — nf 


and then pull through on the remainder with (77d’). My surmise is that a more thought- 
out estimate could prove (77c’) all the time. Having proved the character of the function 
these restrictions we can determine the particular nature of the left hand side 
of (77) by first making ¢ = 0 and z large and positive, to get the function of z; and then 
Z 0 and ¢ large and negative to get the function of ¢. The formulas (77) are thus es- 
tablished. If we reverse (77c’) and set instead y — » < 0 the roles of y and 7 are reversed 


unde} 


because the integral is formally symmetric in z and ¢. 
What makes the integral (77) interesting is that one can pass with it from positive 
ive imaginary values provided (77c’) is maintained. The difficulty in getting 


to negatl 


asymptotic expressions is in the location of the saddle points; one saddle point of (77a), 
for instance, lies at w = ¢; this leads to a trivial cancellation and a confirmation of 
50b). The fact. remains, nevertheless, that formula (77) permits us to cross 


formul: 
| axis; this the simpler equations (38) and (49) do not permit us to do. 


the re 
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MIXED BOUNDARY VALUE PROBLEMS IN SOIL MECHANICS’ 


By 
R. T. SHIELD? 


Brown University 


Summary. The stress-strain law for an ideal soil formulated in a recent paper 
[1]* is applied here to obtain the velocity equations referred to the stress characteristic 
lines in plane strain problems. Simple velocity fields associated with families of straight 
characteristic lines are then examined, together with discontinuities in the velocity field. 
The results are applied to obtain the incipient velocity field for the indentation of a 
semi-infinite mass of material by a flat punch or footing, and to solve the problem of 
indentation by a lubricated wedge. 

1. Introduction. In deriving the solutions of two-dimensional problems in soil 
mechanics, it is usual to assume that the soil is a plastic material in which slip or yield- 
ing occurs when the stresses satisfy the Coulomb formula [2] 





ao + 


r, “—s" | ‘ ‘ 
jo-—> oy sin og+ V4 (o, — o,) + 7, —~c cose = 0, (1) 


where c is the cohesion and ¢ is the angle of internal friction of the soil. The stresses 


also satisfy the equilibrium equations 


Oo, OT ry 
= ———— = (), 
Ox OY 
f (2) 
OT cy 00, 


+= 

Ox Oy J 
in which the weight of the soil is neglected. The equations (1), (2) are hyperbolic and 
the two characteristic directions are inclined at an angle 7/4 + ¢/2 to the direction of 
the algebraically greater principal stress. In Fig. 1, the lines 1, 2 are the directions of the 


YA 


1/4 + $/2 
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FIG.!. Failure Lines 
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principal stresses o, , 2 (a; < o2) at a point and the lines a, B are the characteristic 


lines passing through the point. We shall call the characteristic line which lies between 
the 1 and 2 directions the first failure line and denote by 6 the angle of inclination of this 
line to the z-axis. 
If we now put 
p= rt —" > 0, 
f 2msing ~ 


where m is a positive constant which has the dimensions of stress, then, using (1), it 


can be shown [3] that 


o, = —mp|l + sin ¢ sin (26 + ¢)| +c cotg, 

y= —mai— singsin (20+ ¢ i + ¢ COL ¢g, (4) 
T mp sin ¢ cos (26 + ¢) 

7 N ¢ S « r ¢ 


The equations of equilibrium (2) can be replaced by the equations 


1 cot ¢ log P + &@ = const. along a first failure line, ) 
: D) 
2 COL YG log Pp ¢ = const. along a second failure line, 


which were first obtained by Kétter [4]. 
The three stress components, o, , c, and r,, , are determined from the equilibrium 
equations (2) and the yield condition (1); alternatively, they can be determined by 


integrating equations (5) along the failure lines. The usual treatment of plane-strain 
problems assumes that the problems are statically determinate. In general, however, 
the stress boundary conditions are not sufficient to make the problem statically de- 
terminate and a stress-strain law is necessary in order to allow a more complete in- 
vestigation of the problem. In the following we shall use a stress-strain law which is 
derived by assuming that the soil is a perfectly plastic body. It should be remarked that 

(such as the presence of water in the 





soil). The predictions of the theory must be compared with the actual behavior of soil 


in order to obtain an indication of the value of the assumption. 
2. The velocity field. Drucker and Prager considered a proper generalization 


of the Coulomb hypothesis (1) and they showed that if the soil is assumed to be a plastic 


material then, according to the concept of plastic potential [5], the stress-strain law 
for plane strain corresponding to the yi ld function (1) is 
oT r T —_— @¢ - ) 
€ \ \sin ¢ i ( 
Oo F 4 ; a - T ) 
Oo} " \o 2 | ¢, 
€ ‘ sin ¢ (5 ) 
Oo 2 ’ , 0 ~-- T | 
aft T 
“ \ A 
OT ‘ T O,) w+ T 
where ¢. , €, , ¥:, are the plastic strain rates and is a positive factor of proportionality 


which mav assume different values for different particles. Since we assume that there is 
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no deformation of the soil until plastic yielding occurs, the plastic strain rate is equal 


to the total strain rate, and we have 


fry 


Ou Ov Ou Ov 
é.=- = 


’ eo a E x? 
OX OY oY OX 


where u, v are the components of velocity along the 2x, y-axes. Using equations (4), 


equations (6) can be written 


3 Bh Bae a 
é. = 5 (sing — sin (26 + ¢)}, | 


c, = » {sing + sin (26 + 9}, | (7) 
| 
Vey = d cos (26 + ¢). } 
From the stress-strain relations (6) or (7), the rate of dilation is found to be 
é.- +e, = Asing > 0, (8) 


so that an important feature of the relations is that plastic deformation 
must be accompanied by an increase in volume if g ¥ 0. 
If we put 6 = O and 6 = — (xr/2 + ¢) in turn in the first equation of equations (7) 


we find that 
(2) (2) ee: (9) 
Ox @=<0 Ox O=~(4/2+¢) 


These equations express the fact that the rate of extension along the failure lines is 
zero. 

The components u, v of the velocity vector are to be determined from equations 
(7) when the pattern of the failure lines is known for a given plastic stress field. It can 
easily be shown that the characteristics of the velocities coincide with the characteristics 
of the stresses and it is more convenient to refer the velocity equations to the character- 
istic lines. We denote by v; and v, the orthogonal projections of the velocity vector at a 
point on the directions of the first and second failure lines passing through the point 
see Fig. 2). The signs of these velocity projections are chosen so that a counterclock- 
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FIG. 2. Velocity Projections 
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wise rotation of the positive direction along the first failure line through an angle of 
7/2 + @ transforms it into the positive direction along the second failure line. The 


velocity projections v, , v2 are related to the cartesian components u, v of the velocity 
by the equations 


v, = u cos 6+ vsin 8, v, = —usin(@+ ¢) +vcos(@+ ¢), 
, a ee (10) 
v, cos (6 + ¢) — v2. 81n 6 v, sin (6 + ¢) + Vs Cos 6 
==> " | => ~ 
COs ¢ cos ¢ 


The substitution of these equations into equations (9), which state that the rate 
of extension along the failure lines is zero, gives the equations of the velocity field referred 


to the characteristic lines, 


dv, — (v, tan g + v, sec ¢) dd = 0 along a first failure line, 
> (11) 
dv. + (v, sec ¢ + v2 tan ¢g) dé = 0 along a second failure line.) 


These equations can also be obtained as follows. In Fig. 3, A and B represent two neigh- 





FIG. 3. Inextensibility of failure lines. 


boring points on a first failure line separated by an infinitesimal distance ds, . The end 
points of the element AB move with the velocities indicated in the figure. The extension 
of the element may be considered to be caused by the following two circumstances: 
(i) the velocity along the element at A is smaller than the corresponding velocity at B 
by an amount dv, , and (ii) the normal component of the velocity at A and the normal 
component of the velocity at B are inclined at an angle dé. The corresponding rates of 


i 


extension are dv,/ds, and — (v, tan ¢ + vo sec ¢) dé/ds, , since the normal component 
of velo Itv at A is?, tal gO TT Vo SEC ¢g. The condition that the total rate of extension along 
a first failure line must be zero is therefore 

dv, — (v, tan g + v2 see ¢) dé = 0 


along a first failure line. This is the first of equations (11) and the second equation can 


be obtained in a similar manner. 
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Equations (11), together with the velocity boundary conditions and the condition 
that the dilatation must be positive, suffice to determine the velocity field when the 
failure lines are known. 

3. Families of straight failure lines. The simplest pattern of failure lines consists of 
two families of straight lines intersecting at an angle 7/2 + gy. As can be seen from 
equations (5), this pattern corresponds to a region of constant stress and it is usually 
called an active Rankine zone or a passive Rankine zone in the literature of soil me- 





_ oe 





FIG. 4. Region of Constant Stress. 


chanics. In Fig. 4, the x and y axes are taken along the directions of the minimum and 
maximum principal stresses respectively. Since @ is constant along the failure lines, 
equations (11) show that v, is constant along a first failure line and v, is constant along 
a second failure line. If we denote by d, and d, the distances of a current point from one 
of the first failure lines and one of the second failure lines, say the failure lines which 
pass through the origin, then we can write 


v, = f(d,), v, = g(d.), 


W here 


d. = x sin (: ae ¢) + y cos (* ss e), 


and where the functions f and g are such that the dilatation is positive everywhere in 
the region. The cartesian components of the velocity vector are found from (10) to be 


us (sla) cos (* + ¢) iin q(d2) sin (* a e\} / cs ‘ 
ee Oe a 
i= \ f(a) sim 4 + 9 + g(dz) cos 4 _— > cos ¢. 
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When one of the families of failure lines consists of concurrent straight lines, the 


other family is a system of logarithmic spirals which have the point of intersection as 
centre. This stress distribution is usually called a zone of radial shear. In Fig. 5 we have 


J 














FIG. 5 Zone of Radial Shear. 


taken the first failure lines to be the family of straight lines and we denote Dy Tr the 
C | . ] . 4 ] - _— = - 
aistance m the centre of the splrais. Since @ 18 constant along the first failure lines, 
1S ¢ : 110Nng se lines and Wwe have 
7 
Substitu ing’ fol ) 7 the second ot eauatl Ss 11) we obtain 
sec g | 6 | 
along a second failure line. In this equati 1 is a constant of integration which is con- 


1° 1 1 


stant along each second failure line but which may take different values on different 
second failure lines. Since the spirals are given by the equation 
re '* = constant, 


1e form 


the function A can be written in t 


The functions f and g must be such that the dilatation is positive everywhere in the 


regio! When f 8) is zero, i1.e., when v; is zero, we obtain 


r oe: A 4 


so that v, varies exponentially along each second failure line. 


If the second failure lines had been taken to be the family of concurrent straight 
lines then we would have obtained 
6 


A f(@), ve? = gee @ | F( Ode “'e? do + A, 


“6.6 
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where A is constant along each first failure line. As before, f(@) and A must be such 
that the dilatation is everywhere positive. 


4. Discontinuities in the velocity field. It can be shown, as in the theory for a 


perfectly plastic material, that a line separating a region of plastic flow from a region 
vhich remains at rest must be a failure line. This follows from equations (11) because, 
if the line were not a failure line, a region at rest on one side of the line would imply a 


certain region at rest on the other side of the line. This still applies if the velocity field 
discontinuous across the line separating the two regions. A line of discontinuity in 
the velocity field must be regarded as a thin layer producing a continuous transition 
from one velocity field to another. Since a discontinuity in the tangential velocity must 
be accompanied by a separation or discontinuity in the normal velocity, the transitional 
yer must have appreciable thickness for a soil while there is no need for such a layer 

1 a Prandtl-Reuss material (for which g = 0). 
Let P be a point on the median line of such a transitional layer and take the x, y axes 
long the tangential and normal directions at P. Since, at P, du/dx and dv dx must be 
negligible compared with du /dy, the strain rate e, must be very small compared with 


the strain rate y,, . The stress-strain relations (7) show that this can only be so if we 
have @ = O ora, or 6 = 7/2 — ¢ or 32/2 — g, that is, if the line of discontinuity is a 


first or second failure line respectively. Also from (7), for these values of 6, we have 


\ 
Ov . 
€, > A sm ¢, 
oY 
(12) 
Ou 
Vxx * | _ +X COS ¢, 
OY 


the positive or negative sign being taken according as the line of discontinuity is a first 
or second failure line. Equations (12) show that the change in velocity across the line is 
nclined at angle ¢ to the line of discontinuity. Further, if we denote by u, and u,. the 
values of the velocity component u on the sides of the line of smaller and larger values 
of y, then equations (12) show that vu. > u, or u, > u, according as the line of discon- 
tinuity is a first or second failure line respectively. 

The straight line and the logarithmic spiral of angle ¢ are the only lines of discon- 
tinuity which permit rigid motions, translation and rotation respectively, of the regions 
separated by the line. These two types of “sliding” discontinuity have been used by 
Drucker and Prager to obtain upper bounds for the critical height of a vertical bank of 
soil 

5. Indentation by flat punch. In this section the theory developed above is applied 
to the indentation of a semi-infinite mass of soil by a flat rigid punch or footing (under 
conditions of plane strain). A possible plastic stress distribution was determined by 
Prandtl [6] and we shall consider this stress distribution together with an alternative 
solution. The two stress distributions, which give the same value for the bearing capacity 
of the soil, are illustrated in Figs. 6 and 7. They correspond to the two solutions proposed 
by Prandtl and by Hill [7] for the same problem in a perfectly plastic material (for which 
¢ = 0). It seems probable that Prandtl’s solution is more nearly correct when the punch 
is sufficiently rough, while the other solution will apply when the surface of the punch 


is smooth. 
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Referring to Figs. 6 and 7, plastic regions begin to form at A and B as soon as the 
load is applied to the punch, but no indentation is possible until the plastic region extends 
all the way from A to B. We consider only the incipient plastic flow so that the boundary 
conditions are satisfied at the underformed surface. The problem of determining the 
stresses and velocities after the punch has penetrated a finite distance is of greater 
difficulty and would require a study of the successive phases of the plastic flow. 
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FIG. 6 Prandti Solution. 
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FIG. 7 Alternative Solution 


Consider first the solution represented by Fig. 7, and since it is symmetrical about 
the axis of the punch, we need discuss only the left-hand plastic region. The regions 
AOC, ADE are regions of constant stress (active and passive Rankine zones respectively), 
while the region ACD is a zone of radial shear. We take the downward velocity of the 
punch as the unit of velocity so that along AB the downward component of the velocity 
of the material must be unity. The material below the second failure line OCDE remains 
at rest so that OCDE is a line of discontinuity. Thus the velocity along this line is every- 


where inclined at an angle ¢ to the line, i.e., along this line v, = 0. It follows that », 
0 throughout the plastic region since v, must be constant along each first failure line. 

In the region AOC, v, is also constant and the region moves as a rigid body in the 
direction perpendicular to AC. If v is the velocity of the region then the boundary con- 


dition along AO gives 
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and we also have 
T 


Vv. = vcosg = cos ¢gsec (= + 


dNoIs 
ed 


in this region. In the zone of radial shear ACD, we have 


—@tane 
vv», = Ae tan 


where A is constant along each second failure line. Since v, is constant along AC, we 
see that v. is constant along each first failure line of ACD and also 


T p mainte 
v. = cosy see (* + £) Hi . 


along AD. Finally, the region ADE moves as a rigid body with velocity V = sec (x/4 + 
g/2)e"*’** in the direction perpendicular to AD. The velocity field is represented 


by the small arrows in Fig. 7. 
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FIG. 8 Resulting deformation of square grid if incipient velocity field was 


mointained for a short time, according to alternative solution 


Figure 8 shows the distortion of a square grid which would result if the initial velocity 
field was maintained for a short period of time (where ¢ is taken to be 20°). In obtaining 
this diagram a thin transition layer was assumed to exist between the line OCDE of 
Fig. 7 and the material which remains at rest. The initial position of the layer is indi- 
cated by broken lines in the figure. 

The initial velocity field for the Prandtl solution, Fig. 6, can be obtained in an anal- 
ogous manner. In this case the region ABC moves downward as a rigid body and the 
lines AC, BC in addition to the lines CDE, CFG, are lines of discontinuity. As before, 
we take the downward velocity of the punch as the unit of velocity. Referring to the 
left-hand side of Fig. 6, the material below the second failure line CDE remains at rest 
so that the velocity along this line must be inclined at an angle ¢ to the line. Hence 
v, = 0 along this line and since the first failure lines are straight it follows that v, = 0 
in the region ACDE A. The velocity of the material just to the left of the discontinuity 
line AC is perpendicular to AC and its magnitude must be such that the change in 
velocity across AC is inclined at an angle ¢ to AC. By drawing the velocity diagram 
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shown in Fig. 9, we see that the velocity must have the magnitude 
. {7 ¥ | 
= a \4 9 / 
~ r - 
In the zone of radial shear ACD the velocity increases exponentially to the value 
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FIG. 9 Velocity diagram for discontinuity line AC 


n Fig. & 


on the line AD. The region ADE moves as a rigid body in the direction perpendicular 
to AD with this velocity V. 

The small arrows in Fig. 6 represent the velocity field. The distortion of a square 
grid which would result if the material moved with this initial velocity field for a short 
period of time is shown in Fig. 10 (where ¢ is taken to be 20°). Thin transition layers 
were assumed to exist between the lines AC, BC and the material in ABC which moves 
downward with the punch, and a thin transition layer was taken between the line EDCFG 
and the material at rest. These layers are indicated by the broken lines in Fig. 10. 

We notice that the velocity V of the material along AF in the Prandtl solution is 





FIG. 10 Resulting deformation of square grid if incipient velocity field was 


maintained for a short time according to Prandtl solution 
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exactly half the value of the corresponding velocity in the alternative solution and the 
velocities are inclined at the same angle, 7/4 + ¢/2, to the surface AZ. For a given 
width AB of the punch, however, the length of AF in Fig. 6 is twice the length of AE 
in Fig. 7. It follows that the volume of soil raised above the undisturbed surface after 
the punch has penetrated a small distance is the same in both solutions. 

6. Wedge indentation. We now consider the problem of the indentation of a semi- 
infinite mass of soil by a smooth, rigid wedge under conditions of plane strain. The 
solution of the same problem in a perfectly plastic material was obtained by Hill, Lee 
r [8], and the solution for the material that we are considering follows this 


and mnne 
nd Tuppe 


closely. Since the configuration is geometrically similar at each stage of the 


penetration, it is possible to obtain the complete history of the motion without following 


SOLUTLON 


the deformation step by step. 
HAs the wedge is pressed into the soil, the displaced soil will form a raised lip at each 


side of the wedge, and the shape of the lips must be determined as part of the solution 


to the problem. We shall assume that the surfaces of the lips are straight and show that 
vhich satisfies this assumption. The pattern of failure lines is indicated 
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FIG. 11 Indentation by lubricated wedge 
in Fig. 11. AF is the right flank of the wedge, which is of angle 28, BG is the as yet unde- 
formed surface of the soil and AB is the lip. The regions AED, ABC are regions of con- 


stant stress while ADC is a zone of radial shear of angle a. We denote by / the length AB 
of the lip, by 6 the distance of O from AB, and by h the elevation of A above OB. The 
depth of penetration of the wedge is denoted by ¢ and if the downward velocity of the 
wedge is taken as the unit of velocity, we may take ¢ to be the time variable. 

The lip AB makes an angle 8 — a with the undisturbed level OB and it is easily 
shown that we have the following expressions for /, b, h in terms of 1, a, 6: 

lat J \ ) 
t/ se *""* tan (= - ¢) cos 8 — sin (8 — a)?, 


9 
/ / 


h {sin (6 — a) ti mais =: 
( f 
! 


° atan¢g T ° 
{sin (8B — ate , tan ( _ ©) sin B+ cos (gs — a} 


f os aa (= — ) cos 8B — sin (6 — a} | 


nos 


) cos 8 — sin (6 — a}, | (13) 
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The first failure line BCDE is a line of discontinuity in the velocity field so that the 
velocity vector along this line must make an angle ¢ with the line. It follows that v. = 0 
along BCDE and therefore v, is zero throughout the plastic region ABCDE, since the 
second failure lines are straight. In the region ADE, v, is constant and given by 


v, = vcCcOsg, 


where v is the magnitude of the velocity vector in ADE. The boundary condition along 
AE requires that the velocity of the wedge and that of the soil in contact with it must 
have the same projection on the normal to AE, and therefore 


P T Q 
_— cae | 
v= in B see (* + £), 


In the zone of radial shear the velocity increases exponentially along each first failure 
line, and along AC the velocity vector has the constant magnitude 


Prien cts {ft ¢)\ atang 
Vy = in B see (* + £) ; 
At a given instant, the region ABC is moving as a rigid body with velocity V in the 
direction perpendicular to AC. The velocity field in the plastic region is illustrated on 
the left of Fig. 11. 
The projection of the velocity of the lip AB on the normal to AB is 


\ 


A 


V cos (= -- é) = sin @ tan (: — he singh 


while the projection of the velocity of the vertex FE, which is moving downward with 
unit velocity, on the normal to AB is cos(8 — a). Hence, the distance of E from AB 
increases at the (constant) rate which is the sum of these two projections. At a time 
t, i.e., since the beginning of the indentation, the distance of E from AB has therefore 


reached the value 


Tr , ¢ \ xtan¢g 


tisin 6 tan (: + 5 | + cos(B — a)|]. (14) 


From Fig. 11 we see that this distance is also equal to 
b + t cos (8 — a), (15) 


and equating the expressions (14) and (15) gives 


. T Y)\ atang ; 
b = tsin 6 tan (" + 4) ‘ (16) 


The substitution of the expressions (13) into this equation furnishes a relation between 
the angles a, 8, and after some reduction we obtain 


. + 
atan¢ (7, ), rtang Tr _¢ 
COS al ¢€ tan (5 Tone tan \ eat 


aan SERENE ieee (17) 





Qq= 
f ; atan¢g Tv g ~atang T 9 
2sina + e**"* tan (? +- a ¥ 2? tan ( om: )} 
2 4° 2/J 


! 


cos (28 — 
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The variation of the angle a with the angle 8, obtained from (17), is shown in Fig. 12 
for g = 20°. 

The pressure P on the flank of the wedge can easily be obtained from the pattern 
of the failure lines and the boundary condition of zero traction along AB. It is found 


to be given by 
P =c cot A 2atang tan? (z ¢) —— 
died 49 


where, as before, c is the cohesion of the soil. 
90’ r? . 
a lec 
J +7 
60 4 +6 


4 Proc rs 











FIG. 12 Variation of P and a with B 


The variation of P with the angle 8 is shown in Fig. 12 for g = 20°. Since we have 


: . T ee g ,~ atan¢g 
AE = I tan (* e), : 


the total downward force F necessary to drive the wedge into the soil is given by 


eS ee wT €)\ -atang 
F = 2Plsin 8 tan (: 4) . 

The determination of the motion of a particular element of the soil appears to be 
rather complicated since the velocity of the element is influenced by its varying position 
in space and also by the continual expansion of the velocity field. However, the use of 
the unit diagram introduced by Hill, Lee and Tupper greatly facilitates the solution. 
The unit diagram is obtained by transforming the velocity field into a geometrically 
similar field in which the penetration is always unity. This is effected by the trans- 
formation 

r= r*t, 
where r is the actual position vector of the element with reference to O and 1™* is its corre- 
sponding position vector in the unit diagram. The method of determining the trajectory 
of an element is very similar to that used by Hill, Lee and Tupper for the perfectly 
plastic material and for brevity we omit the details of the solution. The unit diagram, 
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with three typical trajectories, is given in Fig. 13, and in Fig. 14 we show the deformation 
of the originally square meshes of a grid, where we have taken 8 = 30°, g = 20°. 

It is also possible to solve the problem when there is a moderate amount of sliding 
friction between the soil and the wedge. The angle of the wedge must not be too large 
in order to make it possible to extend the plastic field in a satisfactory manner below 
In this case the lines of failure do not meet the flank of the 


the vertex of the wedge 
¢/2 but the pressure on the wedge is uniformly distributed. 


wedge at an angle 7/4 
The lip AC is still straight although it is inclined at a smaller angle to the undisturbed 


surface. 
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INFINITE MATRICES ASSOCIATED WITH DIFFRACTION BY AN APERTURE* 


By 
WILHELM MAGNUS 


New York University 


1. Introduction and summary. As an example of their ‘variational method”’, 
LEVINE and SCHWINGER [1] investigated a boundary value problem which arises 
from the diffraction of a plane scalar (acoustical) wave by a plane screen with a circular 
aperture. It is equivalent to the problem of finding the field of a freely vibrating circular 
disk. A full discussion of the physical problems was given by Bouwkamp [2]. Let 2z, 


p, 9 be cylindrical coordinates and let z = O be the plane occupied by the screen. Let 
Z 0,0 S p < a define the aperture (or the vibrating disk). The diffracted field is 
given by a function u which satisfies V°u + k°u = 0 (with a constant k) everywhere 
except for z = 0 and at infinity satisfies a Sommerfeld radiation condition. For z = 0, 
u must satisfy the “mixed” boundary conditions u = 0 for p > a and du/dz = vp, with 
a given constant value v) for 0 S p < a. These conditions determine u uniquely. For 
Z 0,0 S p < a,u = (p) becomes a function of p only, and if ®(p) is known or even 


if only C,@(p) with an undetermined constant factor C, is known, u can be determined 
everywhere; see formulas (A.1), (A.2), (A.3) in [1]. 

Levine and Schwinger [1] show that the ratio of the energy transmitted through 
the aperture to the energy incident on the aperture is the imaginary part of the complex 
transmission coefficient 7*, which is a quotient of two integrals involving ®({p) quad- 
ratically. As a functional of ®(p), T* becomes stationary for the correct function ® 
which determines uv. Levine and Schwinger find approximate values for 7* by expanding 
first @(p) in an infinite series of auxiliary functions (see 3.1 and 3.2) with coefficients 
D Then 7* becomes a linear form in the D,, (see 3.10), and the unknowns D,, are 
determined by an inhomogeneous system of infinitely many linear equations with a 
coefficient matrix L (see 3.4, 3.5). In [1], these equations are solved ‘‘section wise’, using 
the first ] = 1, 2, 3, --- equations to determine the first / unknowns. All quantities D,, , 
7*, L are power series in 8 = ka/2, and Levine and Schwinger compute the first coeffi- 
cients of the expansion of 7* in a power series in 8 which were determined independently 
by Bouwkamp [2], who used spheroidal wave functions. 

It will be shown that the algebraic properties of the matrix L make it possible not 
only to find approximate values for 7* as in [1] but also to determine #(p). This is due 
to the fact that L factorizes in a product L“’ S,where L“ is the matrix for the static 
case k = 0 and where S can be inverted by solving finite recurrence relations. The 
details are stated in lemma 1 and theorem 1 of section 3. Lemma 2 gives additional 
algebraic relations. Problems of convergence and uniqueness are settled in section 5. 
These depend largely on an investigation of the properties of L“’ which is carried through 
in section 4. There it is shown that in the limiting cases k = 0 and k = ~ the matrices L‘” 
and L‘*’ of the linear equations also arise from a problem of moments. This also makes 
it possible to prove that the variational method for the calculation of the transmission 

*Received March 27, 1952. This work was performed at Washington Square College of Arts and 
Science, New York University, and was supported in part by Contract No. AF-19(122)-42 with the 
United States Air Force through sponsorship of Geophysical Research Division, Air Force Cambridge 
Research Center, Air Materiel Command. 
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~I 
¥ 8) 


coefficient will work even for / where the linear equations for the D,, do not have 


any solution at all 


2. Notations. I ements ol DELETE matrices are denoted by s ibs 5S wv, 
m= UV, 1, Z, lenotes the rows and m denotes the columns. A vect with 
components is de by {z,{. We also use the notations 

l I (a a | a l a = | 2.1) 


M 
bo 
bo 


where I’ denotes the ¢ ima function and F/’ denotes he hypergeometric series. Tor 
results needed here see Whittaker and Wats 3] and Bailey [4]. 


3. Algebraic properties of the linear equations. Let 


be the expal sion of tl field ® p) in the aperture in terms of powers of 1 — p /a Here 
C, denotes an undetermined constant and 


—lar, = D 3.2) 


where the D,, are the knowns used by Levine and Schwinger [1]. The linear equations 
for the z,, as obtained from the variational method can be written as follows: 
Let p, g = 0,1, 2 and let L‘*”’, L‘°*** be infinite matrices with elements 


122*®) defined by 


where, for any values m.t 

A(n, m, t) = I + 3/2)T(m + 3/2)T(n + m + 2? + JD), 

Bin, m, t) = 4T@+ )Tm+it )r(imt+it+ )Vrat+mt+ist 5/2). 
Let L be the matrix 


LD i 4 2 6°" PL 3.5) 


the general element /, ,, = /,.,,(8) of which is a power series in 8 = 3ka. Then 


bh Lt, =(n-+3 a. 3.6) 


0 


Let denote the vector with the components x and let —"’, r = 0, 1, --+ be the vector 


with the components z,,’ where 
Le = bs Bx”. (3.7) 


Let 7°” denote the vector with the components 1/(m + 3/2). Comparing the coefficients 
of B’, r = 0, 1, --- , on both sides of (3.6) we find 


— 


oe ay’, Le = 0, 3.8 
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(3.9) 


ar. = > Im/(m + 3/2), 


(3.10) 
the transmission coefficient 7’ becomes 
T = 8/2 1m T* (3.11) 
where Im denotes the imaginary part. We shall now show that L’” is a common left 
hand factor of all the matrices L at oy 
mati 


such that the right hand factor is a bounded 


el p 1, 2, 3, --- and q = 0, 1, 2, --- . and let S™™ = (s?”) and 
. pe } fhe matrice defined by 
s = @ if n>p+m| 
(3.12) 
S 0 ll n> @ 


(—1)’G(n, m, p)/H(n, 


(3.13) 
: 1 + 3/2), (3.14) 
? ot nm, m, f 


G(n, (—#+ 3/2), r(2t —n + mT(m + 3/2 
Hin, rié¢+ YrOr(it + m —n + IEE + m + 3/2)(3/2) 
iy 
a a ee Ce ah a (3.15 
P The element in the n-th row and m-th column of L“”’S'~” is 
Vr (-1)’ Cin + 3/2) C(m + 3/2 > (3.16) 
’ i <a 2.10) 
1 pi(p — 1)! n! I'(m + p + 3/2 
where, because of (2.1) and simple properties of the Gamma function 
= - (n +r!) T(r + 3/2) (—p + 3/2), (2p + m — 1) - 
> > l 7 - I 7 (3 17) 
hand had y! Tinac+r oe 5 ‘2) (3/2), I'(p _ m—r +- 1) ° 
? TT 2p + m)T(3 2) bs (n 4. 1). (3 2 iam Pp) i—9 _ m), (3 18) 
(p + m)!T(n + 5/2) = r! (n + 5/2), (1 — m — 2p). ties 
The sum in (3.1 


18) can be computed by using Saalschuetz’s formula (cf. Bailey [4] for 
simple proof) which can be written in the form 
(a),(b) (—h) = (e — a),(Cc —_ b), 
—ricd(1+at+b—c—k), ©.le —a — dD), 
sc#0,-1—-2,--» -k—-—1;1+4+ 


(3.19) 


a+b—c#1,2,--:,k) 
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Takinga = n-+1,) —p+3/2.c = 5/2+n,k.= n+ m, (3.19) gives for 7 ie 
] \ f Le 
(3.17) 
sw _ nT (2p + m)T(8/2) (38/2 (ntpt Das, 320) 
5 p+ m)!T(n + 5/2 N + 5/2) mer(D) mop — 
From (3.20) and (3.16) it follows that L er rye ea proof of ” Fa De 5 errs 


follows by the same method. 


The elements of the matrices S‘“**”’ are zero except for those in the first g rows. 
This is not true for the S°” but the following lemma shows that S°” is a polynomial 
in S‘“’ apart from right hand factors which are either the identity or of the type of 
the S°*** 
We have: 
Lemma 2. Let p, t = 1, 2, 3, --- and let R'"’ be the matrix for which the element in 
the first row and m-th column is 
(—1)""  ——s_—sP(m + 3/2) (2t + m+ 1)_ 3 91) 
t — 1/2)t"t — 1)! Tim + ¢ + 3/2)(t + m+ 1)! aa 


=_ > l F ; , 

ss = eee. 3.22) 

LZ 
2¢+2 / 9\#+1 i ‘ y¥(2)) 4 t—u ‘ 
iS = > (-2)**"[(—t + 1/2) yas /(-1 — Dy] {SVR ™, 3.23) 

where, for u= tf, R de notes S i ‘ if Jere ral, 

(¢+ np)! I(38/2)T(—t — p + 3/2) nconser 

feed a E — eer (3.24) 


pit! l(—p + 3/2)r(—t + 3/2) 


1s a matrix in u hich all ¢ le ments are zero € xcept those in the first Pp TOWS. 
The proof of lemma 2 follows again from Saalschuetz’s formula. We have now: 


Theorem : If the equations 


Lt =n (3.25) 
. 7 % . ge (f ° (1 
have a solution, then all the vectors —”’ are determined by —~’ and by the relations — ' = O 
and the recurrence re lations 
c <2 @ — S c a —- se" (3.26) 
In the particular case where 
n=7 = (2/3, 2/5, 2/7 -), (3.27) 
we have 
ge" = (8/e, 0, 6, 0,. -- +), (3.28) 
0 


and at most the first r + 1 components of &°"’ are different from.zero. &, --- , &"’ are the 
solutions of the original system (3.6), if we use the first r + 1 equations for determining the 
first r + 1 unknowns and thereby neglect all terms involving the higher powers of 8 from the 
r-th power onwards. &, --- , €'” also determine the exact values of the first r + 1 coefficients 


of the expansion of T* in powers of B. 
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The proof of theorem 1 follows immediately from lemma 1 and in particular from 
the fact that the S~”’, S***” involve many vanishing elements. The uniqueness of the 
é”’, and the existence of the z,, (at least for sufficiently small values of 8) will be proved 
in section 5. 


4. Limiting cases for the matrix L. Let 
P(t) = T(t + 3/2)/TC + 1), Q(t) = T(t + 5/2)/T(t + 1). (4.1) 


Then Theorem 1 states that the equations 


> In.m(B)tnm =h (n=O, 1, 2, ---) (4.2) 
can be solved by formal (i.e. not necessarily convergent) power series in 6 if the equations 
WE = \ “P(n) p LmP?(m)/Q(n + mb = {4h,} (4.3) 
have a solution x, = 2,,. We shall investigate (4.3) together with the limiting case 
B - . Levine and Schwinger [1] have shown that then (4.2) tends towards the system 
of linea equations 

Lee ( ) r,/(n + m+ 2} = oth! (n = 0, 1, 2, -+*) (4.4) 

Vi here M is a constant. 


We have to define first the linear space of admissible solutions z,, from the nature 
of the problem. Since (3.1) is supposed to define the field in the aperture, and since 
the field cannot have a singularity in the center of the aperture, we must assume that 

lim >> x,(1 — 6” (4.5) 
€ > m=0 
exists. Since the original system (3.6) was set up merely in order to define the trans- 
mission coefficient, we shall assume that 


> x,,/(m + 3/2) (4.6) 


m=0 


converges. This implies, that 
do tut” (4.7) 


converges for | z| < 1 and therefore that the z,, actually define the field in the aperture. 
The hn we prove fi 

Lemma 3. If the vector & with the components x,, satisfies (4.5) and (4.6), then the 
operators L’”’ and L'”’ are defined for & in the sense that the sums in (4.3), (4.4) converge 
forn = 0,1,2 


Proof: Let Q(t) be defined as in (4.1) and let 


rst 


m 


Tm = Q(m)/Q(n + m), Cc. = » z,/(r + 3/2). (4.8) 
r 0 
Then the partial sums of the series in (4.3) are 


m—1 


Zz 7,t,/(r + 3/2) = > (tr, — Trarlor + Tom (4.9) 


r=() 
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where 


aT, — 27, = 8nP(r t+ D/{Qna+nin +r + 5/2]}. (4.10) 
Since the | o, | are bounded and » ee T, — Tra, | converges, the sums in (4.3) also con- 
verge. The proof for the convergence of the sums in (4.4) is even simpler. 
Theorem 2. If the equations L’’& = {h,} or LE = [h*} havea solution é } or 
& = ja, } satisfying (4.5) and (4.6), then the integral equations 
1/2 } 1/2 r € fsb 
| 1 — v)'7(1 — v2)" dv = 407"? DY 2"h,n!/(3/2), , 1.11) 


have analytic solutions 


. m_(0) ry 9/9 . o) = ‘ 
fe) = Yo ocT(m + 3/2)/m!, f*®% = Drv". (4.13) 


m=(0 


The solutions are unique and they also solve the problems of moments 


| f@~ A — Qo" dv = 40" h,n!/(3/2), ; | f*(yo"** dv = h*. (4.14) 


The integrals in (4.11) (4.12) are defined by 
al ni~e 
} =lim] . (4.15) 


“0 € » #0 


Since a formal expansion of the left hand sides of (4.11) and (4.12) leads to the linear 
fh,} and L'’& = {h*}, it has only to be shown that, under the as- 


§ 


equations L“’¢ = 
sumptions made about the z,, , such an expansion is legitimate. It suffices to prove that 


al—e 


lim | fMA — p)'’*y" dv = 8x ‘h,n! (in + 3/2) (4.16) 


where now f(v) is defined by (4.13) and h, by L°’& = {h,}. Since it follows from the 
< 


assumption (4.5) about the z,, that f(v) converges absolutely and uniformly for 0 S 


ri 
) 


v <= 1 — «e, we may integrate term by term in (4.16). Putting Y,, = z,, T(m + 3/2)/m! 
this gives (with v = (1 — «)W) 
yy. f oma dh 

1.17) 


= > Y,(1 — 6’ n+m+1)'F(-1/2,n+m+1,n+m+2;1—.6 (4.18) 
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According to Gauss’s formula (cf. Whittaker-Watson [3]) 
m+ 231) = (n+ m+ IIT(3B/2)/T + m 4+ 5/2), (4.20) 


F l/2,.n+m+i:nc+ 2; 


lemma and from lemma 3 it follows that 
> Y,(1 — « 'T(3/2)(n + m)!/(T(n + m + 5/2) = 4n7'*h,n!/(3/2), . (4.21) 


Now we have to show that the second sum in (4.19) tends towards zero as e — 0. Because 


of (4.5) it suffices to show that 


( € 
r(m + 3/2)[m!]"[n + m + 1) {F(-1/2,n + m+ 1;n+m+2;1— 6 (4.22) 
— F(—1/2,n + m+ 1;n + m + 2; 1)} 
r(m + 3/2)(2m!)"' > [1 — (1 — &)**"] 
id (4.23) 
- (1/2),/ {kK + Din + m+k+ 2)} 50 
as € 0 uniformly in n, m. We can prove that |c,,,,(e)| < ¢€ by observing that 
| 1 — e)**' < (k + 1)e. This and (4.23) gives 
( é < «I'(m + 3/2)(2m!)’ > (1/2),(n + m+ k + 2) "fk" 


2m'(n + m + 2)}7'F(1/2,n + m+ 2;n + m + 331) 


9 


= «F(m+3 


(4.24) 
eP'(1/2)0(m + 3/2)(n + m + IN Qm!IT(n + m + 5/2)" 
ai /2 » <L ) Ss ee , 
2 im - at. a! (m+n + I) < en/2 <e. 
2 (m + 3/2)(m + 5/2) --- (m +n 4+ 3/2) 


The uniqueness of the solution follows from 
Lemma 4: If >°2_, x,,/(m + 3/2) converges, then forO0 S v < 1, (1 — v)*” | f(x) | is 
1. The proof follows from summation by parts with the notation (4.8) and from 


bow 


the remark that 


> r(m + 5/2) | o, ;v"/(m + 1)! S Cll — v)*? — 17, (4.25) 
where c does not depend on ?. 


Now we can show that (4.3) cannot have a null solution. Because then the difference 


of two solutions of (4.11) would satisfy 


ey ee (4.26) 


and therefore: 


. ee (4.27) 


| A(v)\(1 — v)'?701 — vv" dv = 0, n = 0, I, 2, 


would be a function continuous in 0 < v S 1 according to lemma 4 
and therefore (4.27) shows that ¢(v)(1 — v)°” would be identically zero. 
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Conclusions from theorem 1. The equivalence of the equations L“’t = {h,,! and 
Lt = {hx} to a problem of moments shows that these sets of linear equations are 
unstable in the following sense: Not only may these equations have no solution at all, 
but this is certain to happen if we start with a set {h,,} of right hand sides for which 
a solution exists and then change a finite number of the h,, by an amount however small. 
In this case there does not even exist a continuous function f(v) which satisfies (4.11) 
or (4.12) with the modified right hand sides. 

The integral operators in (4.11), (4.12) are extensions of the linear operators defined 
by L™’ or L‘*’, since (4.11) or (4.12) may have a continuous solution f(v) which is 
not analytic. Consequently, a quantity like the transmission coefficient 


al 
‘ big = | f(vy' , q= » & Cm/\ Mm + 3 2) ( 1.28) 
“0 m=O 


can be defined even in cases where the z,, do not exist. An easy example is offered by 


the equat ions 


n+ m+ 2) = u/(n + 3/2), n= 0,1, 2, ---) (4.29) 


' 


which were also investigated by Levine and Schwinger. The corresponding integral 
equation is 
+1 sl 
| fel — oW)' dv = > W"*/(n + 3/2) =u | v'?/1 —oW) ah (4.30) 
* n=0 “0 
which gives 
fo=w”, =z. (4.31) 


In this case no set of x,, satisfying (4.29) can exist. However, it is possible to find sequences 


of constants Y,,’ such that 
2. Y, (m+n+2) =y (4.32) 
exist and 
lim = fun’ — p/(n + 3/2)!° = 0, lim z Y,. /(m + 3/2) = gz. (4.33) 


For this purpose, we can choose the Y,, from 
> VL" = (1 — v)‘(1/2),/k! (4.34) 


The right hand side in (4.34) is a polynomial which approximates v ““, since it is the 


(r + 1)-th partial sum of [1 — (1 — v)]"’”*. Clearly, the Y,7’ ~© asr—-o. 

5. Uniqueness and existence of the solution. Once a vector “” has been determined 
such that L‘’’é‘” = 7, where 7 is the vector of the right hand sides in the original equa- 
tions Li = 7, we can determine é from 


Mt = (5.1) 


1953] INFINITE MATRICES 85 
where, for all values of 8, M is defined by 


M=sg 4. re 4 Ps gpatt gitar (5.2) 
p=l 


q=0 


Here 9 denotes the identity. We shall call a vector — bounded if ) 3: t,, |> <@ and we 
shall call a matrix M bounded if there exists a constant U > 0 such that for all bounded 
vectors &: 

t*M'*Mé < U? > [é,. |? (5.3) 
where MM’ is the transposed matrix of M and an asterisk denotes the conjugate complex 


quantity. U is called an upper bound for M. It is well known that, if U, is an upper 
bound for S‘’(r = 1, 2,3 ---), the matrix M in (5.2) has a bounded inverse M~* if 


> 


p BU, <1 (5.4) 


M~* can be obtained from a Neumann series. We can use this in order to prove: 
Theorem 3. Let L, M, 9°”, — be defined by (3.5), (5.1), (3.27), (3.28). Then MM 
exists and is bounded for sufficiently small values of | 8 | < By and the equations Lt = 9 
have exactly one solution & which satisfies (4.5) and (4.6), namely —& = M7". 
Proof: Let V‘”’ be matrices such that 


{a + 5 as” Mg + , o av} = J. (5.5) 


r=() 


It is easily seen that the V‘” can be obtained from the S‘”’ by recurrence formulas. 
Let U‘” be upper bounds for the S‘” and assume that there exist constants 2, such 
that 


(1 - . U,)(1 + > s"9,) = 1. (5.6) 
This is true if 
1- , i B'U, (5.7) 
r=2 
is convergent and positive for 0 < 6B < 8. Then it can be shown that Q, is an upper 


. rir) as . -1,.(0) 
bound for V‘"’. Since it can also be shown that z,, (the m-th component of § = M™'"") 
is equal to the m-th component of 


be ave, (5.8) 
it follows that 


sce, F (5.9) 


From this it can easily be shown that for | 8 | < 8 condition (4.5) for the 2,, is satisfied. 
This proves the existence of 1/~* and of a bounded £ satisfying (4.5), (condition (4.6) 
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is always satisfied for bounded &) if we can find U. which are sufficiently small. We have 


Lemma 4. The matrices 


te pei - ai (5.10) 
havea nner bound 
2 1/2 1/2;_2 1/2 9.2 w\1/2yt+274 fc / 
r(r — 8) t, 2°"(x — 8) t!, (Qr° — 16)°°°2'**(1/2),/(¢ + 1}, 
(5.11) 
2°°'(2r" — 16)'’?/(q + 1)! 


The proof is elementary but laborious and will be omitted since the upper bounds 


are not the best possible ones. 
In order to prove the uniqueness of the solution g 
(M — 9)é is bounded for every — merely satisfying (4.5); provided that 8 is so small 
from Lemma 4, converges. This can be proved by an elementary 


that (5.4), with the U 
investigation of the S’’’. Now if there is a &* satisfying (4.5) and (4.6) such that Lé* = 0, 
we would have J/¢? é* + ¢ where ¢ is bounded and L‘“’ é* + L“’¢ = 0. Now it follows 
from the equivalence of the operator L’’ to the operator of a moment problem (cf. 
0. Therefore £* is bounded, and since M™~' is bounded, £* 


= M~':’ we observe first that 


Theorem 2) that &* + ¢ = 


must be zero since J/é* = &* + ¢ = 0. 

No numerical values for the permissible ranges of 6 are given since it is entirely 
possible that the inverse M/~* exists for all values of 8. This seems to be indicated by 
a result of Sommerfeld and Perron [5] who showed that for the related problem of the 
freely vibrating disc the real part of a resulting set of linear equations can be solved 


explicitly and without restrictions, 
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A GENERALIZATION OF MODULATION SPECTRA* 


BY 
HAN CHANG anp V. C. RIDEOUT 


University of Wisconsin 


I. Introduction. A general theory of modulation spectra may be developed by the 
use of Fourier analysis. It may be applied to frequency as well as to amplitude modu- 
lation and is particularly valuable in the study of modulation products resulting when 
nonlinear devices such as rectifiers are used as modulators. In all cases it shows that the 
modulation products are harmonics of the highest common factor among the carrier and 
the modulating frequencies. Also, this approach yields some new results and some 
clarification of concepts. 

Of course the methods of Fourier can only be used where there is an integral re- 
lationship between the carrier and each modulating frequency so that the modulated 
wave may be treated as a periodic function. When this is not strictly true Bohr’smethod 
for almost periodic functions may be used. 

II. Outline of Theory. 

1. Modulation products. When two or more waves are combined in a nonlinear 
circuit such as a diode rectifier, a reactance-tube oscillator, or the human ear, new 
frequencies appear as a result of some characteristic (such as amplitude or frequency) 
of one wave being modified by another. Mathematically the process may be expressed as 


= Fle, » €2 “thy 5 Gels (1) 


The new waves, which may include waves of the same frequency as the original waves, 
are the modulation products. 

The principal ways in which modulation may be achieved for the simple case of two 
input waves are: 

(a) Mixing in a nonlinear circuit whose characteristic is representable by a finite 
number of terms of a power series. 


m=k 


Fi(e: , 2) = >> (ae, + be)”. (2) 


m=) 
b) Mixing in a nonlinear circuit consisting of a biased ideal rectifier whose forward 
characteristic is representable by a finite number of terms of a power series. 


mak 
F(e, , @2) = 1fae, + ber — E] >> (ae, + be, — E)”. (3) 
m=0 
Here 1| |] is the Heaviside unit function. The summation is zero unless the term 
In Square brackets is positive.) 
(c) Amplitude Modulation 
P(e, ,@2) = >> dX (an + bner)"(a, + b,e2)”. (4) 
*Received April 1, 1952. This paper includes material from a thesis submitted by Han Chang in 
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(d) Angle Modulation 


F4(e, ,€2) = K or {file; , ee) + foler , €o)}. (oO) 


COs 


Particular values of the constants in these expressions reduce them to ones which 
are more familiar in engineering practice. Thus if in (4) a, = 0,m = n = 1, one has ordi- 
nary amplitude modulation. 


F(e; , @2) = Anb,(1 + keseo , (6) 
where K = b,,a, is the modulation factor. 
In (5) if fi(e, , : k, sin e, é = k,wt, and f.(e, , €2) = ke e. then one has 
e l 1 1 1 e es 2 2 


ordinary phase modulation, or 


sin 
Fife, , @2) = K or (kywt + hve). (7) 


COs 


If the frequencies of the waves which are combined in a nonlinear circuit are com- 
mensurable (in the language of electrical engineering) or contain a common factor or 
factors, then the modulation products will have a common period which is given by the 
highest common factor among the original frequencies. In this case ordinary Fourier 
analysis may be used to find the spectral components. 

It is possible that the frequencies of the combined waves may be incommensurable 


or have no common factor. In such cases (1) may be treated as an ‘‘almost periodic 
} 


function”, the theory of which was first advanced by H. Bohr in 1925 [1]. Such a wave 
would never repeat itself exactly, but for any small quantity ¢ there is always an approxi- 
mate period 7 at the beginning and end of which the amplitudes of the wave differ by less 
than e. There are actually infinitely many such periods. 

The expansion of the almost periodic function is called a generalized Fourier series 
whose coefficients are found by a limiting process as follows: 


7 


. 


bo 


a, = lim = | f(t) cos X,¢ dt, (8) 
T— ] “0 5 
oa a 

b, = lim T | J() sma t dt. (9) 


It is to be noted that it is no longer necessary to ascertain , beforehand. If one 
replaces \,, in (8) by some variable z, the limit will be in general zero. The X,’s are then 
the values of x which could render these limits not identically zero. 

Therefore, to sum up, we see that the modulation products can always be analyzed 
into systematic spectral components by finding the Fourier or generalized Fourier 
coefficients. The result is often more revealing than the conventional trigonometrical 


expansion used in engineering. 


An important theorem on Fourier coefficients known as Parseval’s theorem [2] will 
be of use in the development of energy changes due to modulation. In its simplest form, 


this theorem states that if a function f(z) has its square summable in (— 7, 7) and if its 
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Fourier coefficients are a)/2, a, , dz, «+: 0, , b2 , bs, --+ then 
1 ® a. o 
ba 2 io 2 2 
tf U@ar= $+ D@+ wo. (10) 


This theorem fits our problem because it is well known that the average power or 
energy of an electric wave is proportional to the sum of the squares of the amplitudes 
of its Fourier components. 

Bohr’s “Fundamental Theorem” has proved that the generalized Fourier series for 
almost periodic functions still satisfies Parseval’s formula, providing again a theoretical 
basis for considering the energy changes in modulation. 

III. Application of the Method and New Results. The method will now be applied 
to the analysis of several typical examples of modulated waves whose conventional 
sideband expansions are known. The new findings and clarification of concepts will be 
evident where they occur. 

1. Modulation products with pure period. 

a. Simple Amplitude-Modulated Wave. This is usually given in the form 


e = A(1 + m cosw,t) cos w,t, (11) 
where A is the carrier amplitude, 
w, is the modulation frequency, 
w, is the carrier frequency, 
m is the degree of modulation. 
By a trigonometrical identity (11) can be written as 


e = A cosw,t + (mA/2) cos (w, + w.)t + (mA/2) cos ww. — w,)t. (12) 


It can be shown that these sidebands are actually harnonics of a certain fundamental 
in (11) as follows. Let the highest common factor between w, and w, be wo , such that 
W, = NW , We = Ng wW. Then T = 27/w. , and 


T/2 
a, =/ A(1 + mM COS N,wol) COS N.wot COS Nwol dt, 
0 


(13) 
b, = 0. 


The cosine coefficient in (13) gives 


> Tr 
a= A cos Not cos Niwot dwt 
FT Jo 


~ [ mA[cos (n, + N.)wot + cos (N, — Nq)wWol] COS Nwot dwot. (14) 


Thus a, is not zero only when n = n, and when n = n, + n, giving 
a,, = A, Q,,o,, = mA/z. (15) 


Therefore the three terms in (12) are the (n,.)th, the (n, + n,)th and the (n, — n,)th 
harmonics of a wave of period 2x/w, whose fundamental and other harmonics are zero. 
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b. Simple Frequency-Modulated Wave. Assume again the simplest form with con- 
ventional notatio1 


e = Asin (wt + m,;sin w,!). (16) 


Since this is an odd function, one can safely ignore all the cosine coefficients. Then if 


@ is again the highest common factor between w and "i 


( - | A sin (n.wyt + mM; SIN N,wot) SIN Nwol daot. (17) 


iw 
Mee . 
s m, sin 2 ao (m,) sin (2k + I)a ISb 
, = 1 } ee ] . 
then (17) can be exp ded and integrated to give 


b, = A[J,.(m,) — J,,(n | (19 


whi I S } has values w hich are positive or ihe gvative inte vers including 
, “a er ‘ ; : 

zero and n, has values which are positive integers and J,(m,) is the 

Bessel coefficient the first kind of order s and argument m,; . When s is negative, 


J .(m —1)"J_,(m,). 20 


There are sper lues for n, for which (19) will actually involve two terms as 
given but otherwise there will be but one term. This can be seen as follows: n, and n, are 
prime to each oth d therefore both are odd numbers or one is odd and the other even. 
Because ? t 2 where s, and s, are integers then for n, and n, both odd it is 
only ] ssible to |} 1. If n, iseven and nN, is odd it is only possible to have nw = Zz 
while if 7 is odd : i even It is again only possible to have n, i 

In view of th equation 19) can be written as follows: 

If ) # Vi het he . is even or odd, 


b. = (A/2)[J,...(m;) — Je.a.(m)(—1)*""'”. (22) 


If n l or 2, either even or odd, 


b, = (A/n,)|J,.-.(my,) |(— I) sale (23) 
To see what these coefficients really mean, take n, = 26,n, = 7. Ifn = 1, (i.e. consider 
the fundamental of the wave), then (n. + n)/n, = (26 + 1)/7 is not an integer indi- 


1 


cating that the frequency component at w, is zero. 


Next let nm = 2; then (n., — n)/n, = (26 — 2)/7 ¥ an integer, and (n. + n)/n, = 
(26 + 2)/7 = 4 and 1)* = 1 indicating that the second harmonic at 2w, of magnitude 


—AJ,(m,) exists 
A continuation of this process will show that the fifth harmonic exists and is of 
magnitude —AJ;(m-). The ninth harmonic exists and is of magnitude AJ,(m,) and the 


l- 


if 
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twelfth harmonic exists and is of magnitude AJ.(m,), and etc. This result is shown in the 





It is interesting to note that this is the same spectrum as if in the ordinary expansion 
those side-bands of negative frequencies were reflected at the zero frequency axis with 
sions reversed. This diagram also shows that in general there will be additional frequency 

a . . 1: “Sano ye 
components sand hed between the ordinary sideband spaces, for; example, those 
between the carrier J, and the first sidebands J, 
3 Js J, Jy. J, 
a 
ov Spe ete 
-" Ys vs I, Jf ff requency 
Fic. 1. FM sideban ls, Ne = 26, Na = 7, (magnitudes not to scale 

\nother defect of the ordinary expansion appears in the special case when n, = lor 

when 7 2 and n, is an odd number. Thus, if n, = 1, n, = 9, say, then application of 


19) will show that the fundamental is of magnitude 
A[Js(m,) — Jio(m,)], 
and the second harmonic is of magnitude 
— A[J-(m,) ae J 33h )]. 


And if n, 2,7 9, the application of (19) will show that even harmonics do not 
exist and the odd harmonics each involve two terms. These are shown in Fig. 2. 





Sidebands 
-_~+ 


J,~ Ji, b Jig J,° Jis J-Si2 frequency 





(a) 





3 J,- Js Jz- J; J.- Jo 

9 | 

2 = < - : 

3 | | 

DN Js~ J, U-Je frequency 





(b) 
Fic. 2. (a) FM sidebands, ne = 9, me = 1. (b) FM sidebands, ne = 9, ma = 2. 
It may be seen that the ordinary expansion is inadequate for these particular cases, 
for it gives the correct frequency components but not the correct amplitude for each 


component. 








92 HAN CHANG AND V. C. RIDEOUT [Vol. XI, No. 1 


c. Multi-tone Modulation. The same procedure applied to multi-tone modulation 
shows that for an AM wave of the form 


k \ 
e= A(1 + 2. mM, cos wt} COS wt, (24) 


there exist the Fourier coefficients 


(1 + 2 m, COS 2.) COS wt COS Nwyt d(wol). (25) 


rT J0 s= 1 


| 
Li) 
3 — 
ee, 


where w, again is the highest common factor among a , w2 --* o , and each sideband 
can be identified as one of these Fourier harmonics. 


For the FM case 
e = Asin (wet + bs m, COS w ‘). (26) 
s=1 


The Fourier coefficients, 


9 { 7 ‘a , 
a, = - | sin (nt + a m, COS n wnt) COS Nwot Awol), (27) 
Tr Jo s=1 
2a if - ‘ 
b. = | si |N.wot + > Mm, COS N,wWof} SIN Nwol d(wot), (28) 
T Jo s=1 


can be obtained and one has, 


a, = | = | I] Ju.lm,)(—1)""") + ry (II J,..(m,)(— 1)" ‘Ia (29) 


22 ke =| 
where 
RK K 
> kan, | = |ne +n |, > kaon, | = |n. — ni], (29a) 
s= s=1 
and 


i 
ky, = >» Ka or bis ae (29b) 
s= s=] 


when they are odd numbers. 


> IT J,,,(m,)(— 1)"" ) _ > ( 


\ \ 
kas 8 


5. = J,.,,(m,)(— 1)"* ‘) Ja (30) 


where 
K K 
> k.n,| =|n—ni, pa k,on,| =|n, +n |, (30a) 
| s=l1 
and 
K K 
k=D> ki or Dike, (30b) 
s=1 s=1 


when they are even numbers. 
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Therefore, we have 


e = > a, cos nat + >> b, sin nat, 
n n 


which is considerably different from the ordinary result, 


an a » 
Pica 2 {Tl J.dm)} cos (> kat). (31) 


ke =0 s=1 
d. Modulation Products from a Linear Rectifier. The subject of heterodyne detection 
has been investigated by many. Engineering practice assumes that detector can follow 
the envelope ideally so that higher harmonics can be neglected and the difference fre- 
quency taken as the fundamental. W. R. Bennett [3] gives a double Fourier series de- 
velopment of the output of such a rectifier for any amplitude and frequency ratio which 
seems to be the only exact analysis which has so far appeared. The present method of 
analysis appears to be an interesting and useful alternative. 
Express Heaviside’s Unit Function by, 


Wy =54+2f Pde. (32) 
“ T Jo Ww 
Assume an input wave of the form 
e = A cosa,t + B cos wol. (33) 
From (3) and (32) the output wave from a zero-bias linear rectifier is 


J l re? si an 
@, = 15 + 8 | nwt ay 


“a0 w 


= : (A cos wf + B eos Wot) 


~ 





4 A cos wt + B Cos wot | sin w(A cos w,t + B cos wot) he. (34) 


T /0 WwW 
Let wot = x, where w, is again the highest common factor between w, and w, , such 
that w, = %; @» , @ = Ng». Then the Fourier series of (34), of cosine terms only, can be 


specified completely by 


2 oe 
a, = | E (A cos n,x + B cos np2) 
T Jo “a 
A cosn,x + Beosn x [* sin (A cosn,x + B cos nox) - 
4 an <n dw | cos nx dz. (35) 
us “0 ® 


The first part of (35) gives 
a,, = A/2, a,, = B/2. (36) 


The second part can be written as 


I, = cd [ [ (4 cos n,x + B cos na) ( Foxe. ox+1(Aw)(—1)* cos (2K + Dn) 
K=0 


“0 “0 


: (EoxJ ox(Bw)(—1)* cos 2Kn.x2) + (> EoxJ2x(Aw)(—1)* cos 2Kn.z) 
K=-0 


; (> Eoxs1J ox+1(Bw)(—1)* cos (2K + mz) cos nx dx “ 
K=0 
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39) 

D 2, ~ 0) 
re in (37) infinitely many terms in the integrand, each term involving fou 
cosine functio s multiplied together associated with a product of two Bessel coefficients 
which are const t if we integrate with 1 spect to x first. Integration with re spect to x 
shows that each term in the integrand is not zero o1 1} when the sum and difference of 
two of the foul cosine oles are equal re specti ely to the sum or difference ol the other 


two al le Ss. Follo ving this the formula below is obt: ined. 


If 7 S eve nd ? 1s odd, thei for n odd 


B f = : . dw 
+=—}| Dd Eoxei(—! *Jox+:(Aw)J2s,(Bo) — (40a) 


where 


tn + 2K 
28, +1 — (40b) 
4i 
28, = tn+tn. + (2K + I)n, (400) 
n 
and for n even 
{ j : 4 de 
a= 7 | th | —]) J 1c J Be 


ay ’ 4 dw : 
to—] DO Ensies(—1)****J2x(Aw)J25,+1(Bo) —, (41a) 
where, 
28, = | = = a (41b) 
+n + n, + 2Kn, 


(41¢) 


No 


If n, and m2 are both odd, all odd harmonics are missing in the output because then 
in (42) there will be no possibility of combining the angles such that the whole integral 


0. 
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is not zero. For the even harmonics a, is given by (41a) with (41b) and (41c) replaced by, 


+n + nm, + (2K + 1)m, | 


28S, = 
;' “me (42a) 
and 
9 
8.41 +n + MN. + 2Kn, | (42b) 


No 


In these formulas, the summations should run over all possible integral values of 
S and K that may satisfy the Diophantine equations in absolute value form. 

Without loss of generality, we can assume that A > B (The case A = B will be dis- 
cussed later). Then (40) and (41) can be integrated as a special case of the infinite dis- 
continuous integral of Weber and Schafheitlin [4]. Thus if n, is even and n, odd, then for 


n odd, 


E,m sx-0 \A/ 2K — 28, +1 
I 9 r(28, -+- 2 
f2K +2S8,+1 —2K +258, +1 2) 
F4 ,) 9 28; + 2 ay 
\ - “ ) 


4 2k — 28, + 3 
; r( a )res, +1 
2K + 28,41 —2K + 28, —1.. B’\ 
‘ re ; 9 ] : +. ; ; 2S. + Ss >, (43) 
\ 2 é A ) 
Here equations (40b) and (40c) have to be satisfied by 2S, + land 28, . 
For nmeven 
9K + 298 — 
’ a (24 oz 1 + Ny 1 
- A i (2 \ 2 2 ' d ‘ 
me Eur s“t-0 \A 9K — 28 2 
(24 St = \res, +1) 
J2K + 2S, +1 —2K + 2S, —1 | 
r(24 + 28, + ) 1)K+8s 
B a. 2 ail 
+ re “8s A 2K — 28S,+ 1 ss 
J2K + 28, 1 —2K 28, | - 3° 
* ri? * r 9 : + a aie Se ok : + 28, + x ” (44) 
2 2 A‘) 


where 2S, and 2S, + 1 satisfy (41b) and (41c). 
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If n, and n, are both odd, only even harmonics exist in the output and these are 
given by (44) with the Diophantine equations replaced by those of (42a) and (42b). 
In the foregoing equations, I'(#) is the gamma function of argument z and F(a, b, c, x) 
is the hypergeometric function of parameters a, b, c and argument z. 
If A = B, (40) and (41) still exist, but the hypergeometric functions simplify to gamma 
functions, i.e., 
re)T(c — a — b) (45) 


F(a, 6, ¢, 1) = I(c — a)I(e — b)° 


This Fourier series analysis of the output shows two interesting results. First, it 
shows that the output may have important components of frequency lower than the 
difference frequency. When the difference frequency is not too much smaller than the 
beating frequencies and if it is not the highest common factor, there will be beat tones of 
considerable amplitude at frequencies lower than the difference frequency. Secondly, 
for a particular frequency component in the output, these formulas give the amplitudes 
of all components provided the Diophantine equations are solved for all the possible S 
and K. For illustration, suppose f, = 800, f. = 1400 such that the highest common factor 
is 200, and n, = 4, n. = 7. Then instead of the difference frequency 600, one would have 
the series of frequencies of values 200, 400, 600, 800, 1,000 etc. in the output. The fre- 
quencies 200 and 400 are both lower than the difference 600. On the other hand, suppose 
fi; = 210, f. = 330, so that the highest common factor is 30 and n, = 7, n2 = 11. We 
would expect the series of frequencies 30, 60, 90, 120, 150 etc. in the output. Here, how- 
ever, since n, and n, are both odd, the odd harmonics 30, 90, 150 ete. will be missing 
except 210 and 330 which will be present, (Eq. 36). 

e. Arbitrary Wave-Shape Modulation. It is evident that the same analysis for fre- 
quency modulation by any arbitrary wave-shape would exhibit the same reflected side- 
band phenomena as the case of sinusoidal modulation. A particular case of rectangular 
wave modulation has been analyzed in detail and the result compared with the ordinary 
spectrum to verify this conclusion. The result, however, does not seem to deserve more 
space here. 

2. Modulation Products With No Exact Period. When a common factor among the 
component frequencies present in the modulation products does not exist, the wave 
as a whole is not periodic. The theory of almost periodic function then relieves us of any 
possible logical confusion. Here a definite spectrum still exists; the frequency compo- 
nents, however, are no longer related by a multiple of a common component but are 
determined by certain characteristic values as explained in connection with (8) and (9). 

To illustrate the principle, consider again the amplitude-modulated wave e = A 
(1 + m Cos wt) Cos w,t. If w, and w, have no factor in common, then the conventional 
Fourier analysis no longer applies. However by, (8), 

aT 


An lim A | (1 + m cos w,t) cos w.t cos At dt 
T-@ T “0 


eT 


lim 4 | [cos (w. + A)t + cos Ww, — A)é] dt 
T-@ r “0 


I] 





T 
+ lim m4 [cos @, + w. + At + cos (w + ow, — A) 
T-@ “0 


+ cos ww. — w, + A)t + cos, — w, — At] dt. (46) 





na 


5) 


ie 


of 


ifs ¥ 


Fd 
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After integration each term in (46) will be of the form 


lim A sin (K + Ajt ig (47) 
fee © K+2X |’ 
and will be identically zero except when \ = + K. Therefore, since we are dealing with 
real or positive frequency only, the characteristic values are \, = w, , A2 = w. + w., and 
As = w, — w, . It can easily be shown that 


a, = A, a, = mA/2, a; = mA/2. (48) 


These are the same as the sideband amplitudes obtained in the periodic case. 

The other kinds of modulation, which so far have been considered only for the case 
where the carrier and modulating frequencies have a common factor, may similarly be 
handled by this method for the almost periodic case. Also the results may be inferred for 
the case where w, and w,, have a very small highest common factor and almost have a 
larger highest common factor. Here we would expect that the true modulation products 
would be large for the cases where they most nearly coincided with the products for the 
period based on the large “highest common factor’. 

IV. Findings on Energy Relationships. It is well known that in amplitude modula- 
tion the modulated wave has its energy increased by an amount corresponding to that 
in the sidebands. From the point of view of Fourier’s series, this finding is nothing but 
an application of Parseval’s theorem, since the energy or average power per cycle of an 
oscillation is proportional to the average square of the wave, and hence equal to the 
sum of the square of its Fourier coefficients. Thus, for the wave in equation (11) the 
energy E is, 

E= >) @= A’? + . A’m + J A’m = A? + 5 A’m? (49) 
res 4 4 2 
which agrees with the usual result. 

For an FM wave, the coefficients obtained from (21), (i.e. when n, = 1, n, either 

even or odd) give 


E= B= A? DY [J2._.(m) + J2.4n(m,) — 2Wa,-n(m)Snesn(m,)], (50) 


n=1 


and the coefficients from (22), i.e. for n, = 2, n, odd, give 


E= B= A? DY [Vi km) + Treraxes(my) — Wa-—x(m) In. eei(m)] (51) 
—- K=0 
where n’ = (n. — 1)/2, and 2K +1 =n. 

Equations (50) and (51) can be simplified (Appendix) to 


E = A*[1 — J2,.(2m,)], (52) 


and a 


Il 


A?|1 — J,,(2m,)] respectively. (53) 


The general case when n, ¥ 1 or 2 as in equation (23) can be easily analyzed by 
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he fact that as n runs from 1 to ~, the Fourier coefficients J;(m,) will each occur 
= 0. Therefore Parseval’s formula gives, with 


using t 
0 and only once when 


twice when 
the aid of the relation 


co 


J> (my, + 2>0 Jz (m,) = ¥ 
1 


: (54) 


E = aa m;) +2 a Ji(m | = A’. 


(52) and (53) show that in those two 


Since the energy before modulation is A’, 
particular cases frequency modulation decreases the wave energy because Jx(x) 1s 
54), the energy remains unchanged. 


always positive if x < K. In the case of 


Appendix 


53) may be derived by starting with equation (50) in the text, 


Equations (52) and ( 


E = A’ > [J _ (my) + Je.an(ms) — RI n.-n(mg) In. +n(my)]. (1) 
Let n. + n = m, then n, — n = 2n, — m. 
When n = 1. m = n. + 1, and when n = ~, m = o~; therefore, 
E = A’ Yo [.J2,..—m( my) + Salm) — Won.—m(Mys)J m(m;)]. (3") 
Now 
> ; = > J one—-m( My) + = . 
= 3 J5(m,) + i J (md), (3") 
” m l 
and 
. b J2(m;) = pe J2(m,) — 2 J (my), (4’) 
so that 
> m,) + J2(m,)] = Fe I2(m) — Jim). (5’) 
Since J_,,(m, —1)” J,,(m,), pe P ids can be written as 
(6’) 


® 


~ 5 





cur 


ith 


By the Addition T 


If 


MV hich i 
By 


hich is equation 
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Equation (5’) becomes, therefore, 


be [Jone-m(™m,) + J,(m,)] = 1 — J5.(m,). (73 
The product sum 


+s 2 on. —m( My) J (My) 


m 


can be treated as follows: 


>  2Wp.-m( Mz) J m( My) 


= > [Jon.-m(m)Jn(m)] + >, 


m=ner+l1 


J one—n( My) oF m( My) 


- | } J one—m( My) mM) + . ip 


+[> J on.—m( My) J »( My -> J on,-m( My) J (m, | 





> Jn(my;) 


m(My)d one-m( My) +  : (my) J one—m( My )| 


m= - 


+ " J ong—m( My) J 


m=0 


m(mM;) — >» J on sm) Jf) 


Zz: Jon.—m(My)eJ »(Mm,) — J2.(m,). 


(8’) 
Combining (7’) and (8’) 


one has (1’) in the simple form 


E= ai _ > : oe sm) Samy |, 


Theorem of Neumann and Schlafli 


(9’) 


}) 


JAY +2) = Zz. oe) © ee 


(10’) 
) = z, therefore, J,(2y) = )o2_. J..(y)Jn—m(y) and (9’) can be replaced by 

E = A’*(lI J >,,(2m;)] (11’) 
is eq uation 52). 


an entirely similar process, it can be shown that (51) can be simplified to 


E = A’[{l — J,.(2m,)] 


(12’) 
53). 
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id 1. Introduction. A steady compressible non-viscous flow is conical if it contains a 
vertex P, such that on every half line through P the velocity components, pressure, 
density, and entropy are constant. Various linearized conical flows have been discussed 
by numerous authors. However, only three examples of non-linearized conical potential 
flow fields are known to exist mathematically: Prandtl-Meyer flow around an edge [4]’ 
with or without sweep; Taylor-Maccoll flow about a non-yawing circular cone [5]; and 
an axisymmetric flow through a convergent nozzle discussed by Busemann [1]. In this 
paper the construction of two new examples will be considered. Both contain two regions 
of swept Prandtl-Meyer flow. In the first, the boundary has been chosen to prevent 
them from interacting, and the hodograph is one-dimensional. From it can be obtained 
a flow, with attached plane shock, over an object resembling an airplane with a swept- 
forward wing of positive dihedral and with a thick fin. In the second, the boundary has 
been chosen to permit interaction and the hodograph is two-dimensional. It was studied 
originally in the hodograph space by one of the authors [2]. In the present treatment, the 
need to consider possible difficulties in constructing the flow field from a knowledge of 
the hodograph has been avoided by confining the discussion to the physical space. It 
should also be remarked that in both examples the second order partial differential 
equation for conical potential flow is of mixed type. 

A pair of these examples could conceivably be used to study a particular, atypical 
case of wing-body interference. Numerical results can easily be calculated, if necessary, 
with the aid of a characteristics table and by means of standard techniques for numerical 
integration of ordinary differential equations and numerical solution of characteristic 
initial value problems for second order hyperbolic partial differential equations in two 
independent variables. 

2. Fundamental Ideas. The velocity potential function of a steady irrotational 
non-viscous isentropic flow satisfies 


(a’6,; — uju;)d°¢/dx,; dx; = 0 (2.1) 
where x; (¢ = 1, 2, 3) are rectangular coordinates, 
u; = d/dz; (2.2) 


is the velocity component parallel to the x,-axis in units of the maximum speed of flow, 
a’ = 1/2Ay — 1)(1 — uu,) = 1/2 — DA —- ¢) (2.3) 


is the square of the velocity of sound, Kronecker’s delta 6;, = 1(0) if i = (#) j, and the 
convention has been adopted that repeated indices imply summation over their range. 


If g = u,u; > a, the flow is supersonic, and there exist real characteristic surfaces, 
which are envelopes of the Mach cones 
[u(x — 2,)]° = (q° — a*)\(xt — 2,)(x% — 2) (2.4) 


*Received April 17, 1952. 
1Numbers in brackets designate references listed at the end of the paper. 
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where 2* are 1 oordinates. At these surfaces the partial derivatives of ¢ of the 


second or higher order may have discontinuities while ¢ and 0¢/0x; remain continuous, 
as may occur when two solutions are patched together. 

The image of a flow under the mapping 2x; — u, is defined to be its hodograph. In 
sional flows have three-dimensional hodographs, but not the ex- 


le! 
amples in the following sections. Flows with one (two)-dimensional hodographs have 


le (double) waves. The relevant properties of simple waves, discussed 


veneral. three-din 


been called sim 
in [3], will be summarized briefly. For some function p(x), u; = u,;(u). Now (2.1) implies 


aul’, = (u,u’) = (qq’) (2.5) 
where u{ = du;/du. The hodograph is any curve obtainable from a Prandtl-Meyer epi- 
cycloid by deforming its plane into a cone with vertex at u; = 0. To specify the curve 
completely an additional equation and initial conditions are required. For example, in 
a swept Prandtl-Meyer flow the velocity component parallel to a fixed unit vector \, 


must be constant, that is, 


u;,\; = constant. (2.6) 
In the physical space the prototype of the hodograph point u;(u), a surface 1 = constant 


which also bears constant pressure, density, and entropy, is a plane normal to u/ . The 
exact locations of the prototype planes in the physical space are as yet undetermined. 
If all are forced to pass through a common point, the simple wave will be a conical flow. 
Finally, note that a prototype plane is a characteristic surface, all of whose Mach cones 
are congruent and have parallel axes. 

Now let 


A. = 2,/adea = 1, 2), 2= 2, w= Us (2.7) 


Later it will also be convenient to use X = X,, Y = X.,u = u;,andv = u.. Note 
that to a point (curve) in the X,X,-plane there corresponds in the z,-space a line through 
(cone with vertex the origin. Consider a general conical flow with vertex at the 
origin. Since u deo dx, are homogeneous of degree zero in the 2, s, g may be assumed 


with no loss of generality to be homogeneous of degree one. Then for some ® 


gy = 20(X, , X, (2.8) 
00/02, = Ob/dX, , w= de/dz = $— XU, (2.9 
0°9/OX, O Ob/0X, 0X, z 0°y/dx, 02 = —X,0°6/0X, 0X5 

200/02 = X,A 66/0X. OX (2.10) 

and (2.1) becomes 
a > oe 4 — (u, — wX,)(Ug — wXz,)] 0°B/OX, OX, = O. (2.11 
Characteristics for (2.11) are curves X, X ,(t) in the X,X,.-plane on which the partial 
derivatives of ® of seco d or higher order may have discontinuities, while @ and d® OX. 
remain continuous. To these curves there correspond cones Z,/2 = Ae) on which, by 
(2.10) or its analogs the partial derivatives of 0a) of second or higher order will have 


correspondi ig liscol tl! uities, while by S 5) and 2.9 y and 00 OX; remain continuous. 
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Thus «,/z = X,(¢) are conical characteristic surfaces for (2.1). A plane tangent to one 


of them along the ray z,/z = X, must be tangent to the Mach cone (2.4) based on 
u, = u;(X, , X,) with vertex at the origin and must correspond to the tangent to the 
characteristic XY, = X,(t). Hence the characteristic directions for (2.11) must be those 


of the lines through (X, , X,) tangent to the conic 


(u,X* + w)? = (q° — a’)(X*X* 4+ 1) (2.12) 


a- 


with running coordinates X*. The type of (2.11) is hyperbolic, parabolic, or elliptic 
accordingly as (X, , X,) is outside, on, or inside (2.12). For a conical simple wave (2.11) 
is of hyperbolic type (except possibly on a curve on which the type is parabolic), with 
straight characteristics. Furthermore, a region of hyperbolic type adjacent to a region 
of uniform flow must be a conical simple wave. 

3. Flow along a conical wall. It is well known that swept Prandtl-Meyer flow over a 
dihedral angle can be generalized into simple wave flow past a curved cylindrical wall. 
Now a generalization to simple wave flow in a trough with boundaries composed of plane 
and conical segments will be discussed. 

Consider uniform supersonic flow along a plane wall on which / is a line inclined with 
respect to the direction of flow by more than the Mach angle. Let P be any point on 1. 
Extend the boundary beyond / as a cone through / with vertex P, and eventually join 
the cone to another plane segment. Attempt to fit a conical simple wave to this boundary. 
Let the origin of coordinates be at P. Describe the boundary by 


x; = ru,(y) (3.1) 
where r and uw are independent. With no loss of generality assume 
vv, = l, vvi = 1. (3.2) 
Then 
vel =v’ = 0, vw!’ = —1. (3.3) 


uv; = 0 (3 .4) 


and the boundary condition on the cone becomes u; = Av, + Bv! for some scalar func- 
tions A(u) and B(y). By (3.3) and (3.4) B = A’, so 


u(u) = A(p)vis(u) + A’(u)v’(u). (3.5) 
Now (2.5) becomes a?A’”(v/’vt’ — 1) = (A” — a’) (A + A”)’, where 

2 l \ 2 re ) > 

a = =(y — 11 — A — A”). (3.6) 


For any curve on the unit sphere the curvature v{’ v/’ = 1. Also, as in Prandtl-Meyer 


flow with leading edge v,;(u), the normal component of velocity A’ 2 a 2 0. Assume 


that the simple wave is an expansion. Then A’ (A + A’’) > 0, and 


aA'(ui’i’ — 1)'? = (A + A”)(A” — a)”. (3.7) 
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The initial values of A and A’ depend, of course, on the original velocity and orientation 
of 1. 
The prototype planes through the origin satisfy 


ui(uz; = 0. (3.8) 


Join the velocity field (3.5), (3.8) to uniform flows at both ends. Let s be any streamline 
not on (3.1) which does not intersect the envelope of the prototype planes, defined by 


= 7n,(u), nn, =nu;’ = 0, nm; = 1. (3.9) 


Such streamlines exist if (3.1) and (3.9) do not intersect. As the second wall of the trough, 
also composed of plane and conical segments, choose the cone through s with vertex P. 
A special example of this type of flow can easily be derived from a Prandtl-Meyer flow. 
Retain as one side of the trough the plane walls of the original boundary, and use a 
conical stream sheet for the other wall. 

4. A conical simple wave. Insert a half plane into a uniform supersonic flow at a 
moderate angle of attack, and make the angle between the leading edge and the un- 
disturbed velocity greater than the Mach angle. On one side there will be a swept Prandtl- 
Meyer flow around the leading edge. On the other there will be an attached plane shock 
wave behind which there will be uniform flow (supersonic if the angle of attack is not 
too large) parallel to the half plane. Throughout the entire flow the component of velocity 
parallel to the leading edge is constant. Introduce a coordinate system with origin on the 
leading edge, z-axis parallel to the undisturbed velocity, yz-plane parallel to the uniform 
flow behind the shock, and hence normal to the shock. Discard that part of the flow on 
the side of the yz-plane that contains the downstream half of the leading edge. Reflect 
the remainder with respect to the yz-plane. So far the boundary, shown in Fig. 1 together 
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with the shock, consists of a dihedral angle with congruent sectors removed from each 
face. Since the two halves of the shock are coplanar, the uniform flows on the compression 
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side join continuously. On the expansion side u, is double valued in the yz-plane. This 
difficulty can be avoided by modifying the upper side of the boundary as follows. 
First, examine the hodograph of the swept Prandtl-Meyer flow. In Fig. 2 let ON* 
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be the undisturbed velocity in the yz-plane, and let N be the projection of N* onto the 
plane II that passes through the origin and is normal to the leading edge. During the 
expansion, the component of velocity normal to II remains equal to NN* while the 
component parallel to II traces an arc NE of a Prandtl-Meyer epicycloid, shrunk by a 
factor (1 — NN**)'”. Fig 1 shows that the z-axis is under the plane boundary, so the 
angle AOE is acute. Hence the epicycloidal are NE cannot have another intersection 
with the line NS, parallel to OA. Now construct the hodograph by subjecting every 
point of NE to the displacement NN.* Clearly, during the expansion from ON* to OE* 
the angle between the velocity vector OU* and the yz-plane (II) steadily increases 
(decreases). 

Now let A be the intersection of x = 0 and the plane Mach surfaces that pass through 
the leading edges and are based on the undisturbed velocity. Let P be any point on A. 
From the nature of the hodograph it is clear that the streamline through P for the simple 
wave in x = O turns immediately into x > 0 and stays there. Hence an entire conical 
stream sheet through A bends into x > 0. Thus it is possible to separate the two regions 
of swept Prandtl-Meyer flow by means of a symmetrical conical fin, the thickness of 
which increases with increasing sweep. Note that at the junctions of the fin and the 
original boundary wu; is parallel to x; . Accordingly, near the corresponding points X, 
(2.11) is of elliptic type. On the other hand, for very large values of XY, near the boundary 

2.11) is of hyperbolic type. 
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To obtain from this boundary a finite obstacle resembling an airplane, symmetrically 
terminate the wing, as in Fig. 3, at a trailing edge which is supersonic with respect to the 
uniform flows adjacent to both sides of the wing. On the lower side there will be an ex- 
pansion around the trailing edges. The flow field just described will be unaffected up to 
the first Mach surface in the expansion fan emanating from the trailing edges. In par- 
ticular, the leading edge shock will cease to be plane at its intersection with this surface. 
On the upper side there will be a shock attached to the trailing edges, ahead of which 
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Flow Thick Fin 





Figure 3. Schematic Representation of Finite 
Obstocle in Simple Wave Flow 


the original flow will be unaltered. If the trailing edges are straight, as in Fig. 3, the 
shock and immediately following flow for either half of the wing will be conical with 
respect to the corresponding wing tip. Possible trailing edges for the fin would be its 
intersections with the trailing edge shocks. If a thick wing is desired, the upper surfaces 
need not be parallel to the lower surfaces. Finally the upper surfaces need not even be 
plane, but may be cylindrical or even conical, with vertices at the wing tips. However, 
after such changes the simple wave flows cease to be conical with respect to the origin. 

5. Interaction of simple waves [2]. Return to the stage of the discussion at the end 
of the first paragraph of Section 3. Examination of Fig. 1 shows that cross sections by 
the planes z = + 1 would have the appearances of Figs. 4 and 5. Ez and E, are traces 
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of the leading edges. The traces of the boundary and shock wave have been shown 
mainly to clarify the preceding example. They really have no immediate bearing on the 


interaction example, since it will be necessary to consider the possibility of having to 
reduce the amount of expansion and the extents of the interacting simple waves. The 
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Figure 5 Trace of Tentative Flow by Z=t] 


boundary will actually be constructed later. The essential point is that the discussion 
starts from simple waves in which all of the straight characteristics pass through the 
centers Ep or E,. 

In Fig. 4 the simple waves have not begun to interact so the boundary need not be 
altered for z < 0. In Fig. 5 they begin to influence each other along the characteristics 
AB, and AB, , shown in a magnified view in Fig. 6. Since the type of (2.11) is hyperbolic 
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at A, the characteristics actually extend some distance beyond A as shown. In the 
simple waves uv, and w are known on AB z and AB, . Hence ® and 0@/dX, are known 
by (2.9). The characteristic initial value problem for (2.11) with these initial data will 
have a unique solution in some sufficiently small characteristic quadrilateral AB,CB, , 
where C may be assumed to be on the Y-axis. Since the initial data are symmetrical 
with respect to the Y-axis, so is 6. Extend ® beyond B,C as the velocity potential of a 
conical simple wave. In general, the corresponding straight characteristics will not be 
centered. However, if the double wave AB,CB, is kept small enough, their envelope 
and intersections can be kept arbitrarily close to EF, , well beyond the boundary to be 
constructed later. Between the two simple waves in X > 0 there falls a region of uniform 
flow with the velocity uj.» of Bz . Extend the definition of ® symmetrically with respect 
to the Y-axis. Finally, between the two non-centered simple waves there falls a region 
of uniform flow with the velocity u;, of C. 

14 ‘mains to choose a boundary in 2 > @. First let U > ‘ Une ; User ; and l denote 
the points where rays from (0, 0, 0) parallel to ujo , Uior , Uiex » and U3 intersect z = 1. 
Let Gp, starting from the right edge of the right hand non-centered simple wave, be 
on the line Epl’s,. Then ujoz is tangent to the corresponding plane through the origin. 
Let FG, be an integral curve of the equation dX/(u — wX) = dY/(v — wY) of conical 
stream sheets. The final part of the boundary in X > 0 consists of the segment Fl’; 
which corresponds to a plane to which u, is tangent. Extend the boundary symmetrically 
into X < 0. 

Note that near U’, (2.11) is of elliptic type. 

It is interesting to observe that an alternation of double waves, simple waves, and 
regions of uniform flow similar to that in Fig. 6 also appears in an intersection of simple 
waves which occurs when uniform plane supersonic flow expands into an infinite sector. 


This is shown schematically in Fig. 7. 
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—NOTES— 


INDUCED MASS WITH VARIABLE DENSITY* 
By GARRETT BIRKHOFF (Harvard University) 

The concept of induced mass, and some of its properties, are extended to the case 
of an incompressible fluid of variable density. The proofs parallel closely those given 
recently by the author, for the case of free boundaries’. 

1. Minimum principle. Let the fluid, supposed originally at rest and occupying a 
region R, be given an acceleration a, by the motion of a wall W bounding the fluid in- 
ternally. By continuity, 

a, = f(x) on W, (1) 


where {(x) is the normal acceleration of W. 


THEOREM 1. The acceleration kinetic energy 


T=] pla-a) dk (2) 


“a/R 
is minimized, relative to all other volume conserving flows satisfying (1). 


Proof. Let a + b be any other volume-conserving initial acceleration satisfying 
(1). Then Div b = 0, and, by (1), b, = 0 on W’.. Consider next the expansion 


T=7| 


JF 


a 


p(a + b)-(a +b) dR 


+ 


| p(b-b) dR. 
JR 


Ni- 


= 7+ | (pa) -b aR + 
“R 

Since the last term is positive unless b vanishes identically, it is sufficient to show that 

the middle integral is zero. But, since Div b = 0 and (by the equations of motion, neglect- 

ing gravity’) pa = —Vp, 


Div (pb) = p Div b + (Vp)-b = —pa-b, 


where p is the scalar pressure. Hence, by the Divergence Theorem, letting b, denote the 
outward normal component of b, 


° 


| (pa)-b dR = [ pb, dS = 0, 
JR Jw 
since b, = 0 on W. (The same conclusion will hold if the boundary is partly “free’’, 
since then we can take p = 0.) This completes the proof. 

The case of a free boundary (say, of an internal cavity) should not, of course, be 
confused with the boundary of an incompressible region of zero density. 


*Received April 17, 1952. 

\* Induced mass with free boundaries,” this Quarterly 10, 81-86 (1952). 

*Gravity can be neglected, for very rapid (“impulsive’”’) accelerations; viscosity is without effect in 
the case of initial acceleration from rest. 
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2. Induced mass tensor and momentum. Letting a’, a’, a® and a‘, a’, a° be the 
fields associated with wnt infinitesimal translations and rotations, one 


acceleration 
by 


can define, as previously , the symmetric 6 X 6 induced mass tensor (matrix) || 7, 


. 


Tas = | pa-a dR = 4 ‘ (3) 


“R 
Since the diagonal components 7’,, satisfy Theorem 1, we obtain immediately 
COROLLARY 1. The diagonal components of induced mass are increased, if the 
fluid density is increased in any region (and left unchanged elsewhere). 
Also, much as before, one can prove 
COROLLARY 2. Let T,. || and || Ty be associated with solids, of which the 
second is obtained from the first by replacing a mass Am of fluid by solid. Then the 
(scalar) components of translation induced mass 7’, satisfy the inequality 
Ti, 2 Tin — Am. (4) 
Again, letting p, denote the initial pressure required to produce the initial accelera- 


° h 
tion a’, we have 


—Ty. = | (Vp,):a dk = | Div (p,a’) dR, 
©R “R 
since pa" = Vp, by the equations of motion, and 
Div (p,a‘) = p,-Div a‘ + Vp,-a* = Vpn-a', 
since the flow is incompressible (Div a° = 0). 
By the Divergence Theorem, 
—T', - | pa, ds, (5) 
’W 


where a, is the normal component of a’, and hence (by continuity) the k"” direction 
, and a similar moment producing factor for rotation if k = 1, 
is the total thrust in the k-direction, produced by the initial 


is the moment about the appropriate axis, 


cosine (for translation) 
2, 3. Thus, in this case, 7; 

h-acceleration. Similarly, if k = 4, 5, 6, T 
produced by the h-acceleration. We conclude 


THEOREM 2. The tensor component 7, represents the total h-component of 
initial pressure force required to produce a unit k-acceleration of the missile. 

The arguments given before’ apply without change to prove also, if gravity is negli- 
gible, and if an “‘h-curve” is defined as before’ to be a cylinder of finite cross-section 
bounded by parallels to the h-axis, in case h = 1, 2, 3, and as a circle perpendicular to, 
and centered on, the (h — 3) = axis if h = 4, 5, 6, 


THEOREM 3. Let a rigid solid be given a unit k-acceleration from rest, in an im- 
compressible liquid of variable density. Then the initial rate of increase in the h-compo- 
nent of liquid momentum is 7, , in any region C bounded by h-curves which contains 
the solid. 

It is worth noting that, if the density p(A) is constant on each sphere of a concentric 
family, or on each ellipsoid of a confocal family, then the acceleration fields can be ex- 
pressed as solutions of an ordinary differential equation. The details will be published 


elsewhere. 
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A NOTE ON THE APPLICATION OF SCHWINGER’S VARIATIONAL 
PRINCIPLE TO DIRAC’S EQUATION OF THE ELECTRON* 


By H. E. MOSES (nstitute for Mathematics and Mechanics, New York University) 





Schwinger’s variational principle has been used for a wide variety of problems in- 
volving wave motion in which it is desired to find the amplitude of a scattered wave in 
terms of the incident wave.** Schwinger’s method makes use of the fact that the ampli- 
tude of the scattered wave satisfies a variational principle. We shall indicate this 
variational principle briefly. 

Let us consider a vector space. We may define two different inner products, the 
Hermitian and symmetric inner product, in this vector space. The Hermitian inner 
product (a,b) of two vectors a and 6 is defined by the condition that 


(a. & = (6,.¢)* (1) 
where the asterisk indicates the complex conjugate. The symmetric inner product is 
defined by the condition 

(a, b) = (b, a). (2) 

Let us consider a vector space with a Hermitian inner product and consider a pair 


of equations 


j a= Ky, 
(3) 


te’ = K'y’, 
where AK and K‘ are Hermitian adjoint operators which by definition satisfy the con- 
dition 

(K’u, v) = (u, Kv) (4) 
for any two vectors u and v. If K’ = K, K isa Hermitian operator. It can be shown 


from (3) that 
(a’, y) = (y’, a) (5) 


which is called the reciprocity theorem. Let us define a number A by 
l l , 
A = * / = = 7 ees , (6) 


and the functional A{v,v’} by 


(v’, Kv K'v’, v) vs 
Af{v, v’} = _* ae... _\ ee (7) 
v,aj(a,v) (v , ay(a , v) 
so that 
Aly, y’} = 2. (8) 


*Received March 25, 1952. 
**See, for example, Schwinger’s unpublished nuclear physics notes, or the lectures of N. Marcuvitz 
in the notes “Recent Developments in the Theory of Wave Propagation”, Inst. for Math. and Mech., 
N.Y.U., 1949-50. The problem of the present paper is treated abstractly in the first set of notes. The 
point of view of the present note is close to that of the latter set of notes. 
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~ 


It can be shown that A{v,v’} is stationary for independent variations of v and v’ about 
the values y and »’ respectively and that, therefore, from (8), the stationary value of 
Afv,v’} is X. 

A similar statement holds if the symmetric inner product is used. In this case K’ 
is said to be the symmetric adjoint operator of K if K’ and K are related by (4). 

In Dirac’s theory of the electron, the elements of the vector space are functions 
f(x,y) of the coordinates denoted collectively by the vector x, and of a variable y which 
is restricted to four values which may be taken as 1, 2, 3, 4. These functions are called 
spinor components. The Hermitian inner product of f(z,y) and g(2,y) is given by 


> | ie, y)* g(x, Y) dx 


y=1°6 


and the symmetric inner product by 


d | s(x, vole, y) ae. 


y=l-¢ 
Dirae’s wave equation for the electron in an electromagnetic field is 


ee 2 = His, y; 0 (9) 
where H is an operator which operates on x and y and is given by 
H=H,+4q (10) 
with 
H, = > La = — mp 
- | (11) 
q = ey A ;(x) + eb(x) 


Here a; and 6 are Hermitian operators which operate with respect to the variable only. 
They satisfy the following commutation relations 


a;a + gE; = 26;;1, 
Ba; +a;8 = 0, ? (12) 
(s)° = I, 


where J is the identity operator. 

The operators a; , 8 can be expressed as integral operators with kernels a;(y,y’), 
B(y,y’) which are the well-known Dirac matrices. 

We have taken h/2x = c = 1. The mass of the electron is m and its charge is e. 
The functions A,(z) and ¢(z) are the vector and scalar potentials of an electromagnetic 
field and are taken as real and are assumed to vanish if | x | > r, , for some 7. 

We shall look for solutions of equation (9) which can be written as 


W(x, ¥; 0) = e "x(x, 1; BE) (13) 
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so that equation (9) leads to 

H x(x, y; FE) = Ex(a, 7; E). (14) 

We shall write 
X = Xin T Ase » (15) 

and require that x;, be a solution of 

Hox.,.(2, y¥; E) = Ex,,(2, ¥; £) (16) 
where H, is given by (11). Suitable solutions are the “spinor plane wave’”’ solutions 
which have the form 

pitkl(nz) 


xin(, ¥3 E) = x75 B, 0, 1) 7% 


(on? (17) 


where 7 is a unit vector which specifies the direction of propagation, (nz) is the inner 
product of the vectors x and 7, and | k | is the absolute value of the momentum vector 
and is given by the relation 


kk)? = BE? — m’. (18) 


Here 7 is a variable which is restricted to two values which may be taken as +1 and 
—1. The significance of 7 is that it represents the component of the spin in the direc- 


tion of the momentum. 


By substituting (31) into the equation Hy xi, = Ex;, it is seen that the functions 
x(vy; E, , 7) x(y; E) satisfy the following equation 
Ba , 

s$E+ ik 2 aj.n; + ambxinsB, 7, 7) = O. (19) 


For the purpose of the present note, it is not necessary to give an explicit form for the 
functions x(y; £); these functions can be found in textbooks. However, it will be useful 
to indicate the orthogonality properties of the functions. These orthogonality relations 


9re 
Al 


> x(y; E, n, 7)*x(y; E, n, rT) 


= ) 
| (20) 


Zz x(v; —E, n, 7)*x(v; E, n, 7’) 0, \ 
also 


>. v(y; BE, n, 7x7’; EB, n, 7) + Z x(vy; —E, n, ax(v’; —E, 0, 7)* = (7, vy’). (21) 


From (14), (15) and (16) it is seen that x,, satisfies 


E — Ho)lx..(2, y¥; FE) = ox(2, ¥; 2) (22) 





the solution of which can be written in terms of influence function g(2,7; 2’, y’), 


(7, 7; £) = 7. [ g(a, ¥; 2’, vy’ )ax(2’, vy’; E) (23) 
97’ 4 
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where 
(E — Hy |g(x, y; 2’, vy’) = ix — 2’) by, vy’). (24) 
In (24) the operator [E — H,] operates on the variables x, y. Since the physically in- 
teresting problem is that for which x,, is an outgoing wave, we shall take an outgoing 
wave solution for (24). It will now be shown how this solution which we donote by g 
is obtained. 
From the commutation rules (12) for a; , 8 one has 


E+4H,) =(F —-m+V"°%) =(k/?+ VY’). (25) 





[E + A, |(E — H,] = [2 — A] 
Consider now a solution of the differential equation 
(EB? — m? + V2)s8(2, 2’) = (|k |? + V2)8(2, 2’) 


, 


= 6(x — 2’). (26) 


(The subscript x on the operator V~ indicates that the differentiations are to be carried 
out on the variable x rather than 2’.) 
Any solution s(z,2’) of equation (26) can be used to form a solution of equation (24). 


From (25) and (26) one has 





[(E — H,|[E + H,]s(x, x’) iy, y’) = G(x — 2x’) Oy, vy’). (27) 
Hence a solution g(x, y; x’, y’) of (24) is 
g(x,y; 2’, y’) = [E + A ls(a, x’) 57, vy’). (28) 


This method of obtaining influence functions for the Dirac operator is a well-known 
one. A solution s,(z,2’) of (26) which leads to an outgoing wave is 


= ns (29) 


( (| k v Fr | 
j= ‘ Z=- 1a + 1 
q.(x, Y; zx’. 7’) ee : i, >» a,(y, y')(x, ions x’) 
I—-2 rae 
, i ole , \ ‘ 
mB(y, 7’) + EB ay, vy’) — am. (30) 
4dn | x—2 
Here a,(y, y’), B(y, y’) are the matrices which represent the operators a, 8. The function 
if f 


. . 9 . iEt 2\)\ 
g, represents an outgoing wave because the time factor is ¢ (see (13)). 


For large values of | x |, the expression for x,. becomes 


sik . ( 2 
4 7] € _ , , , 
x.(2, y;£) = —- “ z. iE (yy, y’) —|k > aly, Ym: — mB(y, ¥ | 


1 


a 


. 1 
OAs, ¥ 5B) dz’ i. (31) 


In the above expressions and in those following the result of operating with an operator 


such as q on a function f(z, y) will be a function of x, y which will be denoted by qf (z, y). 
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Here n,; is a unit vector with components 7,,; defined by 
zr=|xrl|m (32) 


We should like to re-write (31) so that the amplitudes of the spherical waves are 
inner products in order that we may ultimately use the variational principle to obtain 
these amplitudes. This line of thought motivates our use of the identity 


E by, 7’) — |k| Daly, y’)m — mBQy, )| 
= 2E ¥ x7; E, 0, Dx’; E, 2, * — (33) 


Though this identity is fundamental in our treatment, we shall not prove it, since it 
follows directly from (19) and (21). 

If the incident wave part of the function x(z, y; Z) has the direction n’ and the value 
of r is r’, we shall denote x(z, y; E) by x(z, y; E, 7’, r’). We proceed to define the ‘spinor 
spherical wave” 6(2, 7; E, n, 7) by 

ieee ' 
Oz,7;E, 7,7) = - — x(7; E, 9, 7). (34) 
2x | 2 
The spinor spherical wave is analogous to the spinor plane wave (17). Furthermore, 


we define T(E; n, 7; ’, 7’) by 


T(E; , 73.7’, 7’) = [ x's, n, t)*e 2" ay(2’, vy’; E, 0’, 7’) da’ 


= | latx’'s B,D xe’, YE, a, rae’. (85) 


The expression for x,, takes the form for large | z |, 
x.(¢, 7; E) = > Oa, y;E,m, DT(Ej m , 737’, 7’). (36) 


The scattered wave may be regarded as the sum of two spinor spherical waves, each 
being characterized by a different value of r. The function T(E; m , 7: ; m2 , 72) may be 
regarded as the amplitude of the spinor spherical wave in the direction 7, and with 
7 = rt, when the incident spinor plane wave has as its direction of propagation , and 
its value of ris r>. We have therefore obtained x,. in a form where T(E; n, 7; 7’, 7’) is 
an inner product to which, as will now be shown, a variational principle can be applied. 

In order to show how the variational principle discussed abstractly earlier may be 
used to find the amplitudes 7(E; 7, r; n’, 7’) we need only show what quantities are to 
be identified with the quantities appearing in equation (3). 

Let us first consider T(E; n, 7; n’, 7’) as being the symmetric inner product of two 
vectors. We identify x(x, y; E, n’, r’) with the vector y of equation, (3). We see that 
the vector a’ is to be identified with g* x(y; E, n, r)* e"''"' (x, n). From (15) and (23) 
we construct the equation corresponding to the first of equations (3). 


gx(y; E, 9’, re" = ax(z, ¥;3 E, 1’, 1’) 


-¢2 / g(x, 7; 2’, v')ax(2’, v5 B, 0’, 7’) dx’. (37) 
a 
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NOTES 


The integral equation corresponding to the second equation (3) will also be given, 
. Accordingly, the integral 
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(38) 


T) 


since it is not difficult to obtain symmetric adjoint operator K’ 
equation for A which is the function corresponding to the vector y’ is 


os ae oS gt A(z; 7; E, n 


ie ie 


ath id pS g.(z, ¥; x’, vy’ )qt A(z’, y’, E, 0, 7) dx’ 


g*x’(7; E, 9, 
(39) 


7’ 


vy’) = g(z’, v'3 2, ¥)- 


where 
g(x, ¥; x’, 4 
Just as the integral equation (37) for x(x, y; E, n, 7) was derived from the differential 
(40) 


tr) = Ex(z, y; E, n, 7) 
A(z, ¥; 


equation 
Ax(z, vy; E, n, 
together with appropriate boundary conditions, the integral equation (38) for A(a 
(41) 


E; n, 7) can be derived from a differential equation with certain boundary conditions. 


The differential equation satisfied by A is 
H* A(x, y; E, n, 7) = EA(a, y; E, 2, 7) 
so that A is an eigenfunction of H*. The boundary condition on A is that it is to be 
expressed as the sum of an incident wave A,, and a scattered wave A,, such that Aj, 
satisfies 
H*A,,(z, y; EZ, 0, rt) = EA,, (2, ¥; E, 9, 7) (42) 
| x x | for large values of | x | . 
* where x;, 1s given by (17.) 


and A,. behaves like the spherical wave (¢ 
We take as a suitable solution A;, of (42) the function x 
The function A,, as can be shown from (40), satisfies the equation 
; E) dz’ (43) 


, 


— ’\ 7~* er 
YE VIC AS TS 


Ais 
where the inverse operator [fF — H*%]~’ is represented by an integral operator with the 
—_ 
ee re) 


kernel g, . 


The functional corresponding to X(v’,v) whose extremal value is 1/[{7'(/ 
is written as 1/7 {v’,v} and is given by 


i dx 


Liv, 
he [ v(x, yiq(z, vy) — g b [ g(x,y; 2’, y’)qu(x’, vy’) dx’} 
ones ax)( > [ q*x(v; E, 9, 1% o(z, dz) 
} and X\L, Y¥; 


‘ 
Y 


(> v(x, yax(y; EB, n’, 7’) 
where v’ and v are the trial functions which approximate A(z, y; E, 7, 7) 
E, n’, 7’) respectively. 

Having considered the case where 7'(E; 7, 7; 7’, r’) is a symmetric inner product, we 
shall now discuss the case in which this amplitude is considered a Hermitian inner product 
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of the two vectors q x(7; EZ, , r) e’ | k | (n, x) and x(a, 7; E, n’, 7’). The identification 
of the vectors a and y and of the operator K are as before in the case of the symmetric 
inner product. The vector a’ is now identified with g x(y; E, n, re’ “'“"”. The 
vector y’ is identified with the function Q(2, y; EF, 7, 7) which satisfies the following 
integral equation which corresponds to the second of equations (3). 


ilki(n.z 


qx(y’; E, n, 7) *" = gQ(z, 7; E, n, 7) 


—-q>, [ g(x, ¥;3 2’, y’)qQ(2’, vy’; E, n, 7) dx’ (45) 
y’ © 


where 
ss : k z-—-7’|-ie< : ‘ 
g(x, 3 2’, v7’) = g(a’, y'5 2, 1)* = << a $4—9 dD, aly, ya: — 2%) 
|\z—- 2 a= 
— mB(y, 7’) + E by, Yt kmlenn 


It can be shown that Q(2, y; EF, 7, 7) as a solution of a differential equation with 
suitable boundary conditions. The function Q(z, y; E, n, 7) satisfies the differential 


equation 
HQ(x, y; E, n, 7) = EQ(2, y; E, 0, 7) (47) 


so that Q, like x, is an eigenfunction of H. However, the function @ has different boundary 
conditions than x. The boundary condition on Q is that it should be expressed as the 
sum of an incident part Q,, and a “‘concentrating”’ part Q,., . The function Q,, satisfies 


H,Q,,(z, y; E, 0, 7) = EQ,,(z, y; E, n, 7) (48) 


and is taken to be x;, as given by (17). The function Q.,.,, is specified by the condition 
that Q..,, is to behave like an inwardly moving spherical wave for large values of | 2 | . 
Using the boundary condition that Q,,,, represent an inward moving spherical wave for 
large values of | z | , the inverse operator [EZ — H,]~' can be expressed as an integral 
operator while the function g, is the solution of equation (24) expressed in the form (28) 
when the solution s of (26) is taken to be 


instead of s, given by (29). 
In the case of the Hermitian inner product the functional 1/7(v’,v) is given by 


‘ , ‘ 49 
> | v'(x, y)*{q(a, y) — ¢ » [ g(a, ¥3 2’, y')qu(2’, vy’) dx’} dx ad 








(x, y)*ax(y; E, 0’, r’)e* “9 da) [ q*x(y; E, n, r)*e7''*'*"" v(a, y) az) 
ey 
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where v’ and v approximate 
Q(z, y; LE, n, 7) and x(z, y; E£, 7’, 7’), 

respectively. 
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A FORM OF NEWTON’S METHOD WITH CUBIC CONVERGENCE* 


By W. M. STONE (Boeing Airplane Co. and Oregon State College) 


For obtaining an approximation to a root of a transcendental equation f(x) = 0 


Newton’s method may well be unsatisfactory because it is only quadratically convergent, 
thus requiring considerable interpolation if the functions involved are scantily tabulated. 
On the other hand, the formulas of Stewart [1], Hamilton [2], Bodewig [3], and others, 
which offer cubic or higher convergence have the serious drawback of requiring the 
evaluation of second or higher derivatives. Formula (4) below, based on the generalized 
Taylor expansion of Hummel and Seebeck [4], offers cubic convergence in terms of f(z) 


and f’(x) evaluated at points on each side of the root. 


Taking n = m in the Hummel-Seebeck expansion, 
Tr ; f’(a) + f(z) , f(a) — f''(@ 
f(z) = f(a) + (@ — a)* Sa + (x — a)" * 2 aoe, ee (1) 


we obtain two approximations to a root, 


—2f(a) —2f(b) 
M 4 — - aa . aa ‘ 1 a = I => — ‘ a 2) 
Z a '@+f@ un x ) f(b) + f(x) ( 


We choose a and b so f(a) and f(b) are opposite in sign, f’(a) and f’(b) same sign. Elimi- 
nation of f’(x) in equations (2) yields 
f(b) f(a) f(b) — f’(a) 0 (3 
= + ———s = (3) 
x—b I— a 2 


or, by an obvious procedure, 


b+a f(b) — f(a) J 
t= ; — : =» ) a 77 7 ~~ 7 ( 
2 2 f(b) -f@ f(b) — f (@) 


“~ 


(To-a, fi)—f@ 2 2b—a)f(d)\” 
7 -5 ep, CAD 


f'(b) — f(a) 


where choice of the ambiguous sign is quite obvious. 


*Received April 22, 1952. 
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Proof of cubic convergence of the method closely follows the more general dis- 
cussion of Bodewig. We take A as the true value of the first order root, 


f(z) = (@ — A)g(z), = g A) 0, (5) 
and expand g(x) in a generalized Taylor series, powers of (x — A). Setting r = a, , b, 


in the series of equations (2) and carrying out the indicated division we obtain 


Tusi— Mn _ stg gy HA) + Wine) , py, _ , | 
=. [ (Tat A) 29(A) + I i(a A) + ’ 


(6) 


Tnsi — bn ty ee gy) BA) + Gus) oo 
b,, ae A >> E (Tn41 A) 29(A) + Pa A) + ’ 
or, finally, 
na ~ de SAKA, H. (7) 
b, — a, 
where P(x — A) represents a power series in (a, — A), (b, — A), (%,4,; — A), quadratic 


terms and higher. 

As a numerical example consider the first root of x tan x — 1 = 0. Taking a = 0.8, 
b = 0.9 equation (4) yields the tabulated value 0.8603. Two or more applications of 
Newton’s method will involve interpolation if one has only a two place table at hand. 
Estimates of the magnitude of error involved in interpolation by means of (4) have been 


found by Hummel and Seebeck [5]. 
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THE SECOND FUNDAMENTAL THEOREM OF ELECTRICAL NETWORKS* 
By CHARLES SALTZER (Case Institute of Technology) 


1. Introduction. This paper will deal with an extension of the work of W. H. Ingram 
and C. M. Cramlet' as discussed by J. L. Synge.’ In addition it will be shown how their 
theories fit into a unified theory. The terminology of Synge’s paper will be used. 

A network may be represented by its Thévenin representation, i.e. by regarding 
its branches as consisting of impedances in series with constant voltage sources; or, it 
may be represented by its Norton representation, i.e. by regarding the branches as 
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consisting of a constant current generator in parallel with an admittance. In the first 
representation the Kirchhoff nodal constraints become homogeneous in the branch currents 
and permit the introduction of mesh currents and the mesh equations, whereas in the 
second, the Kirchhoff mesh constraints become homogeneous in the branch voltages 
and allow the introduction of what P. Le Corbeiller*® calls the basic voltages and the 
nodal equations. The theorem of Ingram and Cramlet which Synge describes as the 
fundamental theorem of electrical network theory is the fundamental theorem for the 
Thévenin representation of a network; the dual Norton theorem which is fully on a par 
with this theorem will be proved independently in this paper, and the relations of the 
two theorems will be discussed. In addition the nodal equations of G. Kron as given 
by P. Le Corbeiller will be deduced by a method analogous to that of Synge’s. It may 
be noted that the nodal theorem is much easier to prove than the mesh theorem. 

In section (2) we will prove the fundamental theorem for nodal networks and in 
section (3) we deduce Kron’s equations. In section (4) we prove that the nodal equa- 
tions have a unique solution; and in section (5) we show the relation of the theorems 
for the Thévenin and Norton representations. In section (2), (3) and (4) we shall consider 
only the Norton representation of a network. 

2. Basic node voltages for connected networks. Let N be the number of nodes: 
J, ,J2, +--+ Jy;and let B be the number of branches. In addition we assign a direction 
to each branch and define F, as the voltage rise from the initial point to the terminal 
point of the k** branch. A path on the network is defined as a finite sequence of branches 
such that each branch has one node in common with the preceding branch and the 
other node in common with the succeeding branch. The first branch is required to have 
only one node in common with the second and the last branch is required to have only 
one node in common with its predecessor. For two nodes J, and J, and a path directed 
from J, to J, we define V,, as the sum of the voltage rises of those branches of the path 
which have the same sense as the path minus the voltage drops of those branches of 
the path which have a sense opposite to the sense of the path. We note that for the 
given path Ve. =- V. , 

If we consider two paths joining J, and J, , these two paths taken in opposite 
senses form a closed path or mesh. Since there are no voltage sources in the network 
the sum of the voltage rises around the closed path is zero and thus Fas for one path 
plus V,, for the other path is zero. It follows from this that V,, is independent of the 
path. If we select one node say J, as a ground node and define V, = V,, (p = 2,3, --- n) 
then V, , the potential relative to the ground, is also independent of the path. We shall 
call V,, V;, --- , Vy the basic voltages. Since the voltage rise across a branch from J, 
to Ja is V, — V, we have: 

Theorem I(a) If the branch voltages of a network satisfy Kirchhoff’s mesh law then 
es such that the branch voltages can be expressed in terms 
of them. These basic voltages are the potentials of the nodes relative to the ground node. 

In addition if we prescribe N — 1 arbitrary basic voltages and define the voltage 
drops across the branches as above then we have: 


there are N — 1 basic voltag 


Therorem I(b) An arbitrary set of N — 1 basic voltages for a network generates a 
unique set of branch voltage rises which satisfy Kirchhoff’s mesh law. 
If E is the column matrix of branch voltage rises E, , FE, , --- Eg then 


i. TAY. Gui, (2.1) 
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where A,, is defined as follows: 


1 if J, is the terminal node of the p* branch 
A,, = ( —1 if J, is the initial node of the p® branch (2.2) 


0 if J, is not on the pt branch. 
In matrix notation eq. (2.1) can be written 
E = AV (2.3) 
Theorems J(a) and J(b) together may now be stated as: 
Theorem II: For any connected network taking any node as the ground node there 
is a matrix A as in (2.2) such that: 
(i) For arbitrary values of the associated basic voltages the branch voltages given 
by E = AV satisfy Kirchhoff’s mesh law on every mesh 
ii) For any set of branch voltages which satisfy Kirchhoff’s mesh law on every 
mesh there is a set of basic voltages V such that 


E = AV. 


This theorem is the dual analogue of the “central theorem” of Synge, and can be general- 
ized to apply to linearly independent combinations of nodes. In addition the restriction 
to networks which are not connected is readily eliminated by the use of ground nodes 
in each subnetwork. 

3. Node equations. If we designate the current source of the r‘” branch by J, we can 
define B branch currents U, , U, --- Us, each of which consists of the corresponding 
branch source current and the branch current due to the branch potentials. In matrix 
notation 


U =I1-YE (3.1) 


where Y is the admittance matrix of the network whose elements are the self and mutual 
admittances of the branches. 
We write Kirchhoff’s node law as 


> F.»U,=0 (¢=2,3,--: N) (3.2) 
1 


where 


1 if J, is the terminal node of branch p 
F., = \—1 if J, is the initial node of branch p (3.3) 


0 if J, is not on branch p. 
In matrix form this becomes 
FU = 0. (3.4) 


If we compare (2.2) and (3.3) we see that F = A, where A, is the transpose of A. From 
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(3.1) we have 


A,(I — YE) = 0. (3.5) 
If we define J’, the column matrix of basic currents by 
l'=A,lI, (3.6) 
we can write (3.5) as 
l'=A,YE, 
or by (2.3) 
I’ = A,YAV. (3.7) 


Further, if we define Y’ the basic node admittance matrix by 


then (3.7) can be written as 


Equations (2.3), (3.6), (3.8) and (3.9) are the basic equations of Kron’s nodal method 
together with (4.1) below. 
4. The existence and uniqueness of the solution of the nodal equations. We show 


here that if 7 = O implies HE = 0 then the matrix Y’ is non-singular, and hence the system 


of equations (3.9) has a unique solution. If J = 0 then by (3.6) I’ = O and (3.9) becomes 

Y’ V = 0. (4.1) 
Since V = 0 is a solution of (4.1) our assertion will be proved if we show that it is the 
only solution. But if J = 0 implies that EF = 0 then we must have V = 0. This is a conse- 
quence of the remark that if H = 0, the basic voltages of the nodes on any path are all 


equal, and since the network is connected and at least one node is connected to the ground 
node, the basic voltages of all the nodes are zero. Therefore Y’ is non-singular and from 
(3.9 we get 

Pail) ff. 


and by (2.1 


E = A(A,YA)" AI. $2) 


This is the fundamental formula for the solution of nodal networks by Kron’s method. 

5. Relation of the nodal theorem to the mesh theorem. For the given network let us 
consider any tree of the network. The voltage drops of the branches of the tree can be 
prescribed arbitrarily since there are no closed circuits on the tree to which Kirchhoff’s 
mesh law applies. If we choose any node as J, , construct V as above (see sec. 2.1), 
define the voltage drop across any two nodes as the difference between the values of V 
at the initial and terminal points of the node pair, and choose any sequence of nodes as 
defining a path, then Kirchhoff’s mesh law will be satisfied for any closed path. Thus 
if we restore the branches, called chord branches, which were deleted to form the tree, 
then Kirchhoff’s mesh law is certainly satisfied for meshes consisting of branch paths. 
Also if we replace these branches by their Thévenin equivalents then the currents in 
these branches which, as Synge has shown, are a complete set of mesh currents, are 
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determined by Ohm’s law. If we now replace the tree branches of the network by their 
Thévenin equivalents then by Synge’s theorem all the branch currents are determined 
and hence all the branch voltage drops are determined. Conversely if we assign the 
chord branch currents in the network regarded as consisting of Thévenin branches, 
then by Synge’s theorem, the branch currents and hence the branch voltages are de- 
termined. This implies that the branch voltage drops of the original network are de- 
termined. This completes the proof of the equivalence of the nodal and mesh methods. 

In addition, if the Thévenin equivalent branches are given then in Synge’s notation 


W=e-—Ziore=W+i (5.1) 


where e is the branch voltage-source matrix, Z the branch impedance matrix and 7 is 
the matrix of branch currents. In the notation of this paper if the Norton equivalent 
branches are given then by 


U=I1— YE. (5.2) 
Since we are dealing with two representations of the same network we may identify 
i and U, and W and E£. If Z is non-singular (5.1) can be written 
i=Z'e —Z'E, 
and by comparison with (5.2) we have 
I = Ze, YE = Z'E. (5.3) 
If we multiply (5.1) by C, on the left and note that e’ = C,e and i = Ci’ then we have 
e’ = C.E + C,ZC1’. 
Since e’ = C,ZC7i’ it follows that 
C.E = 0. 
But by eq. (2.3) this may be written 
C.AV = 0. 
Since V can be prescribed arbitrarily we have 
C,A = 0. (5.4) 


This last relation can be derived on purely topological grounds‘. Also since in (5.3) E is 
arbitrary we have 


Y=2Z". (5.5) 
Equations (5.3), (5.4) and (5.5) show the relations of the two methods. 
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NOTE ON THE ELASTIC DISTORTION OF A CYLINDRICAL. HOLE 
BY TANGENTIAL TRACTIONS ON THE INNER BOUNDARY. 


By Sistrr CHanpra Das (Chandernagare College, India) 


Introduction. Elastic distortion of a cylindrical hole by localized hydrostatic pressure 
has been discussed by H. M. Westergaard’ and C. J. Tranter.” 

In this paper a few problems of elastic distortion of a cylindrical hole by tangential 
traction on its inner boundary are discussed. The first case considered is concerned 
with the problem of an infinite elastic plate having finite thickness with a cylindrical 
hole acted on by localized tangential traction on the inner boundary, the two faces of 
the plate being free. In the second case one face of the plate is supposed to be fixed 
and the other free, while the hole is acted on by uniform tangential traction throughout. 
In the last case a localized tangential traction is supposed to act over a narrow band 
on the inner boundary of an infinitely long cylindrical hole in an infinite elastic solid. 

The solutions in the first two cases are obtained in terms of infinite series while in 
the last case it is expressed as an infinite integral. 


1. Method of solution. We take the axis of the cylindrical hole as the axis of z. 


Using cylindrical co-ordinates and assuming u = w = O and v to be independent of 6, 
we have 
Ol 
esx 2u ] = 
O02 
(1.1) 
1 oO 
Pu = (), Day = , ri 
/ o 
and 
o = ¢ =¢60,>T = {) 
| Ol fav 1 _ oO ‘y 
- = (, oa T.4 >= Gs nal } = t= (") (73) 
0. \or r/ or \r 
Two equations Ol equilibrium are identically satisfied and the third takes the form 
O ; 1 Oo v 07 
ea oe oe ee a se ee (1.3) 
or” yr or r Oz 
One partir ular solution suitable for the problem is 
As 
v= ; 1.4 
r 
Also substituting 
v = V cos kz or V sin kz (1.5) 
we get 
ME EA (1. nie 1.6 
—— oS : ie oe - (1.0) 
or ' 7 Or \ r) 


where V is a function of r only. 
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As a solution to the above we take 
V = ¢,K,(kr) (1.7) 
where A, is the modified Bessel function of the second kind of degree one. 

We now consider the following cases. 

2. The cylindrical hole in a large plate of finite thickness with a distribution of tan- 
gential traction over a narrow band. We take the plate to be of thickness 2Z and the 
surfaces defined by z = +L. Both ends of the plate are free and the tangential traction 
is localized within the zone z = +h. 

The boundary conditions to be satisfied are 


Ov 
| T77 =G — ae 0, when z=2+L (2.1) 
Oz 
and 2 te = —S,, when ei <a 
(2.2) 
= 0), when z|>h 
Assum< 
p= + Zz c.K (kr) cos kz (2.3) 
r 1 
Then from (1.2) we have 
2G Ay - : 
Tp = fe — .% okGK.(kr) cos kz (2.4) 
r k=1 
where A, is modified Bessel function of the second kind of the degree two. 
From the condition (2.1) we get 
wT 
| (2.5) 
L 
where 7 is any integer. 
Also the conditions given by (2.2) will be satisfied if 
{ = a’ S,h 
oe i oe 
and 
; 2LS, sin (nrh/L) = 
C, = 2 2A (2.6) 
nae GK,(nra/L) 
Hence 
S,a?h 4 = 2LS, sin (nth/L)K,(nar/L) cos (nrz/L) (2.7) 
= : : wT 4 — rf 
2GLr —_ na GK.(nra/L) 


which is evidently convergent as for large values of x 


K,(a) _ 
K(x) 
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3. The cylindrical hole with uniform tangential traction in a large plate one of 
whose faces is fixed and the other free. Here we take the plate to be of thickness L. The 
free surface is given by z = L, the face given by z = 0, being fixed. The uniform tan- 
gential traction is supposed to act throughout the hole. 

The boundary conditions to be satisfied are. 


1) v= 0, when z= § (3.1) 
, Ov P 
2) 7, = G-— =O, when z=, (3.2) 
Oz 
and 3) THe = —S, when 0<2z<L (3.3) 
at the cylindrical surface r = a. 
Assuming 
v= >> D,K,(kr) sin kz, (3.4) 


we get from (1.2) 


t15 = — », D,kGK,(kr) sin kz. (3.5) 


We find that the condition (3.1) is evidently satisfied and condition (3.2) will be satisfied 
if 

(Qn + 1)r = 

k= —_—" (3.6) 


where n is any integer. 
The condition (3.3) will be satisfied if 





*~ (2n + 1)’2’GK,((2n + 1)ra/2L) at 
Hence 
_ > 8SLK,((2n + 1)rr/2L) sin ((2n + 1)ra/2L) (3.8) 
a (2n + 1) GK,.((2n + 1)ra/2L) vies 


which is evidently convergent. 

4. The cylindrical hole in an infinite solid under a tangential traction over a narrow 
band. We now consider the case of an infinite solid having an infinitely long cylindrical 
hole acted on by a tangential traction which is operating over a narrow band of breadth 
2h. 

The boundary conditions here are 


Tre = —S, when lz|<h 


(4.1) 


= 0 when lz|>h 


at the surface of the cylinder r = a. 
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The conditions (4.1) can be expressed in the form 








—_ S i si iS 
(Te) rae oo 25) | = ht ms dt (4.2) 
T 0 
We assume 
y= | c(t)K,(tr) cos 2t dt (4.3) 


70 


as a solution of the equation of equilibrium (1.3) where c(¢) is a function of ¢ only. 
The boundary condition (4.2) will be satisfied if 
28, sin ht 


(i) = =. 4) 
" OrGK (ta) (4 


Therefore 
— [ 28, sin ht cos 2tK,(tr) dt (4.5) 
Jo UnGK,(ta) 
In conclusion I offer my sincere thanks to Dr. B. Sen for his help in the preparation 
of this paper. 


ON CERTAIN SOLUTIONS OF A PENDULUM-TYPE EQUATION* 
By GEORGE SEIFERT (University of Nebraska) 

Introduction. In the study of the oscillations of a synchronous motor around its 
average angular velocity, a differential equation of the following type, the so-called pen- 
dulum-type arises [1]: 

16 


+ f(9) = g(6) (1) 
adi 


d°6 
dt 


where {(@) and g(@) are functions of period 2 z in @. 

It has been shown [2] that in the case where f(@) = a > 0, a constant, and g(6) = 
8 — sin 6, where @ is a constant such that 0 < 8 < 1, there exists a constant a, = 
a.(8) > 0 such that if a < a, , eq. (1) will have a solution 6(¢) such that if @’(f) = 
y(@), then y(@) = y(@ + 2r) for all t, while if a > a, , no such solutions exist. Following 
Vlasov [3] and Minorsky [4], we call any such solution of eq. (1) a periodic solution of 
the second kind. Physically, such a solution corresponds to a subsynchronous level of 
performance of the motor described by eq. (1). It is known also [5] that questions of 
stability of solutions of eq. (1) with respect to the points of equilibrium of (1) involve 
questions of existence of such solutions. 

The purpose of this note is to exhibit a set of explicit conditions on f(6) and g(@) 
which insure the existence of periodic solutions of the second kind for (1). Since it has 
already been noted [5] that if f(@) > 0 and either g(@) < 0 or g(6) > 0 for all 6, there 
will exist such solutions for eq. (1), we restrict ourselves to the case where f(@) > 0 and 
the equation g(@) = 0 has a finite number of roots in 0 < @ < 27. 


*Received April 18, 1952. 








128 NOTES [Vol. XI, No. 1 


We also show (see Remark 1) that in the special case where f(@) = a > O and g(6) = 
8 — sin 6 where 0 < 8 < 1, these conditions lead to a somewhat larger lower bound on 
a, that was obtained by Tricomi [2]. 

Theorem. Let {(8) and g(@) be twice-differentiable functions of period 27 in 6, f(@) 
> 0. and the roots {@,},7 = 0, 1, --: , 2n, of the equation g(6) = 0 in the interval 
0. < 0 < 0% + 2m” = 6, be such that g’(@) > O and g’(@;) ¥ 0. Then eq. (1) will have 
a periodic solution of the second kind, provided each of the conditions below holds for 
each integer 7 such that 0 < 7 < 17: 


5 “+ ; ; fh ( G) ; ; ; 4 0 
(i) —f’(A)g’(@) < 2g (0)(: ‘Sia (6) ) for 0 = 62; ; 
~ ; (f(0 : F : 
(ii f(A GA < 2(g'(G2;) — gO) , TG )) for = Aaj-1 < 4 < 43; ; 
(iii) the equation F;(6) = f(@)h;(@) — g(@) = O has m > 2 tons tn, 4. =. 1, 2, 
-++, m, such that 6, <p? <p? < +++ <p)” < 6,;, for which F’(p;"’) ¥ 0 when 
i = 1, 2, and 
| q(@) d@ > | f(0)h;(0) dd + hj(p;"’) | f(0) dé (2) 
(i) for Bas < OS Br 
} ed = § 
(n!(0) for —;-2 < OS Oni, 
f { f?( \1/2 
(s 8) P 
} 4) = | | 9 (: 4 +g @.,)) | ds, 
an; (@) = | f S) ds a h;' (@ ly 7 


and 6% is the largest root of f(@)h; (p;) = g(@) such that 6% < PY 


. 
er 
i es 
= f 
* . 
== . t . 
~~ qiey/ 
% -/ 
R R, ¥ 
\ 
(ji 2d 
a ie \ an 
- AKT 
™—~_ >" 
~ 
/ oe te + 
yy 
—_— a +#* tio - a + 
P p* ¢ 8 
BOS; "| ; = 
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Remark 1. In the case for which g(@) = 6 — sin 0,0 < 8 < 1, and f(@é) = a>Q0O, 
we put 6; = arc sin 6B, 0 < 6, < 7/2, 6 = x — 6,, and & = 6, — 2x. Conditions (i) 
and (ii) of the theorem are trivially satisfied. If we require that 


a a” ‘ie B 
+ 4 + cos 6; < 78, (3) 


2 


then F(@) = ah,(@) — 8 + sin 6 = 0 will have two roots, p;"? < p;” for whichO < p;” < 
6, , and clearly h,(p;"’) < 8/a, which implies that 6% is such that 0, < 6% < —7. Hence 
inequality (2) of condition (iii) holds since 


B(0. — 6*) + cos 6 — cos 6% > B (6, — 6*) > hi(p,”)a(@. — 6%), 
Q 


and thus inequality (3), which can be written as follows: 


8 (B -1/2 
< 6, (4 + cos 0.) : 


a 


insures the existence of periodic solutions of the second kind in this case. The right hand 
member of this last inequality will, then, serve as a lower bound on a, defined in the 


introduction. Tricomi’s lower bound on a, is 


—1/2 
At: + cos 0.) [2]. 
WAT 

Proof of the theorem. If we put d@/dt = y, eq. (1) becomes dy/dt = g(6) — f() y, 
and the equation of the phase trajectories becomes 


dy _ g(0) — f(@y 
do y , i 


The singularities (in the sense of Poincaré) of eq. (4) consist of the points (0; , 0) of the 
(8, y) cartesian phase plane, the 6; being the zeros of g(@). An analysis of these singularities 
shows that the points (@; ,0) i = 0,1, --- , n, for which g’(6.;) > 0, are saddle points, 
and the phase trajectories on these points, the so-called separatrices, have slopes given by 


(42s) (Le y" 
9 +\ y+ 9b) 


respectively. 

We note that in the region R, of the phase plane bounded by the curve of the equa- 
tion y = g(0)/f(@) and the 6-axis all phase trajectories have positive slopes, while in 
remaining part of the phase plane, R, , the trajectories have non-positive slopes, having 
slope zero only at points on the curve of y = g(0)/f (8). 

For fixed j, let us consider the solution y;(@) of eq. (4) corresponding to the phase 
trajectory going into (@,; , 0) from the left with negative slope; i.e. y;(@) is such that 


2 


(0. *(O., ; 
lim y3(0) = awl D2) - (ac + 9'(0s))) 
Clearly, we have either y,;(0) > 0 for 6.;-2 < 6 < 6,; , or there exists a point (0;", 0) 
such that 0,,.. < 6: < 0; for which y;(0‘”) = 0. (See fig. 1). Suppose such a 0;" 
exists: then the solution y‘,(6) for which y;"}(6.;-2) = 0 and lims.e+,,-, y/i-1(8) > 0 
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must be such that y;";(@) > O for 6.;-2 < 6 < 62; ; for if not, its trajectory would have 
to intersect the @ axis in @.;_, < @ < 6; , which would imply the intersection of the 
trajectories oe to y;_;(@) and y;(6), which is impossible. Hence y;"}(@) > 0 for 
Ge;-2 < 8 < Or, 


If this holds: for 7 = 1, 2, «++ , n, it is easily seen that y,'’(@,) > 0, for if not, its 
trajectory would have to intersect that of some one of the y;'’(6), 7 = 1, 2, --+ , n, which 
is impossible. Hence, 0 = y5'?(0.) < yo (O2n) = Yo (4% + 2x). On the other hand, the 


solution y)(@) for which y;”’(@) = M, hee ‘re M is such that g(@)/f(@) < M for all 6, 
has the property that y5”’(@) > yo’ (@. + 2). By an argument used by Amerio [5], which 
we omit, we conclude that there exists for eq. (4) a solution such that y(@) = y(@ + 27) 
for all @. 

We need only show, then, that the existence of some integer 7, 0 < 7 < n, such that 
y;(0) > O for @.;-. < @ < 62; will contradict the conditions of the theorem. To this end, 
we show that y;(0) < h;(@) for p;” < @ < 6; , and y,;(0) < h;(p;”) for 6* < @ < p;” 
where p; and 6* are defined in the statement of the theorem. 

We first note that 6.;-. <p; < 62;-, . Define ¢;(6) = h,;(6) — y;(@); we have, clearly 
/(8)2,) = $;(62;) = 0 and will show that condition (i) implies ¢/’(6@,;) > 0. We have from 


eq. (4), 
—[y3(6)]° — f(0)y%(@) (0) ' 
y''(0.;) = lim { lytOT — fyi) + 9" es re) 
a3 y ;(@) : 
; f—2y%( oy?" 6) — f(@)yi"(6) — f'(Ay(@) + 2a) 
sin (oO ne 
0-+8 53 Y; ‘(0) 
— f'(62;) 
by L’Hospital’s rule. From this and the fact that sini = iis we obtain 
; “(6;) — dh; 
yl,’ (8o;) — g a vai} - ee b; , 


as “ ey 


and by a routine calculation which we omit we note that the assumption h{’(6.;) < 
y;'(8;) would contradict condition (i). Hence ¢/’(@.;) > O and there is an e > 0 such 
that for 6.; — «< @ < 6; , we have@;(@) > 0; i.e. h;(0) > y;(@). 

We show next that for any 6 such that @.,;-, < 0 < 6; , the assumption ¢;(6) = 0 
would contradict condition (ii). For suppose for 6 = @, in this interval, ¢;(6,) = 0; then 
there exists a §&, 0, < — < 6; , such that ¢/(&) = 0, ¢/’(é) < 0. From eq. (4) we have again 


_— — yal — f(@yi(@) + 9’(8) 
= y ;( 6) 


y;'(8) — f(9), (5) 


and also 
hha 
, y;(@) + f(@) 


hence eq. (5) becomes 
{= [yP — fy) + gO) + $8) _ 


f’(8). 
g( 6) ¢ 


y"(6) = 

If we put 6 = £ in this last equation, use h/(£) = y/(), and recall that g(¢) < 0, a routine 
calculation will again show that the condition h/’(é) < y/’(£) contradicts condition (ii). 
We next consider the interval p;”” < @ < 6,;_, . Recall that in this interval g(@) > 0, 
the equality being possible only for 6 = 6,;-, . Note also that y/(@) = —f(6) for 0 = @2;-_,. 
Now if for 6 = 6,, p\? < 0, < 62;-; , we have y;(0,) = h;(@,) = h;”(0,), there exists an 
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n such that 6, < » < 6.;-, for which y/(n) < —f(»); but this clearly implies 


un) = 7 0 
as yi(n) + f(n) 
which contradicts the assumption y,;(@) > 0 for @.;-2 < @ < 6; . We thus conclude that 
y;(0) < h,(@) for p;’ < 6 < 6a;-,. 
Finally, we note that for 6* < @ < p;”’, the curve of y = {(6)/g(@) is above the line 
y = h;°(p;”); this clearly implies that y;(0) < h;(p;”) for @ in this interval. 
To sum up, we have that for p;’ < @ < 6; we have y,;(@) < h,(@), while for 6% < 
0 < p;” we have y;(0) < h;(p;”’). 
We now substitute y = y;(@) in eq. (4), multiply by y,(6), and integrate from 06% to 
62; ; we obtain 


<0 


065i 65; 


—[y(@)F/2 = | g(@)ae— | $(0)y,(6) do. 
J 6,* / 6;* 
Hence 
sias jes 
| g(a) do < | f()y,() ae. 
| @e / 6; 
However, clearly, 
| | f(0)y;(@) dd < | f(0)h,(0) d@ + h,(p;”’) | f(@ dé 
«| 85° J pj(? - 6;* 


which, taken with the previous inequality, contradicts condition (iii). This proves the 
existence of 65", such that 6.;-. < 6;" < 62; , for which y,(@;") = 0; hence, the theorem. 
Remark 2. In the special case considered in Remark 1, we note that 


we) < 8 — sin 6 — ah,(@) 
ty\ } a seitin iat aa 


» 2 h,(6) for 0<d9< %. 


This shows that y = h,(@) is not a curve for which each trajectory crossing it from the 
right passes below it; if it were such a curve, conditions (i) and (ii) clearly would have 
been unnecessary in this case. 

Remark 3. Since it is usually a question of imposing conditions on the parameters in 
{(@), it is sometimes more convenient to replace inequality (2) of condition (iii) by the 
simpler but stronger condition: 

a Oa;5 


Oa; 
| g(0) d@ > h,(p;”’) [ f(@) dé. 


| O,5% 
On the other hand, condition (iii) can clearly, at the expense of simplicity, be weakened. 
We omit the details. 
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PULSED SURFACE HEATING OF A SEMI-INFINITE SOLID 


By J. C. JAEGER (Australian National University) 


1. Introduction. In many practical systems heat is supplied to the surface of a solid 
in regular rectangular pulses, so that the flux of heat F at the surface at time ¢ will be 
F=F,, nT <t<nT + T, 
eile (1) 
F = 0, nT +7, <t<(r+t+ DT, 


where F, , 7, and 7; are constants and the solid is supposed to be at zero temperature 
when ¢ = 0. Such problems occur, for example, in ‘‘on-off”’ heating, in heat generation 
by friction over portion of the surface of a rotating cylinder ’, in the rotating anode 
X-ray tube’, in the heating of machine guns, in the heating of the anode of a magnetron, 
and so on. 

It is usually desired to know the surface temperature after steady conditions have 
been attained, and, in particular, the temperature at the end of a heating interval. 
Fourier series are not very suitable for the treatment of such problems (since their 
convergence is very slow at the most interesting values of the time) and alternative 
methods have been given by Weber® and Oosterkamp’, but the most powerful and widely 
applicable method seems to be that of the Laplace transformation* in a form which is 
essentially equivalent to the steady-state operational calculus of Waidelich 

To illustrate the method the most important case, that of the semi-infinite solid 
heated over the whole of its surface, is discussed in some detail in §2: the results of this 
section are frequently useful as approximations for the case of a finite solid similarly 
heated. 

The more difficult problem of the semi-infinite solid with heat supplied over a circular 
area of its surface is of interest in connection with rotating anode X-ray tubes and is 
discussed in §3. Results for pulsed point and line sources in infinite medium are given 
in §4. 

2. The semi-infinite solid heated over the whole of its surface. 

The semi-infinite solid > 0, of conductit ity k and diffusivity a, rnetially at zero tem- 
perature, is heated over the whole of the surface x = 0 by the flux (1). It is required to find 
its surface temperature for large values of the time. 

If v is the temperature in the solid we have to solve 


Ov 1 ov 


— «= —— ow 0, a > 0, t> 0, (2) 
Ox a dtl 
with 
v = 0, when ¢ = 0, 2 > 0, (3) 
Ou : 
—-k— = F, x = 0, t> 0, (4) 
Ox 


*Received May 5, 1952. 
1Jaeger, Phil. Mag., (7), 35, 169 (1944). 
2Oosterkamp, Philips Res. Rep., 3, 161 (1948). 
Weber, Ann. der Phys., 146, 257 (1872). 
‘Carslaw and Jaeger, Operational Methods in Applied Mathematics, Ed. 2 (1948), §129. 
5Waidelich, Proc. Inst. Radio Engrs., 23, 78P (1946). 





1953] J. C. JAEGER 133 
and v finite as x —~, where F is given by (1). Writing 


* 


=f «va (5) 


“0 


v 


for the Laplace transform of v, with a similar notation for that of F, etc., we have from 
(1) 

Pod — ¢*"') > 
ee : (6) 


s7 


s(l — ¢ ) 


re = 


and, by the usual Laplace transformation procedure‘ 





F. 2 ~— -sT; rst/32q1/2 = 
y* = Poa me : ye : ’ (7) 
ks’"(1 — e") 
This gives for the surface temperature v, at z = 0, 
2F vl = ah 
— (24) ear<l,, (8) 
ke \wr, 
2F ya’? — on . 
i, = - oa (¢? —- @ — TI, lr, <t<T, (9) 
kr 


and so on. Also, the inversion theorem for the Laplace transformation gives 


~ QrikT . 


a alas {T — 7Te"*™" — T, + T,e~°" ae t -3/2 7. 


for all values of ¢, where y is a positive constant. The integrand of (10) has a branch 
point at s = 0 and simple poles at 


§ = +2n7i/T, ey ae 9, 2. ee. (11) 


When the integral (10) is evaluated by contour integration in the usual way, the 
residues at these poles give rise to a series of terms with period 7’, and the sum of this 
series, Which we shall denote by v, , is the required steady periodic term. The final result 


18 


or T ‘a OF 1/2 pe ol fy sn Pel _ 7 Tet 
DFT (at QF va ev = Te" —T+Te")du yy 


nO RT beg a i kT J, ia, u?(1 — e*“*”) , 


where the first term of (12) comes from the last term in the brace in (10), and the integral 
in (12) comes from the branch point in the first term in the brace in (10). 

For large values of the time the integral in (12) is negligible, and v» is given by the 
steady periodic part vp superposed on the first term of (12) which corresponds to steady 
heating by the average flux F)7',/T commencing at t = 0. 

For small values of ¢ the result (12) must agree with (8) for 0 < i < 7, , and with 
(9) for T, < t < T, and this gives an integral expression for vp over the whole of a period. 
Writing 


a=T,/T, Q = F.T, , (13) 
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the value of vp at time bT after the beginning of a period is 


: = 2 (5) a)b'” — r""I(a, b)} 0<b<a (14) 
3 ka “T) 5 Sy < 
o =e | | (s) (1 — ab’? — (b-— a)? — "71a, )}, + a<b<1, (15) 
( 
where 





[ eid — ale? — eo -* + a} at (16) 
ae ; 16 
0 (1 — e*) 


These integrals are easy to evaluate numerically. The values of a, the ratio of the 
heating time 7, to the period 7, vary greatly in practice. For on-off switching and 
frictional heating of a rotating cylinder they may be relatively large, but in other systems, 
such as machine guns or the anodes of X-ray tubes, they are usually of the order of 
0.02 or less. Values of vp at time of b7' after the beginning of a period are shown in Curves 
I and II of Fig. 1 for a = 0.5 and 0.1, respectively, for the same quantity of heat 
Q emitted per unit area per cycle. 

Curve III of Fig. 1 corresponds to the case a — 0, that is, to the quantity of heat 
Q per unit area being supplied instantaneously at the beginning of each period so that 
(1) is replaced by 


=Q 2d, (t — nT), (17) 


e, . . . 
where 6(é) is the Dirac delta function. 
In this case (14) and (15) are replaced by 


/ 1/2 A 9 p@ ‘helt —_ (£7 ye") ¢) © 
es 2 | ) oo SO os eee | : 523s ae (18) 
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Fic. 1. Periodic temperature oscillations of the surface of a semi-infinite solid when the same quantity 

of heat is liberated over half the cycle (Curve I), over 1/10 of the cycle (Curve II), and instantaneously 
(Curve ITI). 
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The most important quantities in practice are the temperatures at the beginning 
and end of a heating interval, obtained by taking b = 0 and b = a, respectively in (14). 
The variation of these with the parameter a = 7',/T for the same total quantity of heat 
Q supplied per unit area per cycle is shown in Fig. 2, Curves I and II. 
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Fic. 2. Periodic Temperature at the beginning (Curve I) and end (Curve II) of a heating interval. 
Finally it may be remarked that there is no difficulty in writing down expressions 


for the temperature at depth x. For example, for the surface flux (17) the temperature 
at depth z at a time ¢ after the beginning of heating which is so long that the starting 


transient has disappeared is 


2Q(at)'”” £ 
a lerf< Vat)” + up, (19) 
where the value of the periodic part at time b7' after the beginning of a heating interval is 
) 1/2 = side aie : 
Vp = Go 172 ¢ ee a le ierfe 5 _ 
k(rbT) kT 2b 





1/2 o —bg? es 2 —§* alk 
wv i [ u Al ss s z — di, (20) 
wkT /0 a6 _— ie ) 
where 
ierfe r = [ erfe & dé, 
c = 2(aT)~*, and a and Q are defined in (13). 
It follows from (19) that, for large values of the time ¢, the temperature is accurately 
the sum of the periodic oscillation vp whose mean value is zero and a term corresponding 
to heating for time ¢ by the constant flux Q/7, the mean value of the applied flux. The 
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same remark applies to all the other cases discussed. There seems to be no justification 
for the suggestion of Comenetz® that the temperature behaves on the average as if the 
constant flux Q/T had been supplied for time ¢ + 37. 

3. The semi-infinite solid heated over a circular area. 

Suppose the semi-infinite solid is heated by flux (1) applied over a circle of radius R 
in its surface, there being no loss of heat from the remainder of the surface. When steady 
conditions have been attained, the temperature at the centre of the circle at time bT 
after the beginning of a period is found to be 


2 
7 (21) 
kT T Up, «= 
where, now, the steady periodic part vp is given by 
Qa’? {2b ; cr/aby 4 © aug c 
Vp = ws(l —e ) + -—erlec a7 —C 
: kT” \an'” a 2b 
2 £Pr- er fe“ 8-H — oe ©1141 — cos ct} dt 
+=} ——— A=" "—=?, O<b<a, (22) 
7a . E"( l1—e*) 
where 
Q=FT,, a = T,/T, ec = R/(aT)"”. (23) 


As. before, the most interesting quantity is the maximum value of vp which is attained 
when } = a. Values of this, plotted against 7,/7,, are shown in Fig. 3 for various values 
of the parameter c. For c < 0.2 the integrals in (22) do not make an important contribu- 
tion, while for c = 1 they are already close to the limiting case c > © of Fig. 2, Curve II. 


Vp k T*/Q a’? 
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Fic. 3. Values of the periodic temperature at the end of a heating interval for heat supply over a circle 
of radius R. The numbers on the curves are the values of R(aT7’) *. 


‘Quart. App. Math., 5, 503 (1947). 
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4. Pulsed point and line sources in an infinite medium. 
For the point source, suppose that heat is supplied at the origin at the rate of q units 
per unit time in the time intervals 


nT <t<nT+T,, n = 0,1, 2, --> (24) 


no heat being supplied at other times. Then, when steady conditions have been attained, 
the temperature at distance r from the origin at a time bT after the beginning of a heating 
period is, for 0 <b < a, 








ga . 
i + vp, (25) 
darkr 
where the periodic part vp is 
( ie ee 2 f° ee OF — e*] sin CE 
lp = : Leet Ia at é as ~ | " dE ¢, (26) 
irkr 2b T Jo t(1—e*) 
s 
and 
a=T,/T, C =ral)™. (27) 


This is the limiting case of the problem of §3. 

For the line source, suppose that heat is supplied along an infinite line at the rate q 
units per unit length per unit time during the time intervals (24), no heat being supplied 
at other times. Then, for large values of the time the temperature at distance r from 
the line is, for 0 < b <a, 


_—£ F; “ta a 28 
dak Bi dat, ees aia 
where the periodic part vp is 
q(1 — a) C? r= oP CEL — ae — O79" 
oe (1 — a Bi - a | é J(Cé)[(1 — aje* | a Mate. S: (29) 
dork: 1h 2rk Jo #1 — e*) 


and a and C are defined in (27), Ei is the exponential integral, and J, the Bessel function 
of the first kind of order zero. 


BOOK REVIEWS 


Tensor analysis for physicists. By J. A. Schouten. The Clarendon Press, Oxford, 1951. 
x + 275 pp. 50s. 


The first five chapters of this extremely well written book are devoted to an exposition of tensor 
analysis and Riemannian geometry. The author uses his well known expository skill to give a clear 
geometric interpretation to the various quantities discussed. The last five chapters are devoted to appli- 
cations of this material to various physical problems. The two groups of chapters are separated by a 
summary of the salient points of the theory described in the first part of the book. This section should 
prove to be very useful. 

Chapters I and II deal with the algebra of affinors (tensors) and contain a detailed discussion of 
p-vectors and p-vector densities (anti-symmetric tensors and tensor densities with p indices). The 
identification of various affinors obtained by restricting the group of the space is made extremely clear. 
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This discussion is used in Chapter IV to make clear the meaning of Stokes’ Theorem which is the 
name he assigns to the theorem relating integrals over a g + 1 dimensional subspace of a space and 
integrals over the g dimensional] boundary of the space. 

Chapter V is devoted to a discussion of manifold which have an affine connection and contains in 
particular a discussion of those parts of Riemannian geometry needed in subsequent applications. 

Chapter VI is devoted to a discussion of dimensions of physical objects. The relation between 
dimensions of physical quantities and the choice of the underlying group is discussed. 

The applications of tensor calculus that the author discusses are: Elasticity and peizo-electricity 
(Chapter VII), Classical dynamics (Chapter VIII), Relativity (Chapter IX). Chapter X on Dirac’s 
matrix calculus gives a brief exposition of the calculus introduced into quantum mechanics by Dirac. 
Unfortunately the spinor calculus is not discussed. 

The space allotted by the author to these applications is relatively small. He has ‘“‘endeavoured 
to avoid an incoherent enumeration of interesting facts— and has tried to make each chapter a short 
but systematic introduction to some branch of theoretical physics.’ In the main he has succeeded 
admirably. However, the author devotes more than twice the space allotted to a discussion of graviation 
to a phenomenological discussion of hydrodynamics. The reviewer feels that the book would have been 
more useful if this allocation of space were reversed. 

In spite of this minor lapse, the book is to be highly recommended for use in a course on tensor 
analysis and as a reference book. 


A. H. Taus 


The rise of the new physics. By A. D’Abro. Dover Publications, Inc., New York, 1951. 
ix + 426 pp. (Volume I), 553 pp. (Volume IT). $8.00 the set. 


This is a corrected and enlarged edition (two volumes) of The Decline of Mechanism in Physics 
(1939); it is essentially unchanged except for the addition of 36 portraits of leading physicists and 
mathematicians. The book is an account of the development of classical physics and approximately 
the first decade of quantum mechanics. The aceount has a strong philosophical] flavor especially in the 
first volme. The book is certainly worthwhile for the student of physics and perhaps even more so for 


the mathematician interested in physics. 
Roun TRUELL 


Mémoire sur l’intégration graphique des équations aux dérivées partielles. Par J. Massau. 
Edition du Centenaire par les soins du Comité National de Mécanique, Palais des 
Académies, Bruxelles, Belgique, 1952. x + 544 pp. 200 Belg. frs. 


As a matter of general editorial policy, space in these columns is considered too valuable to be 
available for the review of reprint editions. In the present case, however, there is ample justification 
for deviating from this policy. 

In reprinting this treatise, the Belgian National Committee on Mechanics has rendered a significant 
service, not only to the memory of an outstanding Belgian scientist, but also to all those who are inter- 
ested in the application of the method of characteristics to the approximate integration of partial differ- 
ential equations of hyperbolic type. The volume contains a sequence of papers which originally appeared 
in the “Annales des Ingenieurs sortis des Ecoles de Gand” around the turn of the century, and which 
do not seem to be sufficiently well known among scientists working on problems of this kind. 

Pages 1 to 58 are devoted to the general theory of characteristics and pages 59 to 120 to the one- 
dimensional treatment of problems of unsteady flow in canals and rivers. By far the greatest portion of 
the volume (pages 121 to 391) is concerned with the study of limiting states of equilibrium in a co- 
hesionless soil. A 152-page note on the wave equation concludes the volume. 

There is no doubt that acquaintance with this material would have saved much time and effort to 
the scientists who, during the last twenty years, applied the method of characteristics to numerous 
problems in hydro-and aerodynamics, detonics, plasticity, soil mechanics, and other branches of me- 


chanics of continua. To mention only a few examples, in the section on soil mechanics the author 
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establishes the basic geometric properties of the net of failure lines (p. 259), recognizes the possibility of 
limiting lines (p. 163 ff.), and gives an exhaustive discussion of stress discontinuities (p. 267 ff.). In 
the theory of plasticity, the corresponding results have been developed ab initio whereas they could have 
been taken almost without modification from Massau’s work. The other sections of the book contain 
similar cases where Massau has anticipated results usually considered to be of much more recent origin. 


W. PRAGER 


Matter and motion. By J. Clerk Maxwell (Reprinted with notes and appendices by 
Sir Joseph Larmor). Dover Publications, Inc., New York, 1951. xii + 163 pp. $2.50 


clothbound ($1.25 paperbound). 


The first edition of this famous little book appeared in 1876 under the imprint of the Society for the 
Promotion of Christian Knowledge. An American edition, published by D. Van Nostrand Co., was 
brought out in 1875. The original edition was reprinted with notes and appendices by Sir Joseph Larmor 
in 1920 and it is this reprinting which has been reproduced in the edition under review. 

The book was originally intended by the author to provide in terms of the most elementary mathe- 
matics an exposition of the fundamentals of mechanics sufficient to serve as a foundation for the general 
study of physical science. As such it is still eminently readable for the elementary student. The notes 
added in the later edition by Larmor help to place the contents in proper perspective with respect to 
twentieth century revisions of classical physics, notably the theory of relativity. 

R. B. Linpsay 


Dialogues concerning two new sciences. By Galileo Galilei. Translated from the Italian 
and Latin into English by Henry Crew and Alfonso de Salvio. Dover Publications, 
Inc., New York, 1951. xxi + 300 pp. $1.50 (paperbound). 


The publishers are to be congratulated for making this celebrated book once more readily available 
in the fine translation of Crew and de Salvio, which appeared originally under the imprint of the Mac- 
millan Company in 1914 and was reprinted without change in 1939 by Northwestern University. 

The original edition of the Dialogues was published in Leiden in 1638, four years before the author’s 
death, at a time when his difficulties with the ecclesiastical authorities made publication in Italy im- 
possible. It is unnecessary to stress here the importance of the book in the history of physical science. 
In it the man, acknowledged by all to be the founder of modern physics, sets forth in simple terms the 
basic concepts of mechanics essentially in the physical form in which they have been used ever since. 
Every student of elementary physics should read the words of Galileo on the motion of falling bodies 
and on the mechanical properties of solids and fluids. The present inexpensive edition brings them 
within the reach of all who can read English. 

R. B. Lrnpsay 


Vektor-und Dyadenrechnung fiir Physiker und Techniker. By Erwin Lohr. Walter de 
Gruyter & Co., Berlin W. 35, 1950. xv + 488 pp. DM 24.00. 


This is the second edition of a book on vector and dyadic calculus which has apparently enjoyed some 
popularity in Europe since its publication in 1939. The notation is essentially that of Gibbs, and, as the 
title indicates, the book is of interest principally to physicists. No extensive alterations of or additions 
to the first edition are made. The changes consist largely in the addition of a final chapter of miscellaneous 
applications of vector and dyadic calculus. 

The book consists of three parts: I) The arithmetic and algebra of extensive quantities, II) The 
calculus of extensive quantities, and III) Physical applications. These topics, especially the first two, are 
covered extensively, but the author’s enthusiasm for his subject leads him to apply vector methods to 
some topics where little is to be gained and dyadics to some where the use of tensors is more satisfactory. 

The book will be useful to applied mathematicians and, in particular, to physicists. 

WituraM PELL 
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Differential equations. By Robert C. Yates. McGraw-Hill Book Company, Inc., New 
York, Toronto, London, 1952. vii + 215 pp. $3.75. 


This is a text book suitable for an introductory course for engineering students. Equations treated 
include the various standard types of first order, some special forms of first order, (mostly non-linear), and 
linear second order with constant coefficients. There is a brief discussion of solutions by numerical methods 
and by power series, the latter illustrated by the Legendre and Bessel equations. 

Partial differential equations are introduced by means of problems in vibrations and heat flow with 
solutions by separation of variables. 

Finally, there is a brief treatment of Fourier series. 

The book contains several review chapters inserted at appropriate stages and an abundance of prob- 
lems with answers. 

The illustrative applications of the various equations range over a great variety of physical and 
geometrical problems, emphasis being placed on first analyzing the problem, then formulating it in terms 
of a differential equation and finally solving it. 

The treatment is non-rigorous throughout and the writing is clear. 

GeorGcE W. Morcan 


The scientific papers of James Clerk Maxwell. Edited by W. D. Niven. Dover Publica- 
tions, Inc., New York. xxix + 801 pp. (two volumes bound as one). $10.00. 


This collection of papers is bound as a single volume and contains the papers formerly published in 
1890 in two volumes. This edition contains the published papers, lectures, and addresses of Maxwell and 
probably includes all of his work outside of his texts on ‘““Theory of Heat’’, “Electricity and Magnetism”, 
and ‘‘Matter and Motion”. 

While it may be superfluous to list the contents or to review something published in 1890, it is, I 
believe, not out of place to recommend as very interesting reading almost all of these papers by the first 
Cavendish Professor of Experimental Physics. 

20HN TRUELL 


Proceedings of the Second Canadian Mathematical Congress, Vancouver, 1949. University 
of Toronto Press, Toronto, Canada, 1951. xxxi + 255 pp. $6.00. 


This volume contains the papers presented at the Second Canadian Mathematical Congress, held at 
Vancouver, B. C., from September 5-9, 1949. The eighteen papers presented cover a wide variety of 
topics. Historical and survey papers are presented as well as technical work on pure and applied mathe- 
matics. 

The papers included are: A. Zygmund, Polish Mathematics between the two wars; P. A. M. Dirac, 
The relation of Classical to Quantum Mechanics; A. W. Conway, Hamilton, his life, work, and influence; 
H. J. Bhabha, Recent scientific developments in India; G. Szegi, Principal frequency, torsional rigidity and 
electrostatic capacity; L. Schwartz, Les Mathématiques en France pendant et apres la guerre; W. H. Watson, 
The role of the National Research Council in the education of Canadian mathematicians and physicists; 
B. Jessen, Mean motions and almost periodic functions; R. L. Jeffery, Non-absolutely convergent integrals; 
G. F. D. Duff, A developme nt in the theor y of the F-Equation; Josephine Mitchell, An example ofa complete 
orthonormal system and the kernel function in the geometry of matrices; P. Mand] and J. R. Pounder, Wind 
tunnel interference on rolling moment of a rotating wing; R. De Vogelaere, Une nouvelle famille d’orbites 
périodiques dans le problem de Stérmer: les ovales; B. Davison, 8. A. Kushneriuk, and W. P. Seidel, /nflu- 
ence of a small black cylinder upon the neutron density in an infinite non-capturing medium; W. P. Thomp- 
son, Thermal convection in a magnetic field; E. Leimanis, The application of infinitesimal transformations to 
the integration of diffe rential equations of exterior ballistics by quadratures; A. E. Scheidegger, On aravita- 
tional radiation; M. A. Melvin, Symmetry and affinity of electromagnetic fields, charges, and poles 


W. H. Pe. 
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Advanced engineering mathematics. By C. R. Wylie, Jr. McGraw-Hill Book Co., Inc., 
New York, Toronto, London, 1951. xiii + 640 pp. $7.50. 


This text book for third or fourth year engineering students treats material of proven value to the 
analytica] engineer. The main mathematical topics are: linear differential equations with constant co- 
efficients, Fourier series, the Laplace transformation, separable partial differential equations, Bessel 
functions, complex variables, conformal mapping, vector analysis, and numerical analysis. In the de- 
velopment of these topics the relation to physical problems is constantly brought out. Thus, analytic 
functions are related to fluid mechanics. A valuable feature is a chapter on electrical and mechanical 
vibrations which explains electrical circuits and how to set up electromechanical analogies. There are 
an extremely large number of diagrams and worked out examples; this should serve to make the book 
easy to read. 


R. J. Durrin 


Tensor analysis. By 1. 8. Sokolnikoff. John Wiley & Sons, Inc., New York and Chapman 
& Hall, Ltd., London, 1951. ix + 333 pp. $6.00. 


This text should furnish a good introduction to tensor analysis for students on the senior-graduate 

level. The exposition is clear and is illustrated by examples. The subject matter is arranged so that the 
formal processes of tensor analysis are developed and then applications are made to geometry, analyti- 
cal dynamics, relativity, and the mechanics of continuous media. In particular, the author’s develop- 
ment of the theory of non-linear elasticity should prove of interest. There are a few minor errors and 
misprints. 
The text covers « considerable amount of material which may be subdivided into: (1) matrix theory; 
the theory of tensor analysis; (3) applications of tensor analysis. In the first part, the author dis- 
cusses: linear vector spaces, linear transformations in terms of matrices, characteristic values and the 
reduction of two quadratic forms. The tensor theory is developed by introducing the group of coordinate 
transformations and studying the induced group of transformations for the covariant and contravariant 
components of tensors. This is followed by a study of tensor algebra, the metric tensor, the Christoffel 
symbols, covariant differentiation of tensors, the e-systems and generalized Kronecker deltas, and the 
Riemann-Christoffel tensor. In the applications of tensor analysis to geometry, the above results are 
interpreted in terms of curvilinear coordinate systems and base vectors in Euclidean three space. 
Parallelism of vectors, in three space, the Frenet formulas for curves, and an introduction to surface 
theory (imbedded in three space and intrinsic) are discussed. In the sections on analytical dynamics 
and relativity, the author considers the Lagrange equations of motion of a system of particles, the 
principle of least action, Hamilton’s equations, some elements of potential theory, the Lorentz-Einstein 
transformation law of restricted relativity, and the Schwarzchild line element. In the final chapter the 
mechanics of continuous media are studied by following the general approach of F. D. Murnaghan, 
Finite Deformations of an Elastic Solid, Am. Journal of Math. vol. 59, 1937. This leads to a unified and 
interesting development of the strain tensor in non-linear elasticity. 


(2) 


=] 


N. Cospurn 


Conformal mapping. By Zeev Nehari. First Edition. McGraw-Hill Book Company, Inc., 
New York, Toronto, London, 1952. viii + 396 pp. $7.50. 


Professor Nehari’s clearly written and thorough monograph on conformal mapping can be read 
with profit by all mathematicians who have an interest in this important subject. His very detailed 
study of special conformal transformations, and their analytical expressions, fills a real gap in the litera- 
ture. Especially to be commended are the elegant discussions of circular polygons and hypergeometric 
functions, distortion theorems, the symmetry principle, and of elliptic functions and the Picard theorem. 
The introductory discussion of multiply connected domains is also useful. 

The style of the book is modern, and all discussions very rigorous and complete. Because of this 
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fact, and the presence of ample exercises, it should serve as an excellent text for a course on conforma] 
mapping. Because of its specialized emphasis, it would seem more useful as a reference than as a basic 
text for a course on complex variable theory. 

Physicists and engineers who are interested in transformations of special domains of the classical 
type, should also find the book an excellent reference. (Other references for this purpose are A. Betz, 
“Konforme Abbildung” and H. Kober’s “Dictionary of Conformal Transformations,” published by 
the British Admiralty Computing Service.) Those who desire a sound basic understanding of the theory 
of conformal transformations, can hardly do better than study it carefully. 

On the other hand, there is little discussion of numerical methods applicable to general domains; 
reductions to tabulated real functions and problems of determining numerical parameters are not con- 
sidered; and the applications to non-Euclidean geometry and to problems of electromagnetism and 
fluid mechanics are hardly mentioned. For these reasons, and because of its general tone (viz., the proof 
on p. 278 of the elementary formula for tan (z +f), and the discussion of z = } (¢ + ¢~) as a special 
case of z = (a? +b¢+.0)/(de + ef +f) ), the book is less suitable as an “applied” mathematics text. 

However, “applied’’ mathematicians already familiar with practical aspects of the subject, and 
prospective mathematical analysts, will find the volume a nearly ideal reference, both for the general 
theory of conformal transformations, and for the treatment of special cases. 

GARRETT BIRKHOFF 


Foundations of high-speed aerodynamics. Facsimiles of nineteen fundamental studies as 
they were originally reported in the scientific journals. With a bibliography compiled 
by George F. Carrier, Professor of Engineering, Brown University. Dover Publica- 
tions, Inc., New York, 1951. 286 pp. $3.50. 


In addition to an extended bibliography, the book contains photo-offset reproductions of nineteen 
basic papers on the dynamics of compressible fluids. Obviously, space limitations have excluded some 
important contributions to the subject; otherwise, the collection illustrates well the rapid development 
of this field. The bibliography is arranged according to the following headings: Hodograph Method. The 
Rayleigh-Janzen Method. The Prandtl-Glauert Method. Supersonic Flow. Shock Waves. Boundary 


Layer. The Oscillating Airfoil and Other Unsteady Flow Phenomena. General. 
W. PraGer 


The theory of electromagnetic waves. (A Symposium). Interscience Publishers, Inc., New 

York, 1951. vii + 389 pp. $6.50. 

This book is a series of papers presented at a symposium at Washington Square College of New 
York University in June 1950. The symposium was sponsored by the Air Force Cambridge Research 
Laboratories and New York University. This book contains eighteen complete papers and three ab- 
stracts of papers as follows: On the Theory of Electromagnetic Wave Diffraction by an Aperture in an 
Infinite Plane Conducting Screen (37pp.), by H. Levine and J. Schwinger. On Systems of Linear Equa- 
tions in the Theory of Guided Waves (18 pp.), by W. Magnus and F. Oberhettinger. Wiener-Hopf 
Techniques and Mixed Boundary Value Problems (16 pp.), by 8. N. Karp. Asymptotic Solutions of a 
Differential Equation in the Theory of Microwave Propagation (12 pp.), by R. E. Langer. Criteria for 
Discrete Spectra (11 pp.), by K. O. Friedrichs. Extension of Weyl’s Intregral for Harmonic Spherical 
Waves to Arbitrary Wave Shapes (10 pp.), by H. Poritsky. Kirchhoff’s Formula, Its Vector Analogue, 
and Other Field Equivalence Theorems (17 pp.), by S. A. Schelkunoff. On the Diffraction Theory of 
Gaussian Optics (14 pp.), by H. Bremmer. Diffraction and Reflection of Pulses by Wedges and Corners 
(20 pp.), by J. B. Keller and A. Blank. Vector Wave Functions (10 pp.), by R. D. Spence and C. P. 
Wells. The W. K. B. Approximation as the First Term of a Geometric-Optica] Series (21 pp.), by H. 
Bremmer. Remarks Concerning Wave Propagation in Stratified Media (12 pp.), by 8. A. Schelkunoff. 
The Theory of Magneto Ionic Triple Splitting (32 pp.), by O. E. H. Rydbeck. An Asymptotic Solution 
of Maxwell’s Equations (38 pp.), by Morris Kline. Field Representations in Spherically Stratified 
Regions (53 pp.), by N. Marcuvitz. Propagation in a Non-homogeneous Atmosphere (34 pp.), by B. 
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Friedman. Reflection of Electromagnetic Waves from Slightly Rough Surfaces (28 pp.), by S. O. Rice. 
The Theory of Scattering of Radio Waves in the Troposphere and Ionosphere (abstract), by H. G. 
Booker. Properties of Guided Waves on Inhomogeneous Cylindrical Structures (abstract), by R. B. 
Adler. Evaluation of Integrals Associated with Wave Motion in Dispersive Media and the Formation 
of Transients (abstract), by M. Cerrillo. Electromagnetic Research in the U. 8S. Air Force Research 
Program, by N. C. Gerson. 

This book is obviously a discussion of mathematical methods applied to special problems in propaga- 
tion, diffraction, dispersion, etc., of electromagnetic waves. The book will undoubtedly be of value to 
the student of advanced electromagnetic theory and to the research worker. 

Roun TRUELL 


Advanced calculus. By Wilfred Kaplan. Addison-Wesley Press, Inc., Cambridge, 
Mass., 1952. xiii + 679 pp. $8.50. 


This text on advanced calculus is an excellent addition to the list of books on this subject, because 
of its laudable emphasis on vector methods in all phases of the calculus in which they can be applied. 

The material contained covers all the usual subjects found in a course of advanced calculus, but 
the emphasis seems definitely slanted toward the advanced engineer or the applied mathematician. 
The topics discussed include differential and integral calculus of functions of several variables, infinite 
series, Fourier series and orthogonal functions, functions of a complex variable, and ordinary and partial 
differential equations. In addition, there is an excellent treatment of vectors and their properties, vector 
differential calculus, and vector integral calculus, both in the plane and in space. The consistent use 
of vectors is well illustrated in the discussion of orthogonal functions, where the functions in question 
are thought of as elements in a vector space. Thus, the engineering student, to whom the concept of 
a general set of orthogonal functions and their properties may appear as an intellectual abstraction, 
will be able to relate these functions to what is to him a more or less familiar concept, vectors. 

There is a large number of problems of which a goodly portion is of the so-called applied type, 
i.e., relating to the fields of mathematical physics, hydromechanics, elasticity, potential theory, etc. 
The proofs of some of the theorems are given as problems. The rest of the problems, though not trivial, 
are of the standard drill type. Answers are given for all problems which are not self-contained. 

Some of the longer and more difficult results are given without proof, but whenever this is done, 
references are given where the proofs may be found. As the author states, ‘A competent teacher can 
easily fill in these gaps, if so desired, and thereby present a complete course in real analysis.”” Beyond 
this, the work is completely rigorous, clearly presented, and concisely written. 

At the end of each chapter there is a long list of suggested supplementary reading. 

Harry J. WEtss 


The principle of relativity. By H. A. Lorentz, A. Einstein, H. Minkowski and H. Weyl. 
Notes by A. Sommerfeld. Translated by W. Perrett and G. B. Jeffrey. Dover 
Publications, Inc., New York, 1951. viii + 216 pp. $3.50 (clothbound) and $1.50 
(paperbound). 


This collection of the classical papers by the above-named authors has been made available to the 
American market at a very attractive price and with the added convenience of being in English through- 
out. This edition will make it possible for graduate students so inclined to read the key papers that 
established relativity, and about as easily as they could get the same material out of a text book on the 
subject. One can only regret that this collection has not been supplemented by the most important 
papers that have appeared since, but such additions would undoubtedly have driven the price up to 
the point where purchase by students would no longer be feasible. 

It is interesting to note that the so-called Lorentz transformation is not due to Lorentz at all (I 
owe this remark to Dr. H. Zatzkis of the University of Connecticut) cf. page 14. Lorentz correctly 
describes the deformation of moving scales, but the transformation law for the time is quite different 
from the one postulated by Einstein (page 48). 

Einstein’s fundamental paper on the general theory of relativity occupies some fifty pages in this 
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small format. In a modern U.S. physics journal it would probably take about twenty pages. In that 


amount of space, Einstein manages to give a comprehensive introduction to the general theory of rela- 
tivity to the extent he had then been able to develop the theory, including an almost complete presenta- 
Very commonly the authors of today’s theoretical 


tion of the foundations of Riemannian geometry. 
papers excuse their lack of comprehensibility on the grounds that the editor will not let them use enough 


space. This paper by Einstein, covering the complete foundations of a new theory, presenting all of 
its physical and mathematical aspects, is “experimental” proof that within the length of a paner perfectly 
acceptable to today’s editors it is possible to write lucidly, provided the author is willing to go to the 
trouble of organizing his material with the reader in mind. It has been this reviewer’s privilege to see 
work on several papers intended for publication and to watch the infinite pains 


The 


Professor Einstein at 
with which he went over the drafts, sentence by sentence, until they were satisfactory to him. 
results are well worth the effort. 

The volume concludes with a paper by H. 
In this paper Wey! presented his generalization of Riemannian geometry, in which the length of a 
Occasionally, one reads in the literature references 


Weyl, which now has mostly historical significance. 


displacement vector is not an invariant quantity. 
to Weyl’s geometry which convey the impression that Weyl’s geometry reduces the length of a vector 
to a non-integrable status in the same way that Riemannian geometry makes the direction of a parallel 
displaced vector non-integrable. This impression is not quite correct, insofar as Weyl’s geometry 
permits the construction of vectors (those of weight-!4) whose norm is invariant. It is true, though, 
that the elementary displacement vector dz has no invariant norm. Unfortunately, in spite of its intrinsic 
conceptual beauty, Weyl’s geometry has not led to a viable physical theory. This may have been the 
apply it, including Wey] himself, and it is conceivable that another attempt 


fault of those trying to 
may succeed. 
Altogether, this collection of original articles may be warmly recommended to those who find the 


German original articles a bit too difficult and who, nevertheless, appreciate the value of studying the 
source ‘material. Undoubtedly, many readers will find it worthwhile to expand their reading beyond the 
papers available here. There are numerous references (such as to the original papers by Michelson and 
Morley, which were, of course, published in English), but no index. 

PeTteR G. BERGMANN 























